Empirical study of correlated survival times for recurrent events with proportional hazards margins and the effect of correlation and censoring.

Background: In longitudinal studies where subjects experience recurrent incidents over a period of time, such as respiratory infections, fever or diarrhea, statistical methods are required to take into account the within-subject correlation. Methods: For repeated events data with censored failure, the independent increment (AG), marginal (WLW) and conditional (PWP) models are three multiple failure models that generalize Cox"s proportional hazard model. In this paper, we revise the efficiency, accuracy and robustness of all three models under simulated scenarios with varying degrees of within-subject correlation, censoring levels, maximum number of possible recurrences and sample size. We also study the methods performance on a real dataset from a cohort study with bronchial obstruction. Results: We find substantial differences between methods and there is not an optimal method. AG and PWP seem to be preferable to WLW for low correlation levels but the situation reverts for high correlations. Conclusions: All methods are stable in front of censoring, worsen with increasing recurrence levels and share a bias problem which, among other consequences, makes asymptotic normal confidence intervals not fully reliable, although they are well developed theoretically.

Spatial statistics applied to soybean production data from Paraná State for 2003-04 to 2009-10 crop-years

In the current study, we performed a soybean production spatial distribution analysis in Paraná State. Seven crop-year data, from 2003-04 to 2009-10, obtained from the Paraná Department of Agriculture and Supply (SEAB) were used to develop a Boxmap for each crop-year, show soybean production throughout this time interval. Moran's index was used to measure spatial autocorrelation among municipalities at an aggregate level, while LISA index local correlation. For each index, different contiguity matrix and order were used and there was a significance level study. As a result, we have showed spatial relationship among cities regarding the production, which allowed the indication of high and low production clusters. Finally, identifying main soybean-producing cities, what may provide supply chain members with information to strengthen the crop production in Paraná.

Exact Tests for Contemporaneous Correlation of Disturbances in Seemingly Unrelated Regressions.

This paper proposes finite-sample procedures for testing the SURE specification in multi-equation regression models, i.e. whether the disturbances in different equations are contemporaneously uncorrelated or not. We apply the technique of Monte Carlo (MC) tests [Dwass (1957), Barnard (1963)] to obtain exact tests based on standard LR and LM zero correlation tests. We also suggest a MC quasi-LR (QLR) test based on feasible generalized least squares (FGLS). We show that the latter statistics are pivotal under the null, which provides the justification for applying MC tests. Furthermore, we extend the exact independence test proposed by Harvey and Phillips (1982) to the multi-equation framework. Specifically, we introduce several induced tests based on a set of simultaneous Harvey/Phillips-type tests and suggest a simulation-based solution to the associated combination problem. The properties of the proposed tests are studied in a Monte Carlo experiment which shows that standard asymptotic tests exhibit important size distortions, while MC tests achieve complete size control and display good power. Moreover, MC-QLR tests performed best in terms of power, a result of interest from the point of view of simulation-based tests. The power of the MC induced tests improves appreciably in comparison to standard Bonferroni tests and, in certain cases, outperforms the likelihood-based MC tests. The tests are applied to data used by Fischer (1993) to analyze the macroeconomic determinants of growth.

Statistics of avalanches in martensitic transformations, II. Modelling

Using a scaling assumption, we propose a phenomenological model aimed to describe the joint probability distribution of two magnitudes A and T characterizing the spatial and temporal scales of a set of avalanches. The model also describes the correlation function of a sequence of such avalanches. As an example we study the joint distribution of amplitudes and durations of the acoustic emission signals observed in martensitic transformations [Vives et al., preceding paper, Phys. Rev. B 52, 12 644 (1995)].

Modeling data with structural and temporal correlation using lower level and higher level multilevel models

Novel imaging techniques are playing an increasingly important role in drug development, providing insight into the mechanism of action of new chemical entities. The data sets obtained by these methods can be large with complex inter-relationships, but the most appropriate statistical analysis for handling this data is often uncertain - precisely because of the exploratory nature of the way the data are collected. We present an example from a clinical trial using magnetic resonance imaging to assess changes in atherosclerotic plaques following treatment with a tool compound with established clinical benefit. We compared two specific approaches to handle the correlations due to physical location and repeated measurements: two-level and four-level multilevel models. The two methods identified similar structural variables, but higher level multilevel models had the advantage of explaining a greater proportion of variation, and the modeling assumptions appeared to be better satisfied.

Theoretical insight into diagnosing observation error correlations using observation-minus-background and observation-minus-analysis statistics

To improve the quantity and impact of observations used in data assimilation it is necessary to take into account the full, potentially correlated, observation error statistics. A number of methods for estimating correlated observation errors exist, but a popular method is a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. The accuracy of the results it yields is unknown as the diagnostic is sensitive to the difference between the exact background and exact observation error covariances and those that are chosen for use within the assimilation. It has often been stated in the literature that the results using this diagnostic are only valid when the background and observation error correlation length scales are well separated. Here we develop new theory relating to the diagnostic. For observations on a 1D periodic domain we are able to the show the effect of changes in the assumed error statistics used in the assimilation on the estimated observation error covariance matrix. We also provide bounds for the estimated observation error variance and eigenvalues of the estimated observation error correlation matrix. We demonstrate that it is still possible to obtain useful results from the diagnostic when the background and observation error length scales are similar. In general, our results suggest that when correlated observation errors are treated as uncorrelated in the assimilation, the diagnostic will underestimate the correlation length scale. We support our theoretical results with simple illustrative examples. These results have potential use for interpreting the derived covariances estimated using an operational system.

Diagnosing observation error correlations for Doppler radar radial winds in the Met Office UKV model using observation-minus-background and observation-minus-analysis statistics

With the development of convection-permitting numerical weather prediction the efficient use of high resolution observations in data assimilation is becoming increasingly important. The operational assimilation of these observations, such as Dopplerradar radial winds, is now common, though to avoid violating the assumption of un- correlated observation errors the observation density is severely reduced. To improve the quantity of observations used and the impact that they have on the forecast will require the introduction of the full, potentially correlated, error statistics. In this work, observation error statistics are calculated for the Doppler radar radial winds that are assimilated into the Met Office high resolution UK model using a diagnostic that makes use of statistical averages of observation-minus-background and observation-minus-analysis residuals. This is the first in-depth study using the diagnostic to estimate both horizontal and along-beam correlated observation errors. By considering the new results obtained it is found that the Doppler radar radial wind error standard deviations are similar to those used operationally and increase as the observation height increases. Surprisingly the estimated observation error correlation length scales are longer than the operational thinning distance. They are dependent on both the height of the observation and on the distance of the observation away from the radar. Further tests show that the long correlations cannot be attributed to the use of superobservations or the background error covariance matrix used in the assimilation. The large horizontal correlation length scales are, however, in part, a result of using a simplified observation operator.

High-precision measurements of the co-polar correlation coefficient: non-Gaussian errors and retrieval of the dispersion parameter µ in rainfall

The co-polar correlation coefficient (ρhv) has many applications, including hydrometeor classification, ground clutter and melting layer identification, interpretation of ice microphysics and the retrieval of rain drop size distributions (DSDs). However, we currently lack the quantitative error estimates that are necessary if these applications are to be fully exploited. Previous error estimates of ρhv rely on knowledge of the unknown "true" ρhv and implicitly assume a Gaussian probability distribution function of ρhv samples. We show that frequency distributions of ρhv estimates are in fact highly negatively skewed. A new variable: L = -log10(1 - ρhv) is defined, which does have Gaussian error statistics, and a standard deviation depending only on the number of independent radar pulses. This is verified using observations of spherical drizzle drops, allowing, for the first time, the construction of rigorous confidence intervals in estimates of ρhv. In addition, we demonstrate how the imperfect co-location of the horizontal and vertical polarisation sample volumes may be accounted for. The possibility of using L to estimate the dispersion parameter (µ) in the gamma drop size distribution is investigated. We find that including drop oscillations is essential for this application, otherwise there could be biases in retrieved µ of up to ~8. Preliminary results in rainfall are presented. In a convective rain case study, our estimates show µ to be substantially larger than 0 (an exponential DSD). In this particular rain event, rain rate would be overestimated by up to 50% if a simple exponential DSD is assumed.

Inferences about the coefficient of correlation in the standard bivariate normal distribution

The study of the association between two random variables that have a joint normal distribution is of interest in applied statistics; for example, in statistical genetics. This article, targeted to applied statisticians, addresses inferences about the coefficient of correlation (ρ) in the bivariate normal and standard bivariate normal distributions using likelihood, frequentist, and Baycsian perspectives. Some results are surprising. For instance, the maximum likelihood estimator and the posterior distribution of ρ in the standard bivariate normal distribution do not follow directly from results for a general bivariate normal distribution. An example employing bootstrap and rejection sampling procedures is used to illustrate some of the peculiarities.

Correlation between oral health perception and clinical factors in a Brazilian community

Objective: The aim of this study was to analyse associations between self-perception of oral health and relevant clinical, personal and socio-demographic factors in a Brazilian community. Material And Methods: Urban adults living in a city in southern Brazil were interviewerd and examined. Individuals with acute pain and who needed multiple extractions of teeth were excluded. Self-perception and the Oral Health Impact Profile (OHIP-14) were applied to measure the impact of oral conditions on the quality of life. Socio-demographic and clinical indicators were also analysed. Results: The clinical examination revealed a high dental caries experience (DMFT = 18.9) and a high prevalence of periodontal disease. Oral condition was considered normal by 42% of respondents. The variables associated with the OHIP-14 were: education, age, self-assessment, dental caries and the DMFT index. Conclusions: Self-perception of oral health was associated with OHIP-14 and the clinical indicators had low influence in the self-perception. Therefore, the development of educational initiatives and preventive strategies for the adult population is recommended. © BASCD 2011.

Spatial correlation between the composition of the clay fraction and contents of available phosphorus of an Oxisol at hillslope scale

Relief is regarded as the abiotic factor most strongly influencing pedogenic processes at a local scale. The spatial correlations between the composition of the clay fraction (iron - Fe and aluminum - Al oxides, kaolinite and organic matter - OM) and contents of available phosphorus (P) of an Oxisol were evaluated at hillslope scale under sugarcane cultivation. A total of 119 samples were collected at intersection points on a 100. ×. 100. m georeferenced grid of regularly spaced points 10. m apart in the 0.2-0.4. m depth in an area consisting of two landform components namely: component I (an area with a linear hillslope curvature), and component II (one with a concave-convex hillslope curvature). Soil OM and available P contents were subjected to descriptive statistics and geostatistical analyses in order to assess their variability and spatial dependence. All attributes studied were spatially dependent. Available phosphorus had positive spatial correlation with high crystalline goethite, hematite and gibbsite. Identifying small hillslope curvatures is useful with a view to better understanding their relationships with soil organic matter and available phosphorus, as well as kaolinite and Fe and Al oxide attributes. A simple correlation analysis by itself is inadequate to relate attributes, which requires a supplemental, geostatistical technique. © 2012 Elsevier B.V..

Patients' expectation and satisfaction with removable dental prosthesis therapy and correlation with patients' evaluation of the dentists

Objective. Considering that patients' satisfaction is one of the most important goals in conducting removable dental prosthesis (RDP) therapy and the fact that there are many factors which influence this parameter, the present study aims to evaluate the expectation before and satisfaction after therapy with RDP in patients who seek such therapy. As a secondary objective, other variables that may be associated with patient satisfaction are also evaluated, such as gender, age, Kennedy's classification of the arch supporting the RDP, the number of RDP adjustments after delivery and patients' evaluation of the dentists' conduct. Materials and methods. A sample of 44 patients who received RDP therapy were assigned visual analog scale scores for their expectation before and satisfaction after therapy regarding chewing, aesthetics, comfort and phonetics. They also completed a questionnaire concerning the dentists' conduct. Results. There was no statistically significant difference among scores concerning different genders, age, number of post-delivery settings and arch involved in the RDP. Regarding patients' evaluation of the dentists' conduct, there was a predominance of positive evaluations, but only different answers to the statement (i.e. 'The dentist I saw thoroughly explained the recommended treatment before it commenced') present statistically significant different scores for chewing (p = 0.040) and phonetics (p = 0.046). Conclusions. The average visual analog scale scores were high for both expectation prior to treatment and satisfaction after treatment; however, the scores for expectations were higher than those for satisfaction. © 2013 Informa Healthcare.

Report on the production of trade statistics in the Caribbean region

Includes bibliography