Asymptotic robust inferences in the analysis of mean and covariance structures

Structural equation models are widely used in economic, socialand behavioral studies to analyze linear interrelationships amongvariables, some of which may be unobservable or subject to measurementerror. Alternative estimation methods that exploit different distributionalassumptions are now available. The present paper deals with issues ofasymptotic statistical inferences, such as the evaluation of standarderrors of estimates and chi--square goodness--of--fit statistics,in the general context of mean and covariance structures. The emphasisis on drawing correct statistical inferences regardless of thedistribution of the data and the method of estimation employed. A(distribution--free) consistent estimate of $\Gamma$, the matrix ofasymptotic variances of the vector of sample second--order moments,will be used to compute robust standard errors and a robust chi--squaregoodness--of--fit squares. Simple modifications of the usual estimateof $\Gamma$ will also permit correct inferences in the case of multi--stage complex samples. We will also discuss the conditions under which,regardless of the distribution of the data, one can rely on the usual(non--robust) inferential statistics. Finally, a multivariate regressionmodel with errors--in--variables will be used to illustrate, by meansof simulated data, various theoretical aspects of the paper.

Distributional equivalence and subcompositional coherence in the analysis of contingency tables, ratio-scale measurements and compositional data

We consider two fundamental properties in the analysis of two-way tables of positive data: the principle of distributional equivalence, one of the cornerstones of correspondence analysis of contingency tables, and the principle of subcompositional coherence, which forms the basis of compositional data analysis. For an analysis to be subcompositionally coherent, it suffices to analyse the ratios of the data values. The usual approach to dimension reduction in compositional data analysis is to perform principal component analysis on the logarithms of ratios, but this method does not obey the principle of distributional equivalence. We show that by introducing weights for the rows and columns, the method achieves this desirable property. This weighted log-ratio analysis is theoretically equivalent to spectral mapping , a multivariate method developed almost 30 years ago for displaying ratio-scale data from biological activity spectra. The close relationship between spectral mapping and correspondence analysis is also explained, as well as their connection with association modelling. The weighted log-ratio methodology is applied here to frequency data in linguistics and to chemical compositional data in archaeology.

Multiple correspondence analysis of a subset of response categories

In the analysis of multivariate categorical data, typically the analysis of questionnaire data, it is often advantageous, for substantive and technical reasons, to analyse a subset of response categories. In multiple correspondence analysis, where each category is coded as a column of an indicator matrix or row and column of Burt matrix, it is not correct to simply analyse the corresponding submatrix of data, since the whole geometric structure is different for the submatrix . A simple modification of the correspondence analysis algorithm allows the overall geometric structure of the complete data set to be retained while calculating the solution for the selected subset of points. This strategy is useful for analysing patterns of response amongst any subset of categories and relating these patterns to demographic factors, especially for studying patterns of particular responses such as missing and neutral responses. The methodology is illustrated using data from the International Social Survey Program on Family and Changing Gender Roles in 1994.

Analysis of matched matrices

We consider the joint visualization of two matrices which have common rowsand columns, for example multivariate data observed at two time pointsor split accord-ing to a dichotomous variable. Methods of interest includeprincipal components analysis for interval-scaled data, or correspondenceanalysis for frequency data or ratio-scaled variables on commensuratescales. A simple result in matrix algebra shows that by setting up thematrices in a particular block format, matrix sum and difference componentscan be visualized. The case when we have more than two matrices is alsodiscussed and the methodology is applied to data from the InternationalSocial Survey Program.

From correspondence analysis to multiple and joint correspondence analysis

The generalization of simple (two-variable) correspondence analysis to more than two categorical variables, commonly referred to as multiple correspondence analysis, is neither obvious nor well-defined. We present two alternative ways of generalizing correspondence analysis, one based on the quantification of the variables and intercorrelation relationships, and the other based on the geometric ideas of simple correspondence analysis. We propose a version of multiple correspondence analysis, with adjusted principal inertias, as the method of choice for the geometric definition, since it contains simple correspondence analysis as an exact special case, which is not the situation of the standard generalizations. We also clarify the issue of supplementary point representation and the properties of joint correspondence analysis, a method that visualizes all two-way relationships between the variables. The methodology is illustrated using data on attitudes to science from the International Social Survey Program on Environment in 1993.

Achievement of at least very good partial response is a simple and robust prognostic factor in patients with multiple myeloma treated with high-dose therapy: long-term analysis of the IFM 99-02 and 99-04 Trials.

PURPOSE: The prognostic impact of complete response (CR) achievement in multiple myeloma (MM) has been shown mostly in the context of autologous stem-cell transplantation. Other levels of response have been defined because, even with high-dose therapy, CR is a relatively rare event. The purpose of this study was to analyze the prognostic impact of very good partial response (VGPR) in patients treated with high-dose therapy. PATIENTS AND METHODS: All patients were included in the Intergroupe Francophone du Myelome 99-02 and 99-04 trials and treated with vincristine, doxorubicin, and dexamethasone (VAD) induction therapy followed by double autologous stem-cell transplantation (ASCT). Best post-ASCT response assessment was available for 802 patients. RESULTS: With a median follow-up of 67 months, median event-free survival (EFS) and 5-year EFS were 42 months and 34%, respectively, for 405 patients who achieved at least VGPR after ASCT versus 32 months and 26% in 288 patients who achieved only partial remission (P = .005). Five-year overall survival (OS) was significantly superior in patients achieving at least VGPR (74% v 61% P = .0017). In multivariate analysis, achievement of less than VGPR was an independent factor predicting shorter EFS and OS. Response to VAD had no impact on EFS and OS. The impact of VGPR achievement on EFS and OS was significant in patients with International Staging System stages 2 to 3 and for patients with poor-risk cytogenetics t(4;14) or del(17p). CONCLUSION: In the context of ASCT, achievement of at least VGPR is a simple prognostic factor that has importance in intermediate and high-risk MM and can be informative in more patients than CR.

Mutational analysis of class A and class B penicillin-binding proteins in Streptococcus gordonii.

High-molecular-weight (HMW) penicillin-binding proteins (PBPs) are divided into class A and class B PBPs, which are bifunctional transpeptidases/transglycosylases and monofunctional transpeptidases, respectively. We determined the sequences for the HMW PBP genes of Streptococcus gordonii, a gingivo-dental commensal related to Streptococcus pneumoniae. Five HMW PBPs were identified, including three class A (PBPs 1A, 1B, and 2A) and two class B (PBPs 2B and 2X) PBPs, by homology with those of S. pneumoniae and by radiolabeling with [3H]penicillin. Single and double deletions of each of them were achieved by allelic replacement. All could be deleted, except for PBP 2X, which was essential. Morphological alterations occurred after deletion of PBP 1A (lozenge shape), PBP 2A (separation defect and chaining), and PBP 2B (aberrant septation and premature lysis) but not PBP 1B. The muropeptide cross-link patterns remained similar in all strains, indicating that cross-linkage for one missing PBP could be replaced by others. However, PBP 1A mutants presented shorter glycan chains (by 30%) and a relative decrease (25%) in one monomer stem peptide. Growth rate and viability under aeration, hyperosmolarity, and penicillin exposure were affected primarily in PBP 2B-deleted mutants. In contrast, chain-forming PBP 2A-deleted mutants withstood better aeration, probably because they formed clusters that impaired oxygen diffusion. Double deletion could be generated with any PBP combination and resulted in more-altered mutants. Thus, single deletion of four of the five HMW genes had a detectable effect on the bacterial morphology and/or physiology, and only PBP 1B seemed redundant a priori.

Injection drug use before and after liver transplantation: a retrospective multicenter analysis on incidence and outcome.

De Gottardi A, Hilleret M-N, Gelez P, La Mura V, Guillaud O, Majno P, Hadengue A, Morel P, Zarski J-P, Fontana M, Moradpour D, Mentha G, Boillot O, Leroy V, Giostra E, Dumortier J. Injection drug use before and after liver transplantation: a retrospective multicenter analysis on incidence and outcome. Clin Transplant 2009 DOI: 10.1111/j.1399-0012.2009.01121.x.Background and aims: Injecting drug use (IDU) before and after liver transplantation (LT) is poorly described. The aim of this study was to quantify relapse and survival in this population and to describe the causes of mortality after LT. Methods: Past injection drug users were identified from the LT listing protocols from four centers in Switzerland and France. Data on survival and relapse were collected and used for uni- and multivariate analysis. Results: Between 1988 and 2006, we identified 59 patients with a past history of IDU. The mean age at transplantation was 42.4 yr and the majority of patients were men (84.7%). The indication for LT was for the vast majority viral cirrhosis accounting for 91.5% of cases, while alcoholic cirrhosis was 5.1%. There were 16.9% of patients who had a substitution therapy before and 6.8% who continued after LT. Two patients (3.4%) relapsed into IDU after LT and died at 18 and 41 months. The mean follow-up was 51 months. Overall survival was 84%, 66%, and 61% at 1, 5, and 10 yr after transplantation. Conclusions: Documented IDU was rare in liver transplanted patients. Past IDU was not associated with poorer survival after LT, and relapse after LT occurred in 3.4%.

Neuropathological analysis of lacunes and microvascular lesions in late-onset depression.

M. Santos, G. Gold, E. Kövari, F. R. Herrmann, P. R. Hof, C. Bouras and P. Giannakopoulos (2010) Neuropathology and Applied Neurobiology36, 661-672
Neuropathological analysis of lacunes and microvascular lesions in late-onset depression Aims: Previous neuropathological studies documented that small vascular and microvascular pathology is associated with cognitive decline. More recently, we showed that thalamic and basal ganglia lacunes are associated with post-stroke depression and may affect emotional regulation. The present study examines whether this is also the case for late-onset depression. Methods: We performed a detailed analysis of small macrovascular and microvascular pathology in the post mortem brains of 38 patients with late-onset major depression (LOD) and 29 healthy elderly controls. A clinical diagnosis of LOD was established while the subjects were alive using the DSM-IV criteria. Additionally, we retrospectively reviewed all charts for the presence of clinical criteria of vascular depression. Neuropathological evaluation included bilateral semi-quantitative assessment of lacunes, deep white matter and periventricular demyelination, cortical microinfarcts and both focal and diffuse gliosis. The association between vascular burden and LOD was investigated using Fisher's exact test and univariate and multivariate logistic regression models. Results: Neither the existence of lacunes nor the presence of microvascular ischaemic lesions was related to occurrence of LOD. Similarly, there was no relationship between vascular lesion scores and LOD. This was also the case within the subgroup of LOD patients fulfilling the clinical criteria for vascular depression. Conclusions: Our results challenge the vascular depression hypothesis by showing that neither deep white matter nor periventricular demyelination is associated with LOD. In conjunction with our previous observations in stroke patients, they also imply that the impact of lacunes on mood may be significant solely in the presence of acute brain compromise.

Evidence against close linkage of the loci for fraXq of Martin-Bell syndrome and for factor IX.

Linkage between the loci for fraXq of Martin-Bell syndrome and factor IX was studied in nine families exhibiting this syndrome by means of a restriction fragment length polymorphism at the factor IX locus. Computer analysis of the data indicates there to be no evidence for close linkage between the syndrome and the factor IX locus.

Were Spanish migrants attracted by industrial agglomerations?: an analysis for the interwar years in the light of the new economic geography

In this paper we examine whether access to markets had a significant influence onmigration choices of Spanish internal migrants in the inter-war years. We perform astructural contrast of a New Economic Geography model that focus on the forwardlinkage that links workers location choice with the geography of industrial production,one of the centripetal forces that drive agglomeration in the NEG models. The resultshighlight the presence of this forward linkage in the Spanish economy of the inter-warperiod. That is, we prove the existence of a direct relation between workers¿ localizationdecisions and the market potential of the host regions. In addition, the direct estimationof the values associated with key parameters in the NEG model allows us to simulatethe migratory flows derived from different scenarios of the relative size of regions andthe distances between them. We show that in Spain the power of attraction of theagglomerations grew as they increased in size, but the high elasticity estimated for themigration costs reduced the intensity of the migratory flows. This could help to explainthe apparently low intensity of internal migrations in Spain until its upsurge during the1920s. This also explains the geography of migrations in Spain during this period,which hardly affected the regions furthest from the large industrial agglomerations (i.e.,regions such as Andalusia, Estremadura and Castile-La Mancha) but had an intenseeffect on the provinces nearest to the principal centres of industrial development.

Predictors of poor response to methotrexate in polyarticular-course juvenile idiopathic arthritis: analysis of the PRINTO methotrexate trial.

OBJECTIVES: To determine whether baseline demographic, clinical, articular and laboratory variables predict methotrexate (MTX) poor response in polyarticular-course juvenile idiopathic arthritis. METHODS: Patients newly treated for 6 months with MTX enrolled in the Paediatric Rheumatology International Trials Organization (PRINTO) MTX trial. Bivariate and logistic regression analyses were used to identify baseline predictors of poor response according to the American College of Rheumatology pediatric (ACR-ped) 30 and 70 criteria. RESULTS: In all, 405/563 (71.9%) of patients were women; median age at onset and disease duration were 4.3 and 1.4 years, respectively, with anti-nuclear antibody (ANA) detected in 259/537 (48.2%) patients. With multivariate logistic regression analysis, the most important determinants of ACR-ped 70 non-responders were: disease duration > 1.3 years (OR 1.93), ANA negativity (OR 1.77), Childhood Health Assessment Questionnaire (CHAQ) disability index > 1.125 (OR 1.65) and the presence of right and left wrist activity (OR 1.55). Predictors of ACR-ped 30 non-responders were: ANA negativity (OR 1.92), CHAQ disability index > 1.14 (OR 2.18) and a parent's evaluation of child's overall well-being < or = 4.69 (OR 2.2). CONCLUSION: The subgroup of patients with longer disease duration, ANA negativity, higher disability and presence of wrist activity were significantly associated with a poorer response to a 6-month MTX course.

Genetic association tests: a method for the joint analysis of family and case-control data.

With the trend in molecular epidemiology towards both genome-wide association studies and complex modelling, the need for large sample sizes to detect small effects and to allow for the estimation of many parameters within a model continues to increase. Unfortunately, most methods of association analysis have been restricted to either a family-based or a case-control design, resulting in the lack of synthesis of data from multiple studies. Transmission disequilibrium-type methods for detecting linkage disequilibrium from family data were developed as an effective way of preventing the detection of association due to population stratification. Because these methods condition on parental genotype, however, they have precluded the joint analysis of family and case-control data, although methods for case-control data may not protect against population stratification and do not allow for familial correlations. We present here an extension of a family-based association analysis method for continuous traits that will simultaneously test for, and if necessary control for, population stratification. We further extend this method to analyse binary traits (and therefore family and case-control data together) and accurately to estimate genetic effects in the population, even when using an ascertained family sample. Finally, we present the power of this binary extension for both family-only and joint family and case-control data, and demonstrate the accuracy of the association parameter and variance components in an ascertained family sample.

Clinical relevance of CD44 cell-surface expression and N-myc gene amplification in a multicentric analysis of 121 pediatric neuroblastomas.

PURPOSE: In contrast to other human tumors, a repression of the cell-surface glycoprotein CD44 on neuroblastoma is a marker of aggressiveness that usually correlates to N-myc amplification. We thus compared the prognostic value of both markers in the initial staging of 121 children treated for neuroblastoma in collaborative institutions. METHODS: Frozen samples were analyzed by a rapid and well-standardized technique of immunostaining with monoclonal antibodies (MoAbs) against epitopes in the CD44 constant region. RESULTS: In this retrospective series, CD44 was expressed on 102 specimens and strongly correlated with favorable tumor stages and histology, younger age, and normal N-myc copy numbers. In univariate analysis, CD44 expression and normal N-myc were the most powerful markers of favorable clinical outcome (P < 10(-6) and chi 2 = 65.40 and P < 10(-6) and chi 2 = 42.56, respectively), but analysis of CD44 affords significant prognostic discrimination in subgroups of patients with or without N-myc-amplified tumors. In the subgroup of stage IV neuroblastomas, CD44 was the only significant prognostic marker (P < .02, chi 2 = 5.76), whereas N-myc status was not discriminant. In multivariate analysis of five factors, ie, N-myc amplification, CD44 expression, age, tumor stage, and histology, the only independent prognostic factors of event-free survival were CD44 expression and tumor stage. CONCLUSION: The analysis of CD44 cell-surface expression must be recommended as an additional biologic marker in the initial staging of the disease.

Genome-wide meta-analysis for serum calcium identifies significantly associated SNPs near the calcium-sensing receptor (CASR) gene.

Calcium has a pivotal role in biological functions, and serum calcium levels have been associated with numerous disorders of bone and mineral metabolism, as well as with cardiovascular mortality. Here we report results from a genome-wide association study of serum calcium, integrating data from four independent cohorts including a total of 12,865 individuals of European and Indian Asian descent. Our meta-analysis shows that serum calcium is associated with SNPs in or near the calcium-sensing receptor (CASR) gene on 3q13. The top hit with a p-value of 6.3 x 10(-37) is rs1801725, a missense variant, explaining 1.26% of the variance in serum calcium. This SNP had the strongest association in individuals of European descent, while for individuals of Indian Asian descent the top hit was rs17251221 (p = 1.1 x 10(-21)), a SNP in strong linkage disequilibrium with rs1801725. The strongest locus in CASR was shown to replicate in an independent Icelandic cohort of 4,126 individuals (p = 1.02 x 10(-4)). This genome-wide meta-analysis shows that common CASR variants modulate serum calcium levels in the adult general population, which confirms previous results in some candidate gene studies of the CASR locus. This study highlights the key role of CASR in calcium regulation.