RESUMEN

OBJECTIVES: Perceived pain is a multi-factorial subjective variable, commonly measured by numeric rating scales, verbal descriptive scales (VDS), or by a position on an analogue line (VAS). A major question is whether an individual's VAS and VDS pain assessments, on the same occasion, could be comparable. The aim was to compare continuous and discretized VAS pain data with verbal descriptive pain datasets from the Oswestry Disability Index (ODI) and the European Quality of Life Scale (EQ-5D) in paired pain datasets. METHODS: The measurement level of data from any type of scale assessments is ordinal, having rank-invariant properties only. Non-parametric statistical methods were used. Two ways of discretizing the VAS-line to VAS-intervals to fit the number of the comparing VDS-categories were used: the commonly used (equidistant VAS,VDS)-pairs and the (unbiased VAS,VDS)-pairs of pain data. The comparability of the (VAS,VDS)-pairs of data of perceived pain was studied by the bivariate ranking approach. Hence, each pair will be regarded as ordered, disordered, or tied with respect to the other pairs of data. The percentage agreement, PA, the measures of disorder, D, and of order consistency, MA, were calculated. Total interchangeability requires PA = 1 and MA = 1. RESULTS: The wide range of overlapping of (VAS,VDS)-pairs indicated that the continuous VAS data were not comparable to any of the VDS pain datasets. The percentage of agreement, PA; in the (equidistant VAS,ODI) and (equidistant VAS, EQ-5D) pairs were 38 and 49%, and the order consistency, MA, was 0.70 and 0.80, respectively. Corresponding results for the (unbiased VAS,VDS)-pairs of pain data were PA: 54 and 100%, and MA: 0.77 and 1.0. CONCLUSION: Our results confirmed that perceived pain is the individual's subjective experience, and possible scale-interchangeability is only study-specific. The pain experience is not possible to be measured univocally, but is possible for the individual to rate on a scale.

Asunto(s)

RESUMEN

BACKGROUND: Life expectancy is a simple measure of assessing health differences between two or more populations but current life expectancy calculations are not reliable for small populations. A potential solution to this is to borrow strength from larger populations from the same source, but this has not formally been investigated. METHODS: Using data on 451,222 individuals from the Clinical Practice Research Datalink on the presence/absence of intellectual disability and type 2 diabetes mellitus, we compared stratified and combined flexible parametric models, and Chiang's methods, for calculating life expectancy. Confidence intervals were calculated using the Delta method, Chiang's adjusted life table approach and bootstrapping. RESULTS: The flexible parametric models allowed calculation of life expectancy by exact age and beyond traditional life expectancy age thresholds. The combined model that fit age interaction effects as a spline term provided less bias and greater statistical precision for small covariate subgroups by borrowing strength from the larger subgroups. However, careful consideration of the distribution of events in the smallest group was needed. CONCLUSIONS: Life expectancy is a simple measure to compare health differences between populations. The use of combined flexible parametric methods to calculate life expectancy in small samples has shown promising results by allowing life expectancy to be modelled by exact age, greater statistical precision, less bias and prediction of different covariate patterns without stratification. We recommend further investigation of their application for both policymakers and researchers.

Asunto(s)

RESUMEN

The marine microalga Nannochloropsis oculata has garnered significant interest as a potential source of lipids, both for biofuel and nutrition, containing significant amounts of C16:0, C16:1, and C20:5, n-3 (EPA) fatty acids (FA). Growth parameters such as temperature, pH, light intensity, and nutrient availability play a crucial role in the fatty acid profile of microalgae, with N. oculata being no exception. This study aims to identify key variables for the FA profile of N. oculata grown autotrophically. To that end, the most relevant literature data were gathered and combined with our previous work as well as with novel experimental data, with 121 observations in total. The examined variables were the percentages of C14:0, C16:0, C16:1, C18:1, C18:2, and C20:5, n-3 in total FAs, their respective ratios to C16:0, and the respective content of biomass in those fatty acids in terms of ash free dry weight. Many potential predictor variables were collected, while dummy variables were introduced to account for bias in the measured variables originating from different authors as well as for other parameters. The method of multiple imputations was chosen to handle missing data, with limits based on the literature and model-based estimation, such as using the software PHREEQC and residual modelling for the estimation of pH. To eliminate unimportant predictor variables, LASSO (Least Absolute Shrinkage and Selection Operator) regression analysis with a novel definition of optimal lambda was employed. LASSO regression identified the most relevant predictors while minimizing the risk of overfitting the model. Subsequently, stepwise linear regression with interaction terms was used to further study the effects of the selected predictors. After two rounds of regression, sparse refined models were acquired, and their coefficients were evaluated based on significance. Our analysis confirms well-known effects, such as that of temperature, and it uncovers novel unreported effects of aeration, calcium, magnesium, and manganese. Of special interest is the negative effect of aeration on polyunsaturated fatty acids (PUFAs), which is possibly related to the enzymatic kinetics of fatty acid desaturation under increased oxygen concentration. These findings contribute to the optimization of the fatty acid profile of N. oculata for different purposes, such as production of, high in PUFAs, food or feed, or production of, high in saturated and monounsaturated FA methyl esters (FAME), biofuels.

RESUMEN

Research on the long-term effects of health in early life has predominantly relied on parametric methods to assess differences between groups of children. However, this approach leaves a wealth of distributional information untapped. The objective of this study was to assess distributional differences in earnings and mental health in young adulthood between individuals who suffered a chronic illness in childhood compared to those who did not using the non-parametric relative distributions framework. Using data from the Panel Study of Income Dynamics, we find that young adults who suffered a chronic illness in childhood fare worse in terms of earnings and mental health scores in adulthood, particularly for individuals reporting a childhood mental health/developmental disorder. Covariate decompositions suggest that chronic conditions in childhood may indirectly affect later outcomes through educational attainment: had the two groups had similar levels of educational attainment, the proportion of individuals with a report of a chronic condition in childhood in the lower decile of the relative earnings distribution would have been reduced by about 20 percentage points. Findings may inform policy aimed at mitigating longer run effects of health conditions in childhood and may generate hypotheses to be explored in parametric analyses.

RESUMEN

This study aimed to perform a genome-wide association analysis (GWAS) using the Random Forest (RF) approach for scanning candidate genes for age at first calving (AFC) in Nellore cattle. Additionally, potential epistatic effects were investigated using linear mixed models with pairwise interactions between all markers with high importance scores within the tree ensemble non-linear structure. Data from Nellore cattle were used, including records of animals born between 1984 and 2015 and raised in commercial herds located in different regions of Brazil. The estimated breeding values (EBV) were computed and used as the response variable in the genomic analyses. After quality control, the remaining number of animals and SNPs considered were 3,174 and 360,130, respectively. Five independent RF analyses were carried out, considering different initialization seeds. The importance score of each SNP was averaged across the independent RF analyses to rank the markers according to their predictive relevance. A total of 117 SNPs associated with AFC were identified, which spanned 10 autosomes (2, 3, 5, 10, 11, 17, 18, 21, 24, and 25). In total, 23 non-overlapping genomic regions embedded 262 candidate genes for AFC. Enrichment analysis and previous evidence in the literature revealed that many candidate genes annotated close to the lead SNPs have key roles in fertility, including embryo pre-implantation and development, embryonic viability, male germinal cell maturation, and pheromone recognition. Furthermore, some genomic regions previously associated with fertility and growth traits in Nellore cattle were also detected in the present study, reinforcing the effectiveness of RF for pre-screening candidate regions associated with complex traits. Complementary analyses revealed that many SNPs top-ranked in the RF-based GWAS did not present a strong marginal linear effect but are potentially involved in epistatic hotspots between genomic regions in different autosomes, remarkably in the BTAs 3, 5, 11, and 21. The reported results are expected to enhance the understanding of genetic mechanisms involved in the biological regulation of AFC in this cattle breed.

RESUMEN

This study aimed to introduce non-parametric tests and guard bands to assess the compliance of some river water properties with Brazilian environmental regulations. Due to the heterogeneity of the measurands pH, Biochemical Oxygen Demand (BOD), manganese molar concentration, and Escherichia coli, which could be wrongly treated as outliers, as well as the non-Gaussian data, robust methods were used to calculate the measurement uncertainty. Next, based on guard bands, the compliance assessment was evaluated using this previous uncertainty information. For these four measurands, partial overlaps between their uncertainties and the specification limit could generate doubts about compliance. The non-parametric approach for calculating the uncertainty connected to the guard bands concept classified pH and BOD as "conform", with a risk to the consumer of up to 4.0% and 4.9%, respectively; in contrast, manganese molar concentration and Escherichia coli were "not conform", with a risk to the consumer of up to 25% and 7.4%, respectively. The methodology proposed was satisfactory because it considered the natural heterogeneity of data with non-Gaussian behavior instead of wrongly excluding outliers. In an unprecedented way, two connected statistical approaches shed light on the measurement uncertainty in compliance assessment of water analysis.

Asunto(s)

RESUMEN

[11C]UCB-J is a PET radioligand that binds to the presynaptic vesicle glycoprotein 2A. Therefore, [11C]UCB-J PET may serve as an in vivo marker of synaptic integrity. The main objective of this study was to evaluate the quantitative accuracy and the 28-day test-retest repeatability (TRT) of various parametric quantitative methods for dynamic [11C]UCB-J studies in Alzheimer's disease (AD) patients and healthy controls (HC). Eight HCs and seven AD patients underwent two 60-min dynamic [11C]UCB-J PET scans with arterial sampling over a 28-day interval. Several plasma-input based and reference-region based parametric methods were used to generate parametric images using metabolite corrected plasma activity as input function or white matter semi-ovale as reference region. Different parametric outcomes were compared regionally with corresponding non-linear regression (NLR) estimates. Furthermore, the 28-day TRT was assessed for all parametric methods. Spectral analysis (SA) and Logan graphical analysis showed high correlations with NLR estimates. Receptor parametric mapping (RPM) and simplified reference tissue model 2 (SRTM2) BPND, and reference Logan (RLogan) distribution volume ratio (DVR) regional estimates correlated well with plasma-input derived DVR and SRTM BPND. Among the multilinear reference tissue model (MRTM) methods, MRTM1 had the best correspondence with DVR and SRTM BPND. Among the parametric methods evaluated, spectral analysis (SA) and SRTM2 were the best plasma-input and reference tissue methods, respectively, to obtain quantitatively accurate and repeatable parametric images for dynamic [11C]UCB-J PET.

RESUMEN

The same intervention can produce different effects in different sites. Existing transport mediation estimators can estimate the extent to which such differences can be explained by differences in compositional factors and the mechanisms by which mediating or intermediate variables are produced; however, they are limited to consider a single, binary mediator. We propose novel nonparametric estimators of transported interventional (in)direct effects that consider multiple, high-dimensional mediators and a single, binary intermediate variable. They are multiply robust, efficient, asymptotically normal, and can incorporate data-adaptive estimation of nuisance parameters. They can be applied to understand differences in treatment effects across sites and/or to predict treatment effects in a target site based on outcome data in source sites.

Asunto(s)

RESUMEN

The provision of high-quality food is a primary factor in ensuring adequate nourishment and preventing malnourishment-related diseases in Pakistan. This study, therefore, aimed to quantify the impact of income on nutrient consumption in Pakistan, with the hypothesis that income has a primary role in reducing malnourishment in the developing world. To do this, we estimated nutrient-income elasticity-defined as the proportion of change in nutrient consumption in response to a change in income-for total calories, macronutrients, and micronutrients, using the nationally representative Household Integrated Economic Survey data (2010-2011) for Pakistan. Nutrient-income elasticity values were derived using several parametric regression approaches. We also assessed the non-linearity and endogeneity of the relationship. Calorie-income elasticity was found to be significantly different from zero, irrespective of the estimation technique used. Income elasticity for macronutrients and micronutrients was also found to be significantly different from zero, ranging from 0.29 to 0.65. This study, therefore, supports the hypothesis that increased household income likely improves nutrient consumption.

RESUMEN

Sustainable development can be achieved through integrations of livestock and crops farm levels that generate higher economic efficiency in saving production costs through combining of crop and livestock. So, this study analyzed the allocative and economic efficiency of smallholder farmers in mixed crop-livestock production in Ethiopia evidence from Horo district using cross sectional data collected from 152 households in 2019/2020 cropping season by using structured questionnaires'. Data envelopment analysis was used to estimate the allocative and economic efficiency score. From the results of data envelopment analysis mean allocative efficiency was (57.0%) and economic efficiency (38.4%). Thus existence of about 61.6% economic inefficiency in the production of mixed crop-livestock production. Moreover, Tobit regression model results show that allocative efficiency affected by extension, off-non-farm and education levels of household positively. While economic efficiency positively affected by credit use, terrace and extension service positively and distance to the market negatively. Hence government should take steps for the improvement in education levels of household, development of terrace, expansion of nonfarm sectors and reform with extension service.

RESUMEN

PET imaging is an important diagnostic tool for management of patients with cancer and other diseases. Medical decisions based on quantitative PET information could potentially benefit from the availability of tools for evaluation of associated uncertainties. Raw PET data can be viewed as a sample from an inhomogeneous Poisson process so there is the possibility to directly apply bootstrapping to raw projection-domain list-mode data. Unfortunately this is computationally impractical, particularly if data reconstruction is iterative or the acquisition protocol is dynamic. We develop a flexible statistical linear model analysis to be used with multi-frame PET image data to create valid bootstrap samples. The technique is illustrated using data from dynamic PET studies with fluoro-deoxyglucose (FDG) and fluoro-thymidine (FLT) in brain and breast cancer patients. As is often the case with dynamic PET studies, data have been archived without raw list-mode information. Using the bootstrapping technique maps of kinetic parameters and associated uncertainties are obtained. The quantitative performance of the approach is assessed by simulation. The proposed image-domain bootstrap is found to substantially match the projection-domain alternative. Analysis of results points to a close relation between relative uncertainty in voxel-level kinetic parameters and local reconstruction error. This is consistent with statistical theory.

Asunto(s)

RESUMEN

BACKGROUND: Asthma is one of the leading chronic illnesses among children in the United States. Asthma prevalence is higher among African Americans (11.2%) compared to European Americans (7.7%). Bronchodilator medications are part of the first-line therapy, and the rescue medication, for acute asthma symptoms. Bronchodilator drug response (BDR) varies substantially among different racial/ethnic groups. Asthma prevalence in African Americans is only 3.5% higher than that of European Americans, however, asthma mortality among African Americans is four times that of European Americans; variation in BDR may play an important role in explaining this health disparity. To improve our understanding of disparate health outcomes in complex phenotypes such as BDR, it is important to consider interactions between environmental and biological variables. RESULTS: We evaluated the impact of pairwise and three-variable interactions between environmental, social, and biological variables on BDR in 233 African American youth with asthma using Visualization of Statistical Epistasis Networks (ViSEN). ViSEN is a non-parametric entropy-based approach able to quantify interaction effects using an information-theory metric known as Information Gain (IG). We performed analyses in the full dataset and in sex-stratified subsets. Our analyses identified several interaction models significantly, and suggestively, associated with BDR. The strongest interaction significantly associated with BDR was a pairwise interaction between pre-natal smoke exposure and socioeconomic status (full dataset IG: 2.78%, p = 0.001; female IG: 7.27%, p = 0.004)). Sex-stratified analyses yielded divergent results for females and males, indicating the presence of sex-specific effects. CONCLUSIONS: Our study identified novel interaction effects significantly, and suggestively, associated with BDR in African American children with asthma. Notably, we found that all of the interactions identified by ViSEN were "pure" interaction effects, in that they were not the result of strong main effects on BDR, highlighting the complexity of the network of biological and environmental factors impacting this phenotype. Several associations uncovered by ViSEN would not have been detected using regression-based methods, thus emphasizing the importance of employing statistical methods optimized to detect both additive and non-additive interaction effects when studying complex phenotypes such as BDR. The information gained in this study increases our understanding and appreciation of the complex nature of the interactions between environmental and health-related factors that influence BDR and will be invaluable to biomedical researchers designing future studies.

RESUMEN

The aim of this study was to investigate the test-retest (TRT) repeatability of various parametric quantification methods for [18F]Flortaucipir positron emission tomography (PET). We included eight subjects with dementia or mild cognitive impairment due to Alzheimer's disease and six cognitively normal subjects. All underwent two 130-min dynamic [18F]Flortaucipir PET scans within 3 ± 1 weeks. Data were analyzed using reference region models receptor parametric mapping (RPM), simplified reference tissue method 2 (SRTM2) and reference logan (RLogan), as well as standardized uptake value ratios (SUVr, time intervals 40-60, 80-100 and 110-130 min post-injection) with cerebellar gray matter as reference region. We obtained distribution volume ratio or SUVr, first for all brain regions and then in three tau-specific regions-of-interest (ROIs). TRT repeatability (%) was defined as |retest-test|/(average (test + retest)) × 100. For all methods and across ROIs, TRT repeatability ranged from (median (IQR)) 0.84% (0.68-2.15) to 6.84% (2.99-11.50). TRT repeatability was good for all reference methods used, although semi-quantitative models (i.e. SUVr) performed marginally worse than quantitative models, for instance TRT repeatability of RPM: 1.98% (0.78-3.58) vs. SUVr80-100: 3.05% (1.28-5.52), p < 0.001. Furthermore, for SUVr80-100 and SUVr110-130, with higher average SUVr, more variation was observed. In conclusion, while TRT repeatability was good for all models used, quantitative methods performed slightly better than semi-quantitative methods.

Asunto(s)

RESUMEN

It is estimated that a regular reader of biomedical publications will have access to 70% of journals with bivariate procedures with basic statistics, that is, any health professional with basic knowledge is capable of understanding the vast majority of published studies. That is why the need arises to write this text that aims to give an overview of the statistics applied to surgery, making reference to how and when to use each statistical procedure and of course how to interpret the results. We will handle definitions of concepts such as variables, hypotheses, confidence intervals, multivariate analysis, contrast of hypotheses, parametric and nonparametric methods and logistic regression. The importance of this guide is that all health professionals can interpret in a simple and concrete way the different statistical studies in their professional development.

Se estima que un lector habitual de publicaciones biomédicas tendrá acceso a un 70% de revistas con procedimientos bivariados con estadística básica, es decir, cualquier profesional de la salud con conocimientos básicos es capaz de entender la gran mayoría de los estudios publicados. Por eso surge la necesidad de escribir este texto que pretende dar una visión general de la estadística aplicada a la cirugía, haciendo referencia a cómo y cuándo utilizar cada procedimiento estadístico y, por supuesto, cómo interpretar los resultados. Manejaremos definiciones de conceptos como variables, hipótesis, intervalos de confianza, análisis multivariados, contraste de hipótesis, pruebas paramétricas y no paramétricas, y regresión logística. La importancia de esta guía es que todos los profesionales de la salud puedan interpretar de manera sencilla y concreta los diferentes estudios estadísticos en su desarrollo profesional.

Asunto(s)

RESUMEN

La diversidad genómica, expresada en las diferencias entre haplotipos moleculares de un conjunto de individuos, puede dividirse en componentes de variabilidad entre y dentro de algún factor de clasificación de los individuos. Para tal partición de varianzas, se usa análisis molecular de la varianza (AMOVA), el cual se construye a partir de las distancias multivariadas entre pares de haplotipos. El AMOVA clásico permite evaluar la significancia estadística de dos o más factores jerárquicos y consecuentemente no existe prueba de interacción entre factores. Sin embargo, existen situaciones donde los factores que clasifican a los individuos están cruzados y no anidados, es decir todos los niveles de un factor se encuentran representados en cada nivel del otro factor. Este trabajo propone una prueba estadística para evaluar la interacción entre factores cruzados en un AMOVA No-Jerárquico. La hipótesis nula de interacción establece que las diferencias moleculares entre individuos de distintos niveles de un factor son las mismas para todos los niveles del otro factor que los clasifica. La propuesta de análisis de interacción de factores a partir de distancias en un AMOVA No-Jerárquico comprende: cálculo de la matriz de distancia y partición de la misma en bloques, posterior cálculo de residuos y análisis de varianza no-paramétrico sobre los residuos. Su implementación es ilustrada en escenarios simulados y real. Los resultados sugieren que la prueba de interacción propuesta para el AMOVA No- Jerárquico presenta alta potencia.

The genomic diversity, expressed in the differences between molecular haplotypes of a group of individuals, can be divided into components of variability between and within some factor of classification of the individuals. For such variance partitioning, molecular analysis of variance (AMOVA) is used, which is constructed from the multivariate distances between pairs of haplotypes. The classical AMOVA allows the evaluation of the statistical significance of two or more hierarchical factors and consequently there is no interaction test between factors. However, there are situations where the factors that classify individuals are crossed rather than nested, that is, all the levels of a factor are represented in each level of the other one. This paper proposes a statistical test to evaluate the interaction between crossed factors in a Non-Hierarchical AMOVA. The null hypothesis of interaction establishes that the molecular differences between individuals of different levels of a factor are the same for all the levels of the other factor that classifies them. The proposed analysis of interaction in a Non- Hierarchical AMOVA includes: calculation of the distance matrix and partition of it into blocks, subsequent calculation of residuals and analysis of non-parametric variance on the residuals. Its implementation is illustrated in simulated and real scenarios. The results suggest that the proposed interaction test for the Non-Hierarchical AMOVA presents high power.

RESUMEN

BACKGROUND: Machine learning methods have gained popularity and practicality in identifying linear and non-linear effects of variants associated with complex disease/traits. Detection of epistatic interactions still remains a challenge due to the large number of features and relatively small sample size as input, thus leading to the so-called "short fat data" problem. The efficiency of machine learning methods can be increased by limiting the number of input features. Thus, it is very important to perform variable selection before searching for epistasis. Many methods have been evaluated and proposed to perform feature selection, but no single method works best in all scenarios. We demonstrate this by conducting two separate simulation analyses to evaluate the proposed collective feature selection approach. RESULTS: Through our simulation study we propose a collective feature selection approach to select features that are in the "union" of the best performing methods. We explored various parametric, non-parametric, and data mining approaches to perform feature selection. We choose our top performing methods to select the union of the resulting variables based on a user-defined percentage of variants selected from each method to take to downstream analysis. Our simulation analysis shows that non-parametric data mining approaches, such as MDR, may work best under one simulation criteria for the high effect size (penetrance) datasets, while non-parametric methods designed for feature selection, such as Ranger and Gradient boosting, work best under other simulation criteria. Thus, using a collective approach proves to be more beneficial for selecting variables with epistatic effects also in low effect size datasets and different genetic architectures. Following this, we applied our proposed collective feature selection approach to select the top 1% of variables to identify potential interacting variables associated with Body Mass Index (BMI) in ~ 44,000 samples obtained from Geisinger's MyCode Community Health Initiative (on behalf of DiscovEHR collaboration). CONCLUSIONS: In this study, we were able to show that selecting variables using a collective feature selection approach could help in selecting true positive epistatic variables more frequently than applying any single method for feature selection via simulation studies. We were able to demonstrate the effectiveness of collective feature selection along with a comparison of many methods in our simulation analysis. We also applied our method to identify non-linear networks associated with obesity.

RESUMEN

This review article presents an assessment of trends in statistical methods and an evaluation of their appropriateness in articles published in the Archives of Plastic Surgery (APS) from 2012 to 2017. We reviewed 388 original articles published in APS between 2012 and 2017. We categorized the articles that used statistical methods according to the type of statistical method, the number of statistical methods, and the type of statistical software used. We checked whether there were errors in the description of statistical methods and results. A total of 230 articles (59.3%) published in APS between 2012 and 2017 used one or more statistical method. Within these articles, there were 261 applications of statistical methods with continuous or ordinal outcomes, and 139 applications of statistical methods with categorical outcome. The Pearson chi-square test (17.4%) and the Mann-Whitney U test (14.4%) were the most frequently used methods. Errors in describing statistical methods and results were found in 133 of the 230 articles (57.8%). Inadequate description of P-values was the most common error (39.1%). Among the 230 articles that used statistical methods, 71.7% provided details about the statistical software programs used for the analyses. SPSS was predominantly used in the articles that presented statistical analyses. We found that the use of statistical methods in APS has increased over the last 6 years. It seems that researchers have been paying more attention to the proper use of statistics in recent years. It is expected that these positive trends will continue in APS.

RESUMEN

In a recent study Eklund et al. have shown that cluster-wise family-wise error (FWE) rate-corrected inferences made in parametric statistical method-based functional magnetic resonance imaging (fMRI) studies over the past couple of decades may have been invalid, particularly for cluster defining thresholds less stringent than p < 0.001; principally because the spatial autocorrelation functions (sACFs) of fMRI data had been modeled incorrectly to follow a Gaussian form, whereas empirical data suggest otherwise. Hence, the residuals from general linear model (GLM)-based fMRI activation estimates in these studies may not have possessed a homogenously Gaussian sACF. Here we propose a method based on the assumption that heterogeneity and non-Gaussianity of the sACF of the first-level GLM analysis residuals, as well as temporal autocorrelations in the first-level voxel residual time-series, are caused by unmodeled MRI signal from neuronal and physiological processes as well as motion and other artifacts, which can be approximated by appropriate decompositions of the first-level residuals with principal component analysis (PCA), and removed. We show that application of this method yields GLM residuals with significantly reduced spatial correlation, nearly Gaussian sACF and uniform spatial smoothness across the brain, thereby allowing valid cluster-based FWE-corrected inferences based on assumption of Gaussian spatial noise. We further show that application of this method renders the voxel time-series of first-level GLM residuals independent, and identically distributed across time (which is a necessary condition for appropriate voxel-level GLM inference), without having to fit ad hoc stochastic colored noise models. Furthermore, the detection power of individual subject brain activation analysis is enhanced. This method will be especially useful for case studies, which rely on first-level GLM analysis inferences.

Asunto(s)

RESUMEN

In a recent study, Eklund et al. employed resting-state functional magnetic resonance imaging data as a surrogate for null functional magnetic resonance imaging (fMRI) datasets and posited that cluster-wise family-wise error (FWE) rate-corrected inferences made by using parametric statistical methods in fMRI studies over the past two decades may have been invalid, particularly for cluster defining thresholds less stringent than p < 0.001; this was principally because the spatial autocorrelation functions (sACF) of fMRI data had been modeled incorrectly to follow a Gaussian form, whereas empirical data suggested otherwise. Here, we show that accounting for non-Gaussian signal components such as those arising from resting-state neural activity as well as physiological responses and motion artifacts in the null fMRI datasets yields first- and second-level general linear model analysis residuals with nearly uniform and Gaussian sACF. Further comparison with nonparametric permutation tests indicates that cluster-based FWE corrected inferences made with Gaussian spatial noise approximations are valid.

Asunto(s)

RESUMEN

This review article presents an assessment of trends in statistical methods and an evaluation of their appropriateness in articles published in the Archives of Plastic Surgery (APS) from 2012 to 2017. We reviewed 388 original articles published in APS between 2012 and 2017. We categorized the articles that used statistical methods according to the type of statistical method, the number of statistical methods, and the type of statistical software used. We checked whether there were errors in the description of statistical methods and results. A total of 230 articles (59.3%) published in APS between 2012 and 2017 used one or more statistical method. Within these articles, there were 261 applications of statistical methods with continuous or ordinal outcomes, and 139 applications of statistical methods with categorical outcome. The Pearson chi-square test (17.4%) and the Mann-Whitney U test (14.4%) were the most frequently used methods. Errors in describing statistical methods and results were found in 133 of the 230 articles (57.8%). Inadequate description of P-values was the most common error (39.1%). Among the 230 articles that used statistical methods, 71.7% provided details about the statistical software programs used for the analyses. SPSS was predominantly used in the articles that presented statistical analyses. We found that the use of statistical methods in APS has increased over the last 6 years. It seems that researchers have been paying more attention to the proper use of statistics in recent years. It is expected that these positive trends will continue in APS.

Asunto(s)