Abstract:
The disclosed invention provides a method for transforming lists to corresponding matrices and optimizing results of data analyses. Lists of enumerable signs (e.g., numbers, points, vertices, categories) that refer to real phenomena are transformed to manipulable matrices of numbers that are employed in data analyses. Manipulation of said matrices augments optimization of results of data analyses.

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
       [0001]    This application claims benefit, under 35 U.S.C. sctn. 119(e), of U.S. Provisional Patent Application No. 60/850,387, filed on Oct. 10, 2006, and U.S. Provisional Patent Application No. 60/850,424, filed on Oct. 10, 2006, which are hereby incorporated by reference in their entirety. 
     
     BACKGROUND OF THE INVENTION 
       [0002]    The present invention relates to data analyses, and, more specifically, to transforming lists to matrices and optimizing results in data analyses. 
         [0003]    Data analyses that employ mathematics and statistics in analysis of real phenomena have proven to be useful. For example, data analyses are useful in prediction and measurement of many things in many areas, including medicine, public health, economy, industry, education, and much more. 
         [0004]    In this invention, a list is defined as enumerable signs that refer to real phenomena. The signs in this definition could take the form of numbers, points, vertices, categories, or other phenomena that refer to real phenomena. In a list, numbers could be exemplified by weights or counts of diverse objects or individuals; points could be exemplified by points on a map of an area; vertices could be exemplified in a graph of a process, structure, or event; and categories could be exemplified by colors of objects or nominal characteristics of individuals. The present invention is restricted to enumerable signs, meaning that the signs are capable of being put into one-to-one correspondence with the system of positive integers, as exemplified by a list that is denoted by x where x=(x 1 , . . . , x n ) for n signs x i  where i=1:n, or as exemplified by a list with a first number, a second number, and so on, or by a list with point 1, point 2, and so on. The mathematical and statistical sense of enumerable lists is that each list of such lists is one variable; that is, each list is univariate (i.e., the list is where the signs could vary, at least in terms of their correspondence with the system of positive integers). 
         [0005]    Multivariate analysis is a data analysis of more than one variable. Therefore, multivariate analyses analyze multiple lists; they do not analyze a single list. Factor analysis, principal components analysis, cluster analysis, and regression analysis are examples of multivariate analyses of multiple lists of numbers (Everitt 1993; Hartigan 1975; Hotelling 1933; Pearson 1901; Rencher 2002; Spearman 1904). 
         [0006]    The present invention innovates by introducing multivariate analyses of a univariate list. 
         [0007]    The present invention, then, provides a method for transforming one list (i.e., a list of any kind of signs) to multiple (i.e., more than one) lists of numbers, and for conducting data analyses with said new multiple lists of numbers. 
         [0008]    In the present invention, a list x where x=(x 1 , . . . , x n ) for n signs x i  where i=1:n, is transformed to a matrix of numbers r jk , where j=1:n and k=1:n, where each number r jk  denotes a relationship between x j  and x k  of said list x, where said matrix is designated as matrix R, and where matrix R contains n 2  numbers r jk . The present invention provides a method for data analyses of a matrix R, that corresponds to a list x, or for data analyses of matrices R 1 , . . . , R q , that correspond to lists x 1 , . . . , x q  where q designates the number of matrices or lists. 
         [0009]    Transformations of univariate lists to matrices are common, and are found in many areas. For example, distance d jk  between numbers x j  and x k  can be computed with d jk =|x j −x k |, and matrix D containing n 2  cells containing distances d jk  can be constructed where j=1:n and k=1:n. Similar matrices are commonly constructed for distances between points in or on a geometric structure or distance between points in a specific space (e.g., n-dimensional Euclidean space) or distance between two vertices of a graph. The geodesic distance between two points is a locally length-minimizing curve (i.e., it is the shortest distance between the two points, and the shortest path that a particle which is not accelerating would follow between the two points). For example, in the plane, the geodesic distance between two points is a straight line, and the geodesic distance between two points on a sphere is the segment that connects them in the great circle of this sphere (e.g., the equator of a spherical planet). Geodesic distance between points is typically computed using geometric theorems, rules, or methods (e.g., the Pythagorean Theorem for computing the geodesic distance between two points in the plain); computations of geodesic distances between points may also be implemented in various global positioning systems, GPS, and diverse computer programs. In graph theory, a geodesic distance between two vertices is the number of edges in the shortest path connecting them. 
         [0010]    A list of numbers can also be transformed to multiple lists of numbers using spectral, wavelets, or polynomial transformations (Fourier 1822; Gheondea and Sabac 2003; Löpez-Gömez 2001; Mathworld 2006; Meyer 1992; Percival and Walden 2000). Spectral, wavelets, or polynomial transformations of lists of numbers typically change the numbers in each list and increase the number of lists. For example, a 2-degrees orthonormal polynomial transformation of list x that is denoted by x=(x 1 , . . . , x n ) for n numbers x i  where i=1:3 such that x=(1, 2, 3) produces two lists of numbers: x′ 1 =(−0.7071068, 0.0000000, 0.7071068) and x′ 2 =(0.4082483, −0.8164966, 0.4082483) that are also two columns of a 3*2 matrix. The orthonormal polynomial transformation in this example illustrates a transformation of one list to two lists and also a transformation of one list to a matrix. Orthonormality in this example indicates that the orthonormal polynomial transformation produced lists wherein each list has specific internal characteristics (e.g., normality) and wherein inter-relationships among lists have specific characteristics (e.g., orthogonality). 
         [0011]    A transformation of a list of numbers to another list of numbers replaces these numbers by other numbers. For example, if an original list that is denoted by x=(x 1 , . . . , x n ) for n numbers x i  where i=1:3 such that x=(1, 2, 3) is transformed to a new list that is denoted by x′=(x′ 1 , . . . , x′ n ) for n numbers x′ i  using x′ i =x i   2  then x′=(1, 4, 9), which is an example of a power transformation of list x using a power of 2. This is also an example of a transformation of one list to another list. Power transformations are one of the most popular methods of transformation of a list of numbers to another list of numbers (Emerson 1983; Emerson and Stoto 1983). 
         [0012]    Multivariate analyses and data mining focus on more than one list, seeking to classify, reduce, make sense of, or handle or manage a plurality of lists. They typically do so by reducing the number of lists. For example, factor analysis or principal components analysis change the numbers in each list and reduce the number of lists (Hotelling 1933; Pearson 1901; Rencher 2002; Spearman 1904), and cluster analysis enables analysts to identify groups within lists (Everitt 1993; Hartigan 1975). 
         [0013]    The foregoing considerations, and studies of transformations of lists, reveal that a list of numbers could change part, or no part, or all of the components of this list. A transformation of a list could also leave intact or increase or reduce the number of lists. For example, a power or logarithmic transformation of one list of numbers typically changes these numbers. and power or logarithmic transformations of more than one list of numbers tend to leave intact the number of lists. An identity transformation of a list of numbers (e.g., a power transformation with a power of 1) does not change the numbers of this list, and identity transformations of lists of numbers do not change the number of lists. 
         [0014]    Transformations of a list of numbers could offer diverse benefits in data analysis. For example, transformations of lists of numbers could linearize data or generate nonlinear data, increase symmetry or normality, increase homoscedasticity or reduce heteroscedasticity, orthogonalize or otherwise reposition variables, deepen insight into the data, simplify analysis, reveal new dimensions, and reveal meanings (Emerson 1983; Emerson and Stoto 1983; Rencher 2002). 
         [0015]    Because they have diverse beneficial uses, transformations can be evaluated, compared, and optimized in order to assess which transformation would produce optimal results. For example, it is possible to compare various power or logarithmic transformations of a list of numbers for the purpose of discovering which of these transformations would produce an optimal result in a specific data analysis (Box and Cox 1964; Emerson 1983; Emerson and Stoto 1983). 
         [0016]    While some transformation methods (e.g., power or logarithmic transformations) are simple or straight forward, or mechanized to diverse degrees, other transformation methods are complicated. For example, spectral or wavelets analyses typically require special transforms and complex calculations that limit the scope and applicability of these methods. Another drawback of the available methods of transformation is that they may not allow the user sufficient control, maneuverability, and tractability in the performance of the transformation or analysis. Additionally, some transformations may lack a theoretical foundation or a methodological justification that would explain why specific transformations ought to be used. 
         [0017]    The multiple lists of numbers that result from spectral, wavelets, or polynomial transformations of one or more lists of numbers could, in turn, be used to multivariate analyses or data mining (Gheondea and Sabac 2003; Löpez-Gömez 2001; Meyer 1992; Percival and Walden 2000). 
         [0018]    The present inventor published a study, cited herein as Epelbaum (2005), showing specifications and transformations of a list of numbers or a list of points to respective matrices of distances or proximities, wherein said distances or proximities are distances or proximities among components within this list. Epelbaum (2005) also shows multivariate analyses with such matrices of distances or proximities. Additionally, Epelbaum (2005) shows that variations of the specifications and transformations of lists to matrices and variations of the specifications and transformations of the matrices to other matrices can be manipulated so as to increase the usefulness of the results of multivariate analyses. 
         [0019]    Epelbaum (2005) presents factor analyses of matrices that have been produced with transformations of lists to matrices. In one study, Epelbaum (2005) incorporates factors that result from such factor analyses in the investigation of effects of spatial configuration on income per capita in U.S. states in the year 2000. In another study, Epelbaum (2005) incorporates factors that result from such factor analyses in the investigation of changes in solar flux during the first 100 days of 1981. In these investigations, Epelbaum (2005) also employs matrices that result from transformations of lists as working correlations structures in Generalized Estimating Equations, GEE, analyses. GEEs and other regression methods, like mixed-effects models, employ working correlations structures in data analyses of dependent or clustered data (Hardin and Hilbe 2003; Heagerty and Zeger 2000; Liang and Zeger 1986; Pinheiro and Bates 2000). 
         [0020]    The present invention formalizes, generalizes, refines, extends, broadens, and adds new uses to what the inventor presented in Epelbaum (2005). 
       BRIEF SUMMARY OF THE INVENTION 
       [0021]    The disclosed invention provides a method for transforming lists to corresponding matrices and optimizing results of data analyses. 
     
     
       BRIEF DESCRIPTION OF THE DRAWINGS 
         [0022]    The present invention can be more fully understood by reading the subsequent detailed description in conjunction with the examples and references made to the accompanying drawings. 
           [0023]      FIG. 1  provides a flowchart for an embodiment of the disclosed invention. 
           [0024]      FIG. 2  illustrates sections of some of the lists that were transformed to matrices in the embodiment of the invention. 
           [0025]      FIG. 3  presents the optimal specifications of the transformations in the embodiment of the invention. 
           [0026]      FIG. 4  presents a depiction of the structure of distribution of income per capita in U.S. states 2000. The data analysis in the embodiment of the invention shows that this structure of income per capita affected infant mortality in U.S. states 2000. This result from the data analyses in the embodiment of the invention is more valid and reliable than results of pervious research; this result is also useful for public health policy and intervention. 
           [0027]      FIG. 5  presents a depiction of effects of income per capita and its distributional structure on infant mortality in U.S. states 2000. This result from the data analyses in the embodiment of the invention is more valid and reliable than results of pervious research; this result is also useful for public health policy and intervention. 
       
    
    
     DETAILED DESCRIPTION 
       [0028]    The disclosed invention is embodied in data analyses of effects of race, gender, income, educational attainment, state of residence, and population size on infant mortality in U.S. states 2000. 
         [0029]    Previous research shows that data analyses of effects of race, gender, income, educational attainment, state of residence, and population size on mortality in U.S. states 2000 are useful for public health policy and intervention. However, the validity and reliability of previous research is limited because few studies allow and accommodate correlated response data, unobserved heterogeneity, overdispersion, and correlated influential factors, and because previous research does not provide mathematical specifications and graphic depictions of multifaceted, multidimensional, nonlinear, and distinct effects on mortality. The embodiment of the invention could improve the validity and reliability of data analyses of mortality, and thus increase their usefulness for public health policy and intervention. 
         [0030]    The present embodiment of the invention employs data that are originally available in lists of numbers, points, and categories. This embodiment follows the steps that are detailed in the flowchart in  FIG. 1 . Selections from some of the original lists for this embodiment are illustrated in  FIG. 2 . The original lists include a list of counts of dead infants and a list of counts of live births by 3 categories of race (i.e., African American, White, or Other), 2 categories of gender (i.e., female or male), and 51 categories of U.S. state of residence (i.e., 51 U.S. states, including the District of Columbia) in U.S. states 2000. Counts of deaths and counts of live births were compiled in two lists of 3*2*51=306 cases of counts of deaths and live births of infants in U.S. states 2000. Thus, the length n of the list of counts of deaths is n=306, and the length n of the list of counts of live births is also n=306, where n denotes the length of each list, which is also the number of units for the data analysis in the embodiment of the invention. Three corresponding lists of categories of race, gender, and U.S. state of residences were constructed, and were also set to lengths n=306. Each position in these lists was designated as a unit of analysis with an integer u, where u=1:n. Each list also constitutes a variable. The count of deaths in u was designated as variable D u , where u refers to a unit of data analysis, race in u was designated with variable R u , gender in u was designated with variable G u , state of residence in u was designated with variable S u , and the count of live births in u was designated with variable L u . These variables correspond, respectively, to lists that were designated as d, r, g, s, and l. 
         [0031]    Additional original lists for data analysis include state-level aggregate data for the year 2000 from the U.S. Census for the year 2000 (U.S. Census 2006), including: income per capita of all persons, designated with IPCAP, in 51 states; race-specific income per capita, designated with RSIPC, in 51 states by three racial categories; proportion of a state&#39;s households in interval q of household income, designated with HA q , where q=1:16 in 51 states; proportion of a state&#39;s race-specific households in interval q of household income, designated with HR q , where q=1:16 in 51 states by three racial categories; proportion of a state&#39;s population aged 25 or older in interval q of educational attainment, designated with EA q , where q=1:16 in 51 states; and proportion of a state&#39;s race-and-gender-specific population aged 25 or older in interval q of educational attainment, designated with RGSEA q , where q=1:7 in 51 states by three racial categories by two gender categories. In preparation for data analyses, the state level lists are compiled in lists of variable lengths. These lists include one list of IPCAP of length n=51, one list of RSIPC of length n=153, 16 lists of HA qu  of length n=51, 16 lists of HR q  of length n=153, 16 lists of EA q  of length n=51, and 7 lists of RGSEA q  of length n=306. 
         [0032]    Previous research shows that mortality models tend to be characterized by heterogeneity and extra-Poisson variation (Brillinger 1986; Manton, Woodbury, and Stallard 1981, 1983; Manton and Stallard 1984a, 1984b). Therefore, the negative binomial was used in the model for infant mortality, allowing and accommodating heterogeneity and extra-Poisson variation through employment of the negative binomial&#39;s alpha coefficient and gamma function. These procedures follow a previous study by Manton, Woodbury, and Stallard (1981) and subsequent considerations (Brillinger 1986; Lawless 1987; Manton, Woodbury, and Stallard 1983; Manton and Stallard 1984a, 1984b; Mullahy 1997). 
         [0033]    To account for dependence among mortality counts, a negative binomial Generalized Estimating Equation, GEE, model was formulated. GEEs and the employment of the working correlations structures in data analyses of dependent or clustered data are described in Hardin and Hilbe (2003). 
         [0034]    A negative binomial GEE model for effects on D, was formulated following previous work that is presented in Epelbaum (2005), to account for correlated influential factors, and to provide mathematical specifications and graphic depictions of multifaceted, multidimensional, nonlinear, and distinct effects on infant mortality. This model is formulated with 
         [0000]    
       
         
           
             
               
                 D 
                 u 
               
               = 
               
                 
                   exp 
                    
                   
                     ( 
                     α 
                     ) 
                   
                 
                 * 
                 
                   
                     ∏ 
                     
                       v 
                       = 
                       1 
                     
                     V 
                   
                    
                   
                       
                   
                    
                   
                     
                       exp 
                        
                       
                         ( 
                         
                           
                             β 
                             v 
                           
                            
                           
                             F 
                             vu 
                           
                         
                         ) 
                       
                     
                     * 
                     
                       exp 
                        
                       
                         ( 
                         
                           ɛ 
                           u 
                         
                         ) 
                       
                     
                   
                 
               
             
             , 
           
         
       
     
         [0000]    for effects of orthogonal factors F v  on infant death counts, D u , where α and β v  are GEE coefficients from the negative binomial GEE model, and where ε u  is the error term for unit u of this model. The model postulates V orthogonal factors F v  that were constructed with the embodiment of the invention. The GEE analysis also employs a working correlations matrix to account for dependency among infant death counts. The working correlations matrix was also constructed with the embodiment of the invention. The embodiment of the invention transformed the univariate lists that are presented above to the multivariate matrices that are presented next. 
         [0035]    Univariate lists were transformed to multivariate matrices for use in the right-hand side of the negative binomial GEE model for effects on infant mortality. Most of the same univariate lists were also transformed to other matrices for use in the working correlations matrix of the same model. One univariate list, list l that is presented next, was transformed to a matrix that was used in the right-hand side of this model and in the working correlation matrix of this model. 
         [0036]    List l of length n=306 for counts of live births (i.e., l=11,109; 261; . . . ; 3,209), parts of which are illustrated in column  8  of  FIG. 2 , was transformed to list l′ using l′ i =l i   φ  for i=1:n, and list l′ was transformed to a square matrix D of n 2  distances d jk  using d jk =|l′ j −l′ k | for j=1:n and k=1:n, and matrix D was transformed to its corresponding matrix P of proximities p jk  using p jk =[1−d jk /max(d . . . )] θ  for j=1:n and k=1:n where max(d . . . ) refers to the maximal element d jk  in matrix D. Matrix P of size n 2 =306 2  was designated as matrix L for proximities among counts of live births. Matrix L has 306 columns of length  306 , which were treated as 306 lists of length  306 . Thus, a univariate list l was transformed to a multivariate matrix L that has two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses. 
         [0037]    List g of length n=306 for gender (i.e., g=Male, Male, . . . Female), parts of which are illustrated in column  4  of  FIG. 2 , was transformed to two identical lists f and m using identity transformations. List f for females, was subsequently transformed to list f′ using Female=0.5 and Male=0 resulting in numbers f′ i  where i=1:n. Similarly, list m, for males, was transformed to list m′ using Male=0.5 and Female=0 resulting in numbers m′ i  where i=1:n. Numerical lists f′ and m′ were subsequently retransformed to two respective columns of a matrix FM, that was subsequently retransformed to matrix G using G=(FM) θ . Thus, a univariate list g of categories was transformed to a multivariate matrix G of numbers having one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0038]    List r of length n=306 for race (i.e., r=African American, African American, . . . White), parts of which are illustrated in column  3  of  FIG. 2 , was transformed to three identical lists ra, ro, and rw using identity transformations. List ra, for African Americans, was subsequently transformed to list ra′ using African American=0.5, Other=0, and White=0 resulting in numbers ra′ i  where i=1:n. Similarly, list ro, for Other, was subsequently transformed to list ro′ using African American=0, Other=0.5, and White=0 resulting in numbers ro′ i  where i=1:n. Additionally, list rw, for White, was subsequently transformed to list rw′ using African American=0, Other=0, and White=0.5 resulting in numbers rw′ i  where i=1:n. Numerical lists ra′, ro′, and rw′ were subsequently retransformed to three respective columns of a matrix AOW, that was subsequently retransformed to matrix R using R=(AOW) θ . Thus, a univariate list r of categories was transformed to a multivariate matrix R of numbers having one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0039]    List s of length n=51 with 51 categories of names of U.S. states of residence (i.e., s=Alabama, Alaska, . . . , Wyoming), is partly illustrated in column  2  of  FIG. 2 . This list was transformed to a corresponding list sp of 51 points that are indicated by centroids of U.S. states (i.e., sp=33.00N&amp;-86.77W, 61.29N&amp;-148.72W, . . . 42.68N&amp;-107.01W), with data from the U.S. Census (2001), and is partly illustrated in columns  5  and  6  of  FIG. 2 . List sp was transformed to a square matrix D of n 2  geodesic distances d jk  for j=1:n and k=1:n using a computer program written by Simons (2004) for computation of geodesic distances among points on earth. Matrix D was transformed to a corresponding matrix D′ of distances d′ jk  using d′ jk =(d jk ) φ  for j=1:n and k=1:n, and matrix D′ was transformed to a corresponding matrix P of proximities p jk  using p jk =[1−d′ jk /max(d′ . . . )] θ  for j=1:n and k=1:n where max(d′ . . . ) refers to the maximal element d′ jk  in matrix D′. Matrix P of size n 2 =51 2  was designated as matrix S for spatial proximities among U.S. states. Six identical matrices S were row-bound and were thus transformed to a new matrix S′ of 51 columns of length  306 , which were treated as 51 lists of length  306 . Thus, a univariate list s of categories was transformed to a multivariate matrix S′ of numbers that has two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses. 
         [0040]    List ipc of length n=51, for income per capita of all persons by U.S. state 2000 (i.e., ipc=18,188.85; 22,659.75; . . . ; 19,133.95), was transformed to list ipc′ using ipc′ i =ipc i   φ  for i=1:n, and list ipc′ was transformed to a square matrix D of n 2  distances d jk  using d jk =|ipc′ j −ipc′ k | for j=1:n and k=1:n, and matrix D was transformed to its corresponding matrix P of proximities p jk  using p jk =[1−d jk /max(d . . . )] θ  for j=1:n and k=1:n where max(d . . . ) refers to the maximal element d jk  in matrix D. Matrix P of size n 2 =51 2  was designated as matrix IPC for proximities among income per capita in 51 U.S. states. Six identical matrices IPC were row-bound and were thus transformed to a new matrix IPC′ of 51 columns of length  306 , which were treated as 51 lists of length  306 . Thus, a univariate list ipc was transformed to a multivariate matrix IPC′ that has two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses. 
         [0041]    List rsipc of length n=153, for race specific income per capita of all persons by U.S. state 2000 (i.e., rsipc=11,664.61; 18,512.93; . . . ; 19,700.16), was transformed to list rsipc′ using rsipc′ i =rsipc i   φ  for i=1:n, and list rsipc′ was transformed to a square matrix D of n 2  distances d jk  using d jk =|rsipc′ j −rsipc′ k | for j=1:n and k=1:n, and matrix D was transformed to its corresponding matrix P of proximities p jk  using p jk =[1−d jk /max(d . . . )] θ  for j=1:n and k=1:n where max(d . . . ) refers to the maximal element d jk  in matrix D. Matrix P of size n 2 =153 2  was designated as matrix RSIPC for proximities among income per capita in 51 U.S. states. Two identical matrices RSIPC were row-bound and were thus transformed to a new matrix RSIPC′ of 153 columns of length  306 , which were treated as 153 lists of length  306 . Thus, a univariate list rsipc was transformed to a multivariate matrix RSIPC′ that has two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses. 
         [0042]    Sixteen lists ha q  of length n=51 for proportion of a state&#39;s households in interval q of household income, where q=1:16 (i.e., ha 1 =0.144, 0.056, . . . 0.091; . . . ; ha 16 =0.014, 0.019, . . . 0.013), were transformed to a matrix HA of size 51*16 using ha q  as respective 16 columns of said matrix HA, where q=1:16. Matrix HA was transformed to matrix HA′ using HA′=(HA) θ , and six identical matrices HA′ were row-bound and were thus transformed to a new matrix HA″ of 16 columns of length  306 , which were treated as 16 lists of length  306 . Thus, sixteen univariate lists ha q  were transformed to a multivariate matrix HA″ that has one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0043]    Sixteen lists hr q  of length n=153 for proportion of a state&#39;s race-specific households in interval q of race-specific household income, where q=1:16 (i.e., hr 1 =0.25, 0.060, . . . 0.087; . . . ; hr 16 =0.006, 0.006, . . . 0.013), were transformed to a matrix HR of size 153*16 using hr q  as respective 16 columns of said matrix HR, where q=1:16. Matrix HR was transformed to matrix HR′ using HR′=(HR) θ , and two identical matrices HR′ were row-bound and were thus transformed to a new matrix HR″ of 16 columns of length  306 , which were treated as 16 lists of length  306 . Thus, sixteen univariate lists hr q  were transformed to a multivariate matrix HR″ that has one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0044]    Sixteen lists ea q  of length n=51 for proportions of a state&#39;s population aged 25 or older in interval q of educational attainment, where q=1:16 (i.e., ea 1 =0.012, 0.008, . . . 0.004; . . . ; ea 16 =0.007, 0.009, . . . 0.008), were transformed to a matrix EA of size 51*16 using ea q  as respective 16 columns of said matrix EA, where q=1:16. Matrix EA was transformed to matrix EA′ using EA′=(EA) θ , and six identical matrices EA′ were row-bound and were thus transformed to a new matrix EA″ of 16 columns of length  306 , which were treated as 16 lists of length  306 . Thus, sixteen univariate lists ea q  were transformed to a multivariate matrix EA″ that has one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0045]    Seven lists rgsea q  of length n=306 for proportions of a state&#39;s race-and-gender-specific population aged 25 or older in interval q of educational attainment, where q=1:7 (i.e., rgsea 1 =0.114, 0.026, . . . 0.025; . . . ; rgsea 7 =0.031, 0.051, . . . 0.059), were transformed to a matrix RGSEA of size 306*7 using rgsea q  as respective 7 columns of said matrix RGSEA, where q=1:7. Matrix RGSEA was transformed to matrix RGSEA′ using RGSEA′=(RGSEA) θ , resulting in a new matrix RGSEA′ of 7 columns of length  306 , which were treated as 7 lists of length  306 . Thus, seven univariate lists rgsea q  were transformed to a multivariate matrix RGSEA′ that has one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0046]    The matrices that were presented above were employed in the right-hand side of the negative binomial GEE model for effects on infant mortality. The matrices that are presented next were employed in the working correlations matrix of this model. 
         [0047]    List g of length n=306 for gender (i.e., g=Male, Male, . . . Female), parts of which are illustrated in column  4  of  FIG. 2 , was transformed to list g′ using Female=2 and Male=1 resulting in numbers g′ i  where i=1:n. List g′ was transformed to a square matrix D of n 2  distances d jk  using d jk =|g′ j −g′ k | for j=1:n and k=1:n. Matrix D was transformed to its corresponding matrix P of proximities p jk  using p jk =0.5 when d jk =1, and p jk =1 when d jk =0, for j=1:n and k=1:n. Matrix P was transformed to its corresponding matrix P′ of proximities p′ jk  using p″ jk =(p′ jk ) θ  for j=1:n and k=1:n. Matrix P′ was designated as matrix G ▪  for proximities among gender categories. Matrix G ▪  is a square matrix of 306 rows, 306 columns, and a diagonal of ones. Thus, a univariate list g was transformed to a multivariate matrix G ▪  with a manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0048]    List r of length n=306 for race (i.e., r=African American, African American, . . . White), parts of which are illustrated in column  3  of  FIG. 2 , was transformed to list r′ using African American=3, Other=2, and White=1 resulting in numbers r′ i  where i=1:n. List r′ was transformed to a square matrix D of n 2  distances d jk  using d jk =|r′ j −r′ k | for j=1:n and k=1:n. Matrix D was transformed to its corresponding matrix P of proximities p jk  using p jk =0.5 when d jk &gt;0, and using p jk =1 when d jk =0, for j=1:n and k=1:n. Matrix P was transformed to its corresponding matrix P′ of proximities p′ jk  using p″ jk =(p′ jk ) θ  for j=1:n and k=1:n. Matrix P′ was designated as matrix R ▪  for proximities among racial categories. Matrix R ▪  is a square matrix of 306 rows, 306 columns, and a diagonal of ones. Thus, a univariate list r was transformed to a multivariate matrix R ▪  with a manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0049]    Column-binding of six identical matrices S′, each with 51 lists of length  306  and two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses, resulted in their transformation to a new matrix S ▪ . Matrix S ▪  has 306 columns of length  306  and one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. Thus, a univariate list s of categories was transformed to a multivariate matrix S ▪  of numbers that has two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses. 
         [0050]    Column-binding of six identical matrices IPC′, each with 51 lists of length  306  and two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses, resulted in their transformation to a new matrix IPC ▪ . Matrix IPC′ has 306 columns of length  306  and two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses. Thus, a univariate list ipc of numbers was transformed to a multivariate matrix IPC ▪  that has two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses. 
         [0051]    Column-binding of two identical matrices RSIPC′, each with 153 lists of length  306  and two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses, resulted in their transformation to a new matrix RSIPC ▪  of 306 columns of length  306  and two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses. Thus, a univariate list rsipc of numbers was transformed to a multivariate matrix RSIPC ▪  that has two manipulable mathematical coefficients φ and θ whose specified values can be varied in data analyses. 
         [0052]    Each of the sixteen lists ha q  of length n=51 for proportion of a state&#39;s households in interval q of household income, where q=1:16 (i.e., ha 1 =0.144, 0.056, . . . 0.091; . . . ; ha 16 =0.014, 0.019, . . . 0.013), was transformed to a square matrix D q  of n 2  distances d qjk  using d qjk =|ha qj −ha qk | for j=1:n and k=1:n, and matrix D q  was transformed to its corresponding matrix P q  of proximities p qjk  using p qjk =[1−d qjk /max(d . . . )] for j=1:n and k=1:n where max(d . . . ) refers to the maximal element d qjk  in matrix D q . Row-binding matrix P q  of 51 columns of lengths  51  six times produces a new matrix P q ′ of 51 columns of lengths  306 . Column-binding matrix P q ′ of 51 columns of lengths 306 six times produces a new matrix P q ″ of 306 columns of lengths  306 . Repeating these transformations 16 times for q=1:16, produces 16 matrices P q ″ of 306 columns of lengths  306  where q=1:16. Adding the 16 matrices P q ″ produces a new matrix that was designated here as matrix HA that has 306 columns of length  306 , which were treated as 306 lists of length  306 . Matrix HA was transformed to matrix HA ▪  using HA ▪ =(HA) θ . Thus, 16 univariate lists ha q  were transformed to a multivariate matrix HA ▪  that has one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0053]    Each of the sixteen lists hr q  of length n=153 for proportion of a state&#39;s race-specific households in interval q of race-specific household income, where q=1:16 (i.e., hr 1 =0.25, 0.060, . . . 0.087; . . . ; hr 16 =0.006, 0.006, . . . 0.013), was transformed to a square matrix D q  of n 2  distances d qjk  using d qjk =|hr qj −hr qk | for j=1:n and k=1:n, and matrix D q  was transformed to its corresponding matrix P q  of proximities p qjk  using p qjk =[1−d qjk /max(d . . . )] for j=1:n and k=1:n where max(d . . . ) refers to the maximal element d qjk  in matrix D q . Row-binding matrix P q  of 153 columns of lengths  153  two times produces a new matrix P q ′ of 153 columns of lengths  306 . Column-binding matrix P q ′ of 153 columns of lengths  306  two times produces a new matrix P q ″ of 306 columns of lengths  306 . Repeating these transformations 16 times for q=1:16, produces 16 matrices P q ″ of 306 columns of lengths  306  where q=1:16. Adding the 16 matrices P q ″ produces a new matrix that was designated here as matrix HR that has 306 columns of length  306 , which were treated as 306 lists of length  306 . Matrix HR was transformed to matrix HR ▪  using HR ▪ =(HR) θ . Thus, 16 univariate lists hr q  were transformed to a multivariate matrix HR ▪  that has one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0054]    Each of the sixteen lists ea q  of length n=51 for proportions of a state&#39;s population aged 25 or older in interval q of educational attainment, where q=1:16 (i.e., ea 1 =0.012, 0.008, . . . 0.004; . . . ; ea 16 =0.007, 0.009, . . . 0.008), was transformed to a square matrix D q  of n 2  distances d qjk  using d qjk =|ea qj −ea qk | for j=1:n and k=1:n, and matrix D q  was transformed to its corresponding matrix P q  of proximities p qjk  using p qjk =[1−d qjk /max(d . . . )] for j=1:n and k=1:n where max(d . . . ) refers to the maximal element d qjk  in matrix D q . Row-binding matrix P q  of 51 columns of lengths 51 six times produces a new matrix P q ′ of 51 columns of lengths  306 . Column-binding matrix P q ′ of 51 columns of lengths  306  six times produces a new matrix P q ″ of 306 columns of lengths  306 . Repeating these transformations 16 times for q=1:16, produces 16 matrices P q ″ of 306 columns of lengths  306  where q=1:16. Adding the 16 matrices P q ″ produces a new matrix that was designated as matrix EA that has 306 columns of length  306 , which were treated as 306 lists of length  306 . Matrix EA was transformed to matrix EA ▪  using EA ▪ =(EA) θ . Thus, 16 univariate lists ea q  were transformed to a multivariate matrix EA ▪  that has one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0055]    Each of the seven lists rgsea q  of length n=306 for proportions of a state&#39;s race-and-gender-specific population aged 25 or older in interval q of educational attainment, where q=1:7 (i.e., rgsea 1 =0.114, 0.026, . . . 0.025; . . . ; rgsea 7 =0.031, 0.051, . . . 0.059), was transformed to a square matrix D q  of n 2  distances d qjk  using d qjk =|rgsea qj −rgsea qk | for j=1:n and k=1:n, and matrix D q  was transformed to its corresponding matrix P q  of proximities p qjk  using p qjk =[1−d qjk /max(d . . . )] for j=1:n and k=1:n where max(d . . . ) refers to the maximal element d qjk  in matrix D q . Repeating these transformations 7 times for q=1:7, produces 7 matrices P q ″ of 306 columns of lengths  306  where q=1:7. Adding the 7 matrices P q ″ produces a new matrix that was designated as matrix RGSEA that has 306 columns of length  306 , which were treated as 306 lists of length  306 . Matrix RGSEA was transformed to matrix RGSEA ▪  using RGSEA ▪ =(RGSEA) θ . Thus, 7 univariate lists rgsea q  were transformed to a multivariate matrix RGSEA ▪  that has one manipulable mathematical coefficient θ whose specified values can be varied in data analyses. 
         [0056]    A 306*306 matrix with a diagonal of ones and 0.5 in all other cells was constructed and designated as matrix HH, in reference to homogeneity and heterogeneity. Matrix HH was then transformed to matrix HH ▪  using HH ▪ =(HH) θ . Matrix HH ▪  is constructed under the assumption that infants&#39; mortality in the 306 units u is partly homogeneous in so far as infants&#39; mortality responses in units u in the U.S. are similar to each other and differ from other phenomena. For example, infants&#39; mortality responses in units u in the U.S. differ from corresponding mortality responses of other age groups in the U.S. or elsewhere, from mortality responses elsewhere, or from mortality responses of other species. 
         [0057]    The 306 lists of matrix L, 2 lists of matrix G, 3 lists of matrix R, 51 lists of matrix S′, 51 lists of matrix IPC′, 153 lists of matrix RSIPC′, 16 lists of matrix HA″, 16 lists of matrix HR″, 16 lists of matrix EA″, and 7 lists of matrix RGSEA′ included specific manipulable mathematical coefficients φ or θ. The resulting specified 621 lists of length  306  were column-bounded in a new broad matrix that was designated as matrix B. Thus, the manipulable mathematical coefficients φ and θ of matrices L, G, R, S′, IPC′, RSIPC′, HA″, HR″, EA″, and RGSEA′ were incorporated into matrix B. Therefore, specified values of manipulable mathematical coefficients φ and θ could be manipulated in GEE analyses with matrix B. 
         [0058]    A 306*306 matrix C was constructed using 
         [0000]    
       
         
           
             
               
                 
                   C 
                   = 
                     
                    
                   
                     [ 
                     
                       ( 
                       
                         
                           
                             w 
                             l 
                             * 
                           
                            
                           L 
                         
                         + 
                         
                           
                             w 
                             g 
                             * 
                           
                            
                           
                             G 
                             • 
                           
                         
                         + 
                         
                           
                             w 
                             r 
                             * 
                           
                            
                           
                             R 
                             • 
                           
                         
                         + 
                         
                           
                             w 
                             s 
                             * 
                           
                            
                           
                             S 
                             • 
                           
                         
                         + 
                         
                           
                             w 
                             ipc 
                             * 
                           
                            
                           
                             IPC 
                             • 
                           
                         
                         + 
                       
                     
                   
                 
               
             
             
               
                 
                     
                    
                   
                     
                       
                         w 
                         rsipc 
                         * 
                       
                        
                       
                         RSIPC 
                         • 
                       
                     
                     + 
                     
                       
                         w 
                         ha 
                         * 
                       
                        
                       
                         HA 
                         • 
                       
                     
                     + 
                     
                       
                         w 
                         hr 
                         * 
                       
                        
                       
                         HR 
                         • 
                       
                     
                     + 
                     
                       
                         w 
                         ea 
                         * 
                       
                        
                       
                         EA 
                         • 
                       
                     
                     + 
                   
                 
               
             
             
               
                 
                   
                     
                       
                         
                             
                            
                           
                             
                               
                                 w 
                                 rsgsea 
                                 * 
                               
                                
                               
                                 RGSEA 
                                 • 
                               
                             
                             + 
                             
                               
                                 w 
                                 hh 
                                 * 
                               
                                
                               
                                 HH 
                                 • 
                               
                             
                           
                           ) 
                         
                         / 
                         W 
                       
                       ] 
                     
                     Ω 
                   
                   , 
                 
               
             
           
         
       
       
         
           where 
         
       
       
         
           
             
               
                 
                   W 
                   = 
                     
                    
                   
                     ( 
                     
                       
                         w 
                         l 
                       
                       + 
                       
                         w 
                         g 
                       
                       + 
                       
                         w 
                         r 
                       
                       + 
                       
                         w 
                         s 
                       
                       + 
                       
                         w 
                         ipc 
                       
                       + 
                       
                         w 
                         rsipc 
                       
                       + 
                     
                   
                 
               
             
             
               
                 
                   
                       
                      
                     
                       
                         w 
                         ha 
                       
                       + 
                       
                         w 
                         hr 
                       
                       + 
                       
                         w 
                         ea 
                       
                       + 
                       
                         w 
                         rsgsea 
                       
                       + 
                       
                         w 
                         hh 
                       
                     
                     ) 
                   
                   . 
                 
               
             
           
         
       
     
         [0000]    The manipulable mathematical coefficients φ, θ, Ω, w, and W that are associated with matrices L, G ▪ , R ▪ , S ▪ , IPC ▪ , RSIPC ▪ , HA ▪ , HR ▪ , EA ▪ , RGSEA ▪ , and HH ▪  were incorporated into matrix C. Therefore, specified values of manipulable mathematical coefficients φ, θ, Ω, w, and W could be manipulated in GEE analyses with matrix C. 
         [0059]    Matrix B was used in the construction of orthogonal factors F vu  in the right-hand side of the negative binomial GEE model for infant mortality. Matrix C was used as the working correlations matrix for dependent data in this GEE model. 
         [0060]    Two optimization criteria were employed in the GEE analyses. One optimization criterion was Lin&#39;s ρ c  coefficient for the concordance correlation between observed and predicted D u  (Lin 1989, 2000; Hardin and Hilbe 2003). The second optimization criterion was the minimization Pan&#39;s quasi-likelihood information criterion, QIC, for assessing appropriateness of the working correlation structure in GEE analyses (Pan 2001; Hardin and Hilbe 2003). Therefore, optimality was achieved with a calibration of the maximization of Lin&#39;s ρ c  coefficient and the minimization of Pan&#39;s QIC 
         [0061]    Diverse types of factor analyses, with diverse rotations, and diverse factor scoring techniques, as well as principal components analysis, were conducted on matrix B, constructing 90 complex new variables of length  360  on the basis of each analysis. General discussions on diverse procedures of factor analysis and principal components analysis are discussed in StataCorp (2001) and Rencher (2002). Each factor analysis or principal components analysis was followed by negative binomial regression analyses of effects of selected complex variables—selected from among the 90 complex new variables—on infant mortality, D u . Lin&#39;s ρ c  coefficient for the concordance correlation between observed and predicted D u  was computed following each negative binomial regression analysis (Lin 1989, 2000; Hardin and Hilbe 2003). Said analyses were repeated with variable specifications of manipulable mathematical coefficients φ or θ in the component matrices of matrix B. These analyses showed that factor analyses with principal factors, varimax rotation, and Thomson scoring maximize Lin&#39;s ρ c  coefficient. Subsequently, these factor analytic methods were used in all further factor analyses. 
         [0062]    Additional factor analyses with principal factors, varimax rotation, and Thomson scoring were conducted on matrix B, using variable specifications of manipulable mathematical coefficients φ or θ in the component matrices of matrix B in each factor analysis. Each of the following manipulable mathematical coefficients in the component matrices of matrix B was varied more than 6 times: φ and θ in matrix L, θ in matrix G, θ in matrix R, φ and θ in matrix S′, φ and θ in matrix IPC′, φ and θ in matrix RSIPC′, θ in matrix HA″, θ in matrix HR″, θ in matrix EA″, and θ in matrix RGSEA′. Variations in specification of one manipulable mathematical coefficient φ or θ proceeded as follows:
       one variation of each manipulable mathematical coefficient φ or θ established a basis; a second variation determined if an increase leads to improved results;   if the increase of the second variation of the coefficient led to improved results then an increase in the coefficient in a third variation determined if a further increase improved results, and so on, leading to further tests of increases in the manipulable mathematical coefficient until no further increase leads to improved results;   if the increase of the second variation of the coefficient did not change results then a smaller increase in the coefficient in a third variation determined if smaller increase improved results, and so on, otherwise decreases are tried, and so on;   if the increase of the third variation of the coefficient led to weaker results then a decrease in the coefficient in a fourth variation determined if a further decrease improved results, and so on, leading to further tests of decreases in the manipulable mathematical coefficient until no further decrease led to improved results.       
 
       At least two such passes were performed through all the manipulable mathematical coefficients φ or θ in the component matrices of matrix B, involving hundreds of factor analyses and negative binomial GEE analyses. 
       [0067]    Repeated negative binomial GEE analyses of effects of ten factors F vu  on D u  were conducted to optimize the calibration of the maximization of Lin&#39;s ρ c  coefficient for the concordance correlation and the minimization of Pan&#39;s QIC for the correlation structure of each respective negative binomial GEE analysis. The optimal calibration was achieved with QIC=67,816.55 and ρ c =0.99050 using the specified values of manipulable mathematical coefficients φ, θ, Ω, w, and W for the matrices that were inputted into matrices B and C and are shown in columns  3 ,  4 , and  5  of  FIG. 3 . 
         [0068]    The factor analysis that preceded the optimal parsimonious negative binomial GEE analysis was conducted on the 621 specified lists of matrix B, using the specified values of manipulable mathematical coefficients φ or θ of the respective matrices that were inputted matrix B and are shown in column  3  of  FIG. 3 . Therefore, a significant factor loading of one of the 621 variables was assessed using λ=γ/Σγ, where γ designates the absolute factor loading of one of the 621 factors loadings and Σγ designates the sum of the 621 γ coefficients. A significant factor loading on a factor corresponding to F vu  is indicated here by λ&gt; 1/621, noting that 1/621=0.00161=1.61E−3. Significant effect of an original variable (e.g., income per capita) on infant mortality is indicated by a {(mean λ)&gt;( 1/621)} in the factor loadings that corresponded to a factor F vu  from the optimal parsimonious negative binomial GEE model. 
         [0069]    To assess effects of income per capita of all persons, IPCAP, on infant mortality, an analysis was made of the factor loadings of the 51 IPCAP variables among the 621 factor loadings of each of the ten factor loadings that correspond to the Factors F v  in the optimal parsimonious negative binomial GEE model of effects on infant mortality. This inspection revealed that only factor loadings on Factor F 6  were significant, as indicated by mean λ=7.98E−3 for the 51 IPCAP variables corresponding to Factor F 6  where expected mean λ= 1/621=1.61E−3. Moreover, factor loadings of the 51 IPCAP variables corresponding to Factor F 6  were also dominant factor loadings with 10 IPCAP variables in the top 10 λ coefficients of factor loadings. The dominance of IPCAP variables in factor loadings on Factor F 6  was strengthened by having 52.9 percent of their corresponding 51 absolute factor loadings characterized by γ&gt;0.5. 
         [0070]    The structure of the distribution of income per capita of all persons in U.S. states, IPCAP, corresponding to factor loadings on Factor F 6  is shown in  FIG. 4 . The smooth line in  FIG. 4  depicts the structure of distribution of IPCAP that affected infant mortality in U.S. states 2000. 
         [0071]    Using the optimal parsimonious negative binomial GEE model of effects on infant mortality and the results of analyses of significance and dominance of factor loadings for IPCAP variables, it is possible to model the effects of IPCAP on infant mortality with 
         [0000]        M   u =exp(β′ 6   F   6u ), 
         [0000]    where M u  refers to the infant mortality that is attributed to effects of IPCAP in unit u. This model for M u , based upon the optimal parsimonious negative binomial GEE model, shows the effects of the size of IPCAP and effects of structure of the distribution of IPCAP on infant mortality in U.S. states 2000. This model is depicted graphically by the smooth line in  FIG. 5 , depicting effects of IPCAP on infant mortality in U.S. states 2000. 
         [0072]    The embodiment of the invention in data analyses of effects on infant mortality in U.S. states 2000 allows and accommodates correlated response data, unobserved heterogeneity, overdispersion, and correlated influential factors. The analyses provide mathematical specifications and graphic depictions of (A) multifaceted, multidimensional, nonlinear, and distinct factor analytic structures of population size, income, educational attainment, race, gender, and spatial configuration, and (B) incorporation of these structures in respective multiple orthogonal standardized constructs that affected infant mortality in U.S. states 2000. These effects reveal multifaceted, multidimensional, nonlinear, multidirectional, and distinct patterns of infant mortality. 
         [0073]    The results that are presented in  FIG. 4  and  FIG. 5  provide the most valid and reliable results currently available on effects of size and structures of income per capita of all persons on infant mortality in U.S. states 2000. These results are useful for public health policy and interventions in the United States. Similar procedures were employed to analyze effects of population size, race, gender, educational attainment, and spatial configuration on infant mortality. These analyses are now being prepared for publication in the areas of epidemiology, public health, medicine, economics, geography, sociology, and political science and are to be used by local and federal governmental agencies, policy makers, insurance companies, and others. 
         [0074]    Finally, while the invention has been described by way of examples and in terms of the above, it is to be understood that the invention is not limited to the disclosed embodiment. On the contrary, it is intended to cover various modifications that would be apparent to those skilled in the art. 
       REFERENCES 
       [0075]    Box, G. E. P. and D. R. Cox. 1964. “An analysis of transformations.”  Journal of the Royal Statistical Society—Series B  26:211-252. 
         [0076]    Brillinger, David R. 1986. “The natural variability of vital rates and associated statistics.”  Biometrics  42:693-734. 
         [0077]    Emerson, John D. 1983. “Mathematical aspects of transformation.” in  Understanding Robust and Exploratory Data Analysis,  edited by D. C. Hoaglin, F. Mosteller, and J. W. Tukey. New York: Wiley. 
         [0078]    Emerson, John D. and Michael A. Stoto. 1983. “Transforming data.” Pp. 97-128 in  Understanding Robust and Exploratory Data Analysis,  edited by D. C. Hoaglin, F. Mosteller, and J. W. Tukey. New York: John Wiley. 
         [0079]    Epelbaum, Michael. 2005. “A new theory and a new method for temporal, spatial, and other correlated phenomena.”  Complexity  11:36-44. 
         [0080]    Everitt, B. S. 1993.  Cluster Analysis.  London: Edward Arnold. 
         [0081]    Fourier, Jean Baptiste Joseph. 1822.  Théorie Analytique de la Chaleur.  Paris: F Didot. 
         [0082]    Gheondea, Aurelian and Mihai Sabac. 2003.  Spectral Analysis and Its Applications: Ion Colojoara Anniversary Volume.  Bucharest: Theta. 
         [0083]    Hardin, James W., and Joseph Hilbe. 2003.  Generalized Estimating Equations.  Boca Raton, Fla.: Chapman &amp; Hall/CRC. 
         [0084]    Hartigan, J. A. 1975.  Clustering Algorithms.  New York: John Wiley. 
         [0085]    Hotelling, H. 1933. “Analysis of a complex of statistical variables into principal components.”  Journal of Educational Psychology  24:417-441, 498-520. 
         [0086]    Heagerty, Patrick J., and Scott L. Zeger. 2000. “Marginalized multilevel models and likelihood inference.”  Statistical Science  15:1-26. 
         [0087]    Liang, Kung-Yee, and Scott L. Zeger. 1986. “Longitudinal data analysis using generalized linear models.”  Biometrika  73:13-22. 
         [0088]    Lin, L. I-K. 1989. “A concordance correlation coefficient to evaluate reproducibility.”  Biometrics  45:255-268. 
         [0089]    Lin, L. I-K. 2000. “A note on the concordance correlation coefficient.”  Biometrics  56:324-325. 
         [0090]    Löpez-Gömez, Julián. 2001.  Spectral theory and nonlinear functional analysis.  Boca Raton, Fla.: Chapman &amp; Hall/CRC. 
         [0091]    Manton, G. Kenneth, and Eric Stallard. 1984a. “Heterogeneity and its effect on mortality measurement.” Pp. 265-301 in  Methodologies for the collection and analysis of mortality data,  edited by Jacques Vallin, John H. Pollard, and Larry Heligman. Liege, Belgium: Ordina Editions. 
         [0092]    Manton, G. Kenneth, and Eric Stallard. 1984b.  Recent trends in mortality analysis.  Orlando: Academic Press. 
         [0093]    Manton, G. Kenneth, Eric Stallard, and James W. Vaupel. 1986. “Alternative models for the heterogeneity of mortality risk among the aged.”  Journal of the American Statistical Association  81:635-644. 
         [0094]    Manton, G. Kenneth, Max A. Woodbury, and Eric Stallard. 1981. “A variance components approach to categorical data models with heterogeneous cell populations: Analysis of spatial gradients in lung cancer mortality rates in North Carolina counties.”  Biometrics  37:259-269. 
         [0095]    Manton, G. Kenneth, Max A. Woodbury, and Eric Stallard. 1983. “A model for cancer mortality rates.”  Biometrics: 531-532. 
         [0096]    Mathworld. 2006. http://mathworld.wolfram.com/Polynomial.html. Accessed September 2006. 
         [0097]    Meyer, Yves. 1992.  Ondelettes et Opérateurs.  Paris: Hermann. 
         [0098]    Mullahy, John. 1997. “Heterogeneity, excess zeros, and structure of count data models.”  Journal of Applied Econometrics  12:337-350. 
         [0099]    Pan, Wei. 2001. “Akaike&#39;s information criterion in generalized estimating equations.”  Biometrics  57:120-125. 
         [0100]    Pearson, Karl. 1901. “One lines and planes of closest fit to systems of points in space.”  Philosophical Magazine  6:559-572. 
         [0101]    Percival, Donald B. and Andrew T. Walden. 2000.  Wavelet Methods for Time Series Analysis.  Cambridge: Cambridge University Press. 
         [0102]    Pinheiro, Jose C., and Douglas M. Bates. 2000.  Mixed - Effects Models in S and S - Plus.  New York: Springer. 
         [0103]    Rencher, Alvin C. 2002.  Methods of Multivariate Analysis.  New York: J. Wiley. 
         [0104]    Simons, Kenneth L. 2004. “globdist,” http://www.rpi.edu/˜simonk/stata. Accessed April 2004. 
         [0105]    Spearman, C. 1904. “General intelligence objectively determined and measured.”  American Journal of Psychology  15:201-293. 
         [0106]    StataCorp. 2001.  Stata  7  Reference Manual.  College Station, Tex.: Stata Corporation. 
         [0107]    U.S. Census. 2001. “Population and Population Centroid by State: 2000.” http://www.census.gov/Press-Release/www/2001/statecenters.txt. Accessed April 2004. 
         [0108]    U.S. CDC. 2006. “CDC Wonder: Compressed Mortality File.” http://wonder.cdc.gov/mortSQL.html. Accessed July 2006.