serological data on circulating influenza viruses generally contain evolutionarily important information about functional ( antigenic ) variation in the bcell antigens of those viruses . 

, 
  because natural selection acts on antigenic variation , serological data can provide important insights into patterns , causes , and epidemiological consequences of influenza viral evolution . 

, 
 , 
 , 
 , 
  nevertheless , there are important challenges to the use of serological data . in particular , serological assays require considerable amounts of time and effort to perform , and this limits the number of potential epidemic and pandemic viruses on which serological data can be routinely collected . 
  in addition ,
serological data are often contaminated by experimental noise ( e.g. resulting from the serial dilution of sera ) , which may cause independently measured data for the same virus and serum sample to vary greatly . furthermore , the data depend on nonantigenic variables such as the red cell avidity and the antibodyinducing capacity of viruses , 
 , 
  which can make it difficult to obtain accurate information about the antigenic variation of influenza viruses .
recently made great progress toward addressing the above challenges , focusing on data obtained from the widely used hemagglutinationinhibition assay .
those data are typically reported in a tabular format , with rows corresponding to viruses and columns to sera .
the entry found in the ith row and jth column of the table ( called the titer of virus i relative to serum j ) represents the reciprocal of the maximum dilution of serum j that can effectively neutralize virus i. smith et  al .
antigenic map of viruses and sera in which the euclidean distance between virus i and serum j was approximated by the negative log ( base 2 ) of the titer of virus i relative to serum j normalized by the maximum titer obtained for that serum .
( use of base 2 in the logtransformation of titers reflects the fact that titers are typically measured using twofold dilutions of sera .
the constructed antigenic map allowed unavailable normalized titers to be accurately predicted from distances measured on the map , increasing the amount of available data . because the map was constructed using data from multiple tables , data recovered from it are , in principle , less noisy than data obtained from individual tables .
here , i introduce computational methods that extend the above antigenicmap approach into novel directions .
first , i show that on average , logtransformed titers obtained in the same experiment ( i.e. reported in the same table ) have only one effective dimension .
then , i introduce a computationally efficient method that exploits this low dimensionality to accurately recover a complete table of titers using only a subset of those titers ; various quantities derived from titers ( e.g. the normalized titers of ref .
noisefree estimates and 95% confidence intervals ( cis ) can be computed for titers and , hence , for related quantities .
in addition , i derive a mathematical relation that sheds light on the mechanistic basis of titers and allows nonantigenic variation to be filtered out from these titers [ parts 1 and 2 of supporting information ( si ) ] .
i introduce a method for quantifying and visualizing antigenic variation using such filtered titers , and i illustrate this method using titers for novel h1n1 viruses .
the above methods were partly motivated by remarkable recent work on the recoverability of data that exist in a lowdimensional space . 


let   be the singular value decomposition ( svd ) of an m    n table h of logtransformed titers .
the columns ( called eigenvectors ) of u ( v ) are orthonormal bases for the column ( row ) space of h , whereas the diagonal entries of s are the n singular values of h , , sorted in decreasing order of magnitude .
the square of the ith singular value ( i.e. the ith eigenvalue of ) quantifies the variation explained by the ith eigenvector ; 
  therefore , the fraction of the variation explained by the first r eigenvectors is given by : ( 1 ) 
  the rank of h is defined as the number of its nonzero singular values . because some eigenvectors associated with nonzero singular values may explain negligible amounts of the variation in titers ( e.g. because they contain noise ) , the rank of h may be greater than its effective rank , defined here as the smallest value of r for which f 
 r 
 =  f 
 r+1 , r  =  1,  ,n1 ( the effective rank is set to n if this condition is not satisfied for any r  < 
n ) . alternatively , assuming that the noise found in logtransformed titers obtained in the same experiment is normally distributed , the eigenvectors containing noise would be associated with eigenvalues that are significantly smaller than the other eigenvalues . 
  in this case , the effective rank can be estimated by finding the smallest eigenvalue that is significantly greater than all the eigenvalues smaller than it , using fisher s variance ratio test .
this is equivalent to finding the maximum value of r for which 
 , 
 : ( 2 ) 
 where   denotes the inverse cumulative function of the fdistribution with 1 and nr degrees of freedom , and  is the level of statistical significance of the test . here , i used , where the factor ( nr ) represents a holmbonferroni correction for multiple comparisons . 
  both of the above methods for determining the effective rank gave similar results
let p 
 u ( p 
 v ) be the orthogonal projection of h onto the first r eigenvectors found in u ( v ) .
the coherence of h ( more precisely , the maximum coherence of u and v ) with respect to the standard bases in r 
  and r 
  is given by 
 : ( 3 ) 
  the lower the coherence , the greater the probability that measured titers , selected uniformly at random from h , will contain enough information about the unmeasured titers to enable their recovery by the method described elsewhere . let h denote an m   
n table of logtransformed titers of rank r , m    n , and let  denote a subset of rm 
 log ( m ) titers randomly selected from h with uniform probability .
( r  =  1 is used in this study ) . the unselected ( or  unmeasured ) titers are recovered by finding an m    n matrix x that minimizes 
 : ( 4 ) 
 where  and  are lagrange multipliers .  
denotes the nuclear norm of x , that is , the sum of the singular values of x. when r is known ( it is not approximated by the effective rank of h ) ,   is replaced by the sum of the r largest singular values of x. note that ( eqn . 
4 ) was solved using the fixed point continuation algorithm 
 ; for the experiments reported here , it took only a few seconds , on average , to recover titers found in each analyzed table using a pentium iv machine running windows xp . because titers recovered by the above method are theoretically exactly equal to their 
noisefree values when the selected titers are free of noise , 
  discrepancies between the recovered and the noisefree titers are necessarily due to noise found in the selected titers .
if there is no systematic bias in the way that noise found in the selected titers induces variation in the recovered titers , then the recovered titers would be randomly distributed about the corresponding noisefree titers .
the distribution of the recovered titers can therefore be used to compute cis for the noisefree titers .
this is the rationale for the following procedure for computing cis for titers found in h : 1 .
3 . repeat steps 1 and 2 until each titer found in h is recovered at least n ( = 1000 ) times ( on average it will take nn/[nm 
 log(m ) ] repetitions of steps 1 and 2 for this to happen ) . let l 
 i be a list of the k    n values recovered for the ith titer , which are sorted in increasing order .
then , the lower ( upper ) 95% ci for the ith titer is given by the 
th (  th ) element of l 
 i , where   ( ) denotes the largest ( smallest ) integer smaller ( larger ) than x. the mean of the recovered values for each titer is used as a noisefree estimate for that titer .
note that if h is incomplete , then 95% cis can still be computed for titers found in any m 
 1    n 
 1 complete subtable of h , provided m 
 1 
 log(m 
 1 )  <  m 
 1 
 n 
 1 . in practice ,
when designing an experiment to measure titers for m viruses relative to n sera , the experimentalist can measure only a random subset of at least m 
 log(m ) titers and subsequently recover the unmeasured titers as well as missing titers .
in addition , noisefree estimates and 95% cis for all titers ( respectively subsets of titers ) found in a complete ( respectively incomplete ) table can be computed and used in place of the original titers , which may be unreliable .
to quantify antigenic differences between viruses found in h , h is projected onto an rdimensional subspace of its row space : , where v 
 r denotes the first r eigenvectors found in v. the antigenic difference between viruses i and j is defined as the euclidean norm of the difference between the ith and jth rows of w. if r    3 , then antigenic differences can be visualized by plotting the rows of w. antigenic differences computed using titers from different tables can be embedded in a common rdimensional space by means of probabilistic multidimensional scaling ( part 2 of si ) .
note that the above svd decomposition of h is only feasible when no titers are missing from h. in practice , missing titers should be recovered before h is subjected to svd ( see part 2 of si for an svd approach that is applicable to tables with missing titers ) .
let   be the singular value decomposition ( svd ) of an m    n table h of logtransformed titers .
the columns ( called eigenvectors ) of u ( v ) are orthonormal bases for the column ( row ) space of h , whereas the diagonal entries of s are the n singular values of h , , sorted in decreasing order of magnitude .
the square of the ith singular value ( i.e. the ith eigenvalue of ) quantifies the variation explained by the ith eigenvector ; 
  therefore , the fraction of the variation explained by the first r eigenvectors is given by : ( 1 ) 
  the rank of h is defined as the number of its nonzero singular values . because some eigenvectors associated with nonzero singular values may explain negligible amounts of the variation in titers ( e.g. because they contain noise ) , the rank of h may be greater than its effective rank , defined here as the smallest value of r for which f 
 r 
 =  f 
 r+1 , r  =  1,  ,n1 ( the effective rank is set to n if this condition is not satisfied for any r  < 
n ) . alternatively , assuming that the noise found in logtransformed titers obtained in the same experiment is normally distributed , the eigenvectors containing noise would be associated with eigenvalues that are significantly smaller than the other eigenvalues . 
  in this case , the effective rank can be estimated by finding the smallest eigenvalue that is significantly greater than all the eigenvalues smaller than it , using fisher s variance ratio test .
this is equivalent to finding the maximum value of r for which 
 , 
 : ( 2 ) 
 where   denotes the inverse cumulative function of the fdistribution with 1 and nr degrees of freedom , and  is the level of statistical significance of the test . here , i used , where the factor ( nr ) represents a holmbonferroni correction for multiple comparisons . 
  both of the above methods for determining the effective rank gave similar results
let p 
 u ( p 
 v ) be the orthogonal projection of h onto the first r eigenvectors found in u ( v ) .
the coherence of h ( more precisely , the maximum coherence of u and v ) with respect to the standard bases in r 
  and r 
  is given by 
 : ( 3 ) 
  the lower the coherence , the greater the probability that measured titers , selected uniformly at random from h , will contain enough information about the unmeasured titers to enable their recovery by the method described elsewhere .
let h denote an m    n table of logtransformed titers of rank r , m    n , and let  denote a subset of rm 
 log ( m ) titers randomly selected from h with uniform probability .
( r  =  1 is used in this study ) . the unselected ( or  unmeasured ) titers are recovered by finding an m    n matrix x that minimizes 
 : ( 4 ) 
 where  and  are lagrange multipliers .  
denotes the nuclear norm of x , that is , the sum of the singular values of x. when r is known ( it is not approximated by the effective rank of h ) ,   is replaced by the sum of the r largest singular values of x. note that ( eqn . 
4 ) was solved using the fixed point continuation algorithm 
 ; for the experiments reported here , it took only a few seconds , on average , to recover titers found in each analyzed table using a pentium iv machine running windows xp . because titers recovered by the above method are theoretically exactly equal to their 
noisefree values when the selected titers are free of noise , 
  discrepancies between the recovered and the noisefree titers are necessarily due to noise found in the selected titers .
if there is no systematic bias in the way that noise found in the selected titers induces variation in the recovered titers , then the recovered titers would be randomly distributed about the corresponding noisefree titers .
the distribution of the recovered titers can therefore be used to compute cis for the noisefree titers .
this is the rationale for the following procedure for computing cis for titers found in h : 1 .
randomly select m 
 log ( m ) titers from h. 2 . recover the unselected titers ( see above ) .
3 . repeat steps 1 and 2 until each titer found in h is recovered at least n ( = 1000 ) times ( on average it will take nn/[nm 
 log(m ) ] repetitions of steps 1 and 2 for this to happen ) . let l 
 i be a list of the k    n values recovered for the ith titer , which are sorted in increasing order .
then , the lower ( upper ) 95% ci for the ith titer is given by the 
th (  th ) element of l 
 i , where   ( ) denotes the largest ( smallest ) integer smaller ( larger ) than x. the mean of the recovered values for each titer is used as a noisefree estimate for that titer . note that if h is incomplete , then 95% cis can still be computed for titers found in any m 
 1    n 
 1 complete subtable of h , provided m 
 1 
 log(m 
 1 )  <  m 
 1 
 n 
 1 . in practice ,
when designing an experiment to measure titers for m viruses relative to n sera , the experimentalist can measure only a random subset of at least m 
 log(m ) titers and subsequently recover the unmeasured titers as well as missing titers .
in addition , noisefree estimates and 95% cis for all titers ( respectively subsets of titers ) found in a complete ( respectively incomplete ) table can be computed and used in place of the original titers , which may be unreliable .
to quantify antigenic differences between viruses found in h , h is projected onto an rdimensional subspace of its row space : , where v 
 r denotes the first r eigenvectors found in v. the antigenic difference between viruses i and j is defined as the euclidean norm of the difference between the ith and jth rows of w. if r    3 , then antigenic differences can be visualized by plotting the rows of w. antigenic differences computed using titers from different tables can be embedded in a common rdimensional space by means of probabilistic multidimensional scaling ( part 2 of si ) .
note that the above svd decomposition of h is only feasible when no titers are missing from h. in practice , missing titers should be recovered before h is subjected to svd ( see part 2 of si for an svd approach that is applicable to tables with missing titers ) .
let h denote an m    n table of titers , with m    n. because it is tedious to measure all mn titers found in h , it will be very helpful if fewer titers can be measured and used to recover the unmeasured titers .
recent theoretical work 
  suggests that if both the rank r and the coherence of h ( see methods for definitions ) are not much greater than one , then unmeasured titers can be recovered exactly if only rm 
 log ( m ) titers ( selected randomly with uniform probability ) are measured .
i investigated the applicability of this method for recovering unmeasured titers by computing both the rank and the coherence of published tables of empirical titers 
  for influenza a and b viruses ( methods ) .
( because missing titers can artificially lower the rank of a table , only complete tables were analyzed ; see table  s1 ) . the analyzed tables , which typically consist of viruses sampled nonrandomly from multiple influenza seasons , should contain more variation on average than tables consisting of viruses of the same subtype sampled randomly from only one season
nevertheless , when the titers found in each table are logtransformed , on average 99% of the variation in those titers is explained by only one eigenvector , which is associated with the largest singular value of the table ( figure  1a ) .
i obtained a similar result by using a statistical ftest for the effective dimensionality of each table ( methods ) .
in addition , the tables have an average coherence of 136    012 .
note that deviations from an effective rank of 1 may be caused by noise and other factors , including the normalization of titers ( part 3 of si ) . also , based on wellknown results on the ranks of partitioned matrices , 
  when multiple tables of titers are combined the effective rank of the resulting composite table can be > 1 , especially if the rank of any one of the original tables is also > 1 ( part 3 of si ) .
the above factors may partly explain the higher dimensionality of normalized titers reported previously . 

 
( a ) the dimensionality of each of 23 tables of empirical titers ( table  s1 ) was investigated by logtransforming the titers found in each table and subsequently determining the fraction of the variation in titers ( denoted f 
 r ) that is explained by the eigenvectors associated with the r largest singular values of each table ( methods ) . f 
 r ( averaged over all 23 tables ) is plotted against r , for r  =  1,  ,5 .
( b ) a table consisting of simulated titers was contaminated by noise ( independently drawn from a normal distribution with mean 0 and standard deviation  ) , and the resulting noisy titers were used to compute estimates of and 95% cis for the uncontaminated ( or noisefree ) titers ( methods ) .
the fraction of noisefree titers that occurred within their corresponding cis is plotted against . ( c ) the mean ratio of the width of the 95% ci for a particular noisefree titer to the absolute value of that titer is plotted against . ( d ) the mean absolute difference between the estimated and noisefree titers , between the estimated and noisy titers , and between the noisy and noisefree titers , is plotted against . bars denote standard deviations .
the fact that both the effective rank and the coherence of the analyzed tables are not much greater than one suggests unmeasured titers can be recovered exactly by the abovementioned method .
note , however , that because measured titers are likely to contain noise , exact recovery of unmeasured titers may not be possible . also , recovery is not possible for rows / columns of h that do not contain any measured titer .
furthermore , accurate ( not necessarily exact ) recovery is only possible if the number of measured titers is not smaller than the number of degrees of freedom of h , which is given by r ( m  +  n    r ) , 
  a condition that is satisfied by all the tables analyzed earlier . in practice , each row / column must contain more than r titers in order for titers missing from that row / column to be accurately recovered .
it is important to keep in mind that exact recovery of unmeasured titers is still theoretically possible if the effective rank of h is > 1 , provided conditions ( i ) and ( ii ) above are satisfied .
in addition to unmeasured titers , the above conditions are also applicable to the recovery of 
titers  titers that are presumed to be either too high or too low to be measured experimentally .
( note that some missing titers may actually be within the limits of experimental resolution , but they could not be measured because of noise ) .
i applied the above method for recovering unmeasured titers to tables of empirical titers for influenza a and b viruses ( table  s1 ) .
the titers found in each table were logtransformed , and 90% of those titers were randomly selected and used to recover the unselected (  unmeasured ) titers ( methods ) .
[ note that it was necessary to limit this analysis to m    n tables for which 09mn    m 
 log ( m ) . ] both the mean absolute difference ( 066 
( 008    010 ) between the recovered and the unmeasured titers were small , suggesting that the recovery of unmeasured titers was accurate , albeit not exact . a closer look at the above results reveals that the mean absolute difference varies across the analyzed tables , ranging from 015 to 250 .
this may reflect the differing amounts of noise found in titers from different tables . to minimize noise found in empirical titers , i developed a novel method for computing noisefree estimates and 95% cis for titers ( methods ) .
the method was tested using 20    10 rank1 tables consisting of simulated , logtransformed titers .
each table was constructed as the product of two matrices ( of dimensions 20    1 and 1    10 , respectively ) , the entries of which were independently drawn from a truncated normal distribution with support ( 15 , 4 ) ; this support was chosen to simulate the dynamic range of logtransformed , empirical titers ( 332 , 1332 ) .
( note that because empirical titers are measured on a geometric scale , it is reasonable to approximate their distribution by a lognormal distribution . 
 ) to simulate noise found in empirical titers , for each constructed table m , another table m * was obtained by adding a perturbation ( independently drawn from a normal distribution with mean 0 and standard deviation  ) to each entry of m. m * was then used to compute estimates of and 95% cis for the noisefree titers in m ( methods ) , for   =  0 , 01 , 02 ,  , 1 .
[ note that when   =  1 there is a 32% chance that , on a geometric scale , a titer will be either multiplied or divided by at least 2 ( the dilution factor often used in serological assays ) , representing a high amount of noise . ]
representative results obtained using one of the constructed tables are shown in figure  1 .
the results show that the computed 95% cis have excellent coverage properties ; they contain their corresponding noisefree titers in > 99% of cases and they also have small relative widths ( figure  1b , c ) . remarkably , although estimates for the noisefree titers were computed using noisy titers , the mean absolute difference between those estimates and the noisefree titers is much smaller than between the estimated and noisy titers and between the noisy and noisefree titers ( figure  1d ) . indeed , the mean absolute difference between the estimated and noisefree titers grows much more slowly with  than does the mean absolute difference between the estimated and noisy titers ( figure  1d ) .
these results suggest that the developed method can , in principle , be used to systematically minimize noise found in tables of empirical titers . as an illustration ,
the average absolute difference between each titer and its noisefree estimate is small ( 066    067 ) .
however , 60% of the titers do not occur in their associated 95% ci , suggesting that those titers may be unreliable .
this is consistent with the high variation in independently measured normalized titers for the same virus and antiserum that is observed in many cases ( table  s1 ) .
in addition to noise , it is also important to minimize other types of nonantigenic variation found in titers .
this is currently not possible because of limited understanding of the fundamental nature of titers . to shed light on the nature of titers
, i developed and used a mechanistic model of hemagglutination inhibition to derive the first explicit mathematical equation for the titer of virus i relative to serum j [ eqn .
( s23 ) of si ] : ( 5 ) 
 where a 
  denotes the concentration of antibodies found in serum j , k 
  the average affinity of those antibodies for virus i , and j 
  a dimensionless quantity that depends on such nonantigenic variables as the avidity of virus i for red cell , the concentration of virus i , etc .
a 
  depends on nonantigenic variables , including the antibodyinducing capacity of the virus against which serum j was raised and the immune status of the organisms in which that serum was raised . 
 
the derived equation predicts that the normalized titer h 
 /h 
 , a commonly used measure of antigenic difference , 
  is approximately independent of a 
 , but it depends on both j 
  and j 
 .
in contrast , a measure of antigenic difference introduced by archetti and horsfall 
  [ h 

 h 
 /(h 

 h 
 ) ]  is predicted to be approximately independent of the nonantigenic variables a 
 , j 
 , and j 
 , suggesting that it may be more accurate .
importantly , the derived equation predicts that by meancentering each row ( resp . each row and column ) of a table of logtransformed , unnormalized ( resp .
i developed a method for quantifying and visualizing antigenic differences between viruses using such meancentered tables ( methods ) .
as an illustration , i used the method to construct a onedimensional map showing the serological relationships among 123 h1n1 viruses , 118 of which are novel h1n1 viruses ( figure  2 ) .
the map shows that the serological responses of novel h1n1 viruses are on average much more similar to each other than to the corresponding responses of earlier h1n1 viruses , particularly a virus isolated in 1930 and another isolated in 1976 ( figure  2 ) .
in addition , the map suggests the existence of a welldefined direction of antigenic evolution of the viruses . because
 serological maps , such as the one described here , have only one dimension , they may facilitate computational prediction of shortterm changes in the serological attributes of epidemic viruses . 
four tables of titers ( supplementary data  s1 ) that were recently collected by the us .
the tables had dimensions 26 by 10 , 36 by 10 , 40 by 12 , and 41 by 11 , respectively .
measured titers found in each table were logtransformed and used to recover the missing titers , except for titers missing from rows / columns in which fewer than two titers were measured .
each completed table was rowcentered and used to quantify differences between the serological responses of pairs of viruses ( methods ) .
these serological ( or  antigenic ) differences were subsequently mapped onto a onedimensional space by means of probabilistic multidimensional scaling ( part 2 of si ) .
the width of each displayed rectangle represents a unit of serological difference . the novel h1n1 vaccine virus ( ca09 )
the arrow indicates the apparent direction of temporal changes in the serological properties of the analyzed viruses .
virus labels : ca09 california 07/2009 , il07 illinois/9(33304)/2007 , nj76 new jersey/8/1976 , sw30 swine / iowa/1930 , tx08 texas/14/2008 , wi98 wisconsin/10/1998 . in summary ,
the computational methods presented in this paper suggest new possibilities for improving the use of hemagglutinationinhibition titers ( and serological data in general ) for influenza surveillance .
in particular , the method for computing estimates and 95% cis for noisefree titers would allow uncertainties associated with titers to be taken into account , for example , when selecting viruses for use in influenza vaccines . also , the method for minimizing nonantigenic variation found in titers may help to improve the estimation of antigenic differences between viruses .
in addition , the method for recovering unmeasured titers may help to reduce the experimental effort required to collect titers ; for example , only 20log(20 )  =  109 titers may need to be measured in order to accurately determine all 200 possible titers for 20 viruses relative to 10 sera .
this method can be applied to other types of serological data , especially when the number of viruses under consideration is greater than or equal to the fifth power of the dimensionality of the data .
the additional experimental capacity made available by the method would increase the number of viruses , circulating in humans and other organisms , for which serological data can be routinely collected and thereby improve the likelihood that potential epidemic and pandemic viruses will be identified before they become widespread . 
 
this is important in light of the recent pandemic spread of an influenza virus whose initial circulation in humans may have gone undetected for several months . 

 
note that the observation , reported here , that tables of empirical titers have an average effective rank of 1  titers exist in a space that has only one effective dimension  suggests that the effective number of independent variables that determine titers is very small .
in other words , while titers depend on many variables [ see equation  ( 5 ) and part 1 of si ] some of which may be mutually independent , variation in titers may be dominated by one or more codependent variables
. in so far as titers are determined by and contain information about viral phenotypes , the demonstrated low dimensionality of the space of titers and the consequent high recoverability of unknown titers suggests that it may be possible to predict biophysically accessible viral phenotypes from limited information about extant viral phenotypes .
welldesigned theoretical and experimental tests of this idea may yield important new insights and also shed light on questions concerning the evolutionary accessibility of epidemic variants of influenza viruses . 
 
in addition to influenza surveillance , the new methods presented in this paper can be used in the surveillance of other pathogens and in the investigation of basic questions concerning pathogen evolution and dynamics . as should be the case for existing methods ( e.g. recently developed methods for elucidating the antigenic structure 
 , 
  and the adaptation potential 
  of pathogen populations ) , continued empirical / experimental assessment of the new methods is necessary to ensure that they will remain useful .

 data s1 . list of tables of empirical titers for influenza a and b viruses . 
 table s1 . supporting text and list of tables of empirical titers .
 
 please cite this paper as : ndifon w. ( 2011 ) new methods for analyzing serological data with applications to influenza surveillance . 
 influenza and other respiratory viruses doi : 10.1111/j.17502659.2010.00192.x . 
 
 background   two important challenges to the use of serological assays for influenza surveillance include the substantial amount of experimental effort involved and the inherent noisiness of serological data . 
 
 results   i show that logtransformed serological data exist in an effectively onedimensional space . 
 i use this result , together with new mechanistic insights into serological assays , to develop computational methods for accurately and efficiently recovering unmeasured serological data from a sample of measured data , for systematically minimizing noise and other types of nonantigenic variation found in the data , and for quantifying and visualizing antigenic variation . 
 the methods can also be applied to data with effective dimensionality greater than one , under certain conditions . 
 
 conclusion   careful application of the methods developed here would enable the collection of betterquality serological data on a greater number of circulating influenza viruses than is currently possible and improve the ability to identify potential epidemic and pandemic viruses before they become widespread . although the focus here is on influenza surveillance , the described methods are more widely applicable . 
Introduction
Methods
Computing the effective rank and the coherence of a table of titers
Recovering titers
Quantifying and visualizing antigenic differences between viruses
Results and discussion
Supporting information