Modeling Motor Learning Using Heteroscedastic Functional Principal Components Analysis

Total Page:16

File Type:pdf, Size:1020Kb

Modeling Motor Learning Using Heteroscedastic Functional Principal Components Analysis Journal of the American Statistical Association ISSN: 0162-1459 (Print) 1537-274X (Online) Journal homepage: https://www.tandfonline.com/loi/uasa20 Modeling Motor Learning Using Heteroscedastic Functional Principal Components Analysis Daniel Backenroth, Jeff Goldsmith, Michelle D. Harran, Juan C. Cortes, John W. Krakauer & Tomoko Kitago To cite this article: Daniel Backenroth, Jeff Goldsmith, Michelle D. Harran, Juan C. Cortes, John W. Krakauer & Tomoko Kitago (2018) Modeling Motor Learning Using Heteroscedastic Functional Principal Components Analysis, Journal of the American Statistical Association, 113:523, 1003-1015, DOI: 10.1080/01621459.2017.1379403 To link to this article: https://doi.org/10.1080/01621459.2017.1379403 View supplementary material Accepted author version posted online: 29 Sep 2017. Published online: 29 Sep 2017. Submit your article to this journal Article views: 513 View Crossmark data Full Terms & Conditions of access and use can be found at https://www.tandfonline.com/action/journalInformation?journalCode=uasa20 JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION , VOL. , NO. , –: Applications and Case Studies https://doi.org/./.. Modeling Motor Learning Using Heteroscedastic Functional Principal Components Analysis Daniel Backenrotha, Jeff Goldsmitha, Michelle D. Harranb,c,JuanC.Cortesb,d, John W. Krakauerb,c,d, and Tomoko Kitagob,e aDepartment of Biostatistics, Mailman School of Public Health, Columbia University, New York, NY; bDepartment of Neurology, Columbia University Medical Center, New York, NY; cDepartment of Neuroscience, Johns Hopkins University School of Medicine, Baltimore, MD; dDepartment of Neurology, Johns Hopkins University School of Medicine, Baltimore, MD; eBurke Neurological Institute, White Plains, NY ABSTRACT ARTICLE HISTORY We propose a novel method for estimating population-level and subject-specific effects of covariates on the Received February variability of functional data. We extend the functional principal components analysis framework by mod- Revised August eling the variance of principal component scores as a function of covariates and subject-specific random KEYWORDS effects. In a setting where principal components are largely invariant across subjects and covariate values, Functional data; Kinematic modeling the variance of these scores provides a flexible and interpretable way to explore factors that affect data; Motor control; the variability of functional data. Our work is motivated by a novel dataset from an experiment assessing Probabilistic PCA; Variance upper extremity motor control, and quantifies the reduction in movement variability associated with skill modeling; Variational Bayes learning. The proposed methods can be applied broadly to understand movement variability, in settings that include motor learning, impairment due to injury or disease, and recovery. Supplementary materials for this article are available online. 1. Scientific Motivation variance and recovery (Krakauer 2006). By taking full advantage of high-frequency laboratory recordings, we shift focus from 1.1. Motor Learning particular time points to complete curves. Our approach allows us to model the variability of these curves as they depend on Recent work in motor learning has suggested that change in covariates,likethehandusedortherepetitionnumber,aswell movement variability is an important component of improve- as the estimation of random effects reflecting differences in base- ment in motor skill. It has been suggested that when a motor line variability and learning rates among subjects. task is learned, variance is reduced along dimensions relevant Section 1.2 describes our motivating data in more detail, and to the successful accomplishment of the task, although it may Section 2 introduces our modeling framework. A review of rele- increase in other dimensions (Scholz and Schöner 1999;Yarrow, vantstatisticalworkappearsinSection 3.Detailsofourestima- Brown, and Krakauer 2009). Experimental work, moreover, tion approach are in Section 4. Simulations and the application has shown that learning-induced improvement of movement to our motivating data appear in Sections 5 and 6,respectively, execution, measured through the trade-off between speed and andweclosewithadiscussioninSection 7. accuracy, is accompanied by significant reductions in movement variability. In fact, these reductions in movement variability may be a more important feature of learning than changes in the 1.2. Dataset average movement (Shmuelof, Krakauer, and Mazzoni 2012). These results have typically been based on assessments of vari- Our motivating data were gathered as part of a study of motor ability at a few time points, for example, at the end of the move- learning among healthy subjects. Kinematic data were acquired ment, although high-frequency laboratory recordings of com- in a standard task used to measure control of reaching move- plete movements are often available. ments. In this task, subjects rest their forearm on an air-sled In this article we develop a modeling framework that can be system to reduce effects of friction and gravity. The subjects are used to quantify movement variability based on dense record- presented with a screen showing eight targets arranged in a cir- ings of fingertip position throughout movement execution. This cle around a starting point, and they reach with their arm to a framework can be used to explore many aspects of motor skill target and back when it is illuminated on the screen. Subjects’ and learning: differences in baseline skill among healthy sub- movements are displayed on the screen, and they are rewarded jects, effects of repetition and training to modulate variability with 10 points if they turn their hand around within the target, over time, or the effect of baseline stroke severity on movement and 3 or 1 otherwise, depending on how far their hand is from CONTACT Daniel Backenroth [email protected] Department of Biostatistics, Mailman School of Public Health, Columbia University, West th Street, th Floor, New York, NY . Color versions of one or more of the figures in the article can be found online at www.tandfonline.com/r/JASA. Supplementary materials for this article are available online. Please go to www.tandfonline.com/r/JASA. © American Statistical Association 1004 D. BACKENROTH ET AL. thetargetatthepointofreturn.Subjectsarenotrewardedfor each target by each subject. The reduction in movement vari- movements outside prespecified velocity thresholds. ability after practice is clear. Our dataset consists of 9481 movements by 26 right-handed Prior to our analyses, we rotate curves so that all movements subjects. After becoming familiarized with the experimental extend to the target at 0◦. This rotation preserves shape and apparatus, each subject made 24 or 25 reaching movements to scale, but improves interpretation. After rotation, movement each of the 8 targets, in a semirandom order, with both the left along the X coordinate represents movement parallel to the line and right hand. Movements that did not reach at least 30% of the between origin and target, and movement along the Y coordi- distance to the target and movements with a direction more than nate represents movement perpendicular to this line. We build 90 deg away from the target direction at the point of peak veloc- models for X and Y coordinate curves separately in our primary ity were excluded from the dataset, because of the likelihood analysis. An alternative bivariate analysis appears in Appendix that they were made to the wrong target or not attempted due C. to distraction. Movements made at speeds outside the range of interest, with peak velocity less than 0.04 or greater than 2.0 m/s, 2. Model for Trajectory Variability were also excluded. These exclusion rules and other similar rules We adopt a functional data approach to model position curves have been used previously in similar kinematic experiments, and P (t).HereweomittheX andY superscripts for notational sim- are designed to increase the specificity of these experiments for ij plicity. Our starting point is the functional principal component probingmotorcontrolmechanisms(Huangetal.2012; Tanaka, analysis (FPCA) model of Yao, Müller, and Wang (2005)with Sejnowski, and Krakauer 2009;Kitagoetal.2015). A small num- subject-specific means. In this model, it is assumed that each ber of additional movements were removed from the dataset due curve P (t) can be modeled as to instrumentation and recording errors. The data we consider ij have not been previously reported. ( ) = μ ( ) + δ ( ) For each movement, the X and Y position of the hand move- Pij t ij t ij t ∞ ment is recorded as a function of time from movement onset to = μ ( ) + ξ φ ( ) + ( ). the initiation of return to the starting point, resulting in bivari- ij t ijk k t ij t (1) X ( ), Y ( ) k=1 ate functional observations denoted [Pij t Pij t ]forsubject i and movement j. In practice, observations are recorded not Here μij(t) is the mean function for curve Pij(t),thedevia- as functions but as discrete vectors. There is some variability in tion δij(t) is modeled as a linear combination of eigenfunctions φ ( ) ξ movement duration, which we remove for computational con- k t ,the ijk are uncorrelated random variables with mean 0 venience by linearly registering each movement onto a common and variances λk,where λk < ∞ and λ1 ≥ λ2 ≥··· , and = k grid of length D 50. The structure of the registered kinematic ij(t) iswhitenoise.Hereallthedeviationsδij(t) are assumed data is illustrated
Recommended publications
  • T-TEST Outline Hypothesis Testing Steps Bivariate Analysis
    Bivariate Analysis T-TEST Variable 1 2 LEVELS >2 LEVELS CONTINUOUS Variable 2 2 LEVELS X2 X2 t-test chi square test chi square test >2 LEVELS X2 X2 ANOVA chi square test chi square test (F-test) CONTINUOUS t-test ANOVA -Correlation (F-test) -Simple linear Regression Outline Comparison of means: t-test Hypothesis testing steps T-test is used when one variable is of a continuous nature and the other is dichotomous. T-test The t-test is used to compare the means of two groups on a given variable. Anova Examples: Difference in average blood pressure among males & females. Difference in average BMI among those who exercise and those who do not. Hypothesis testing steps Comparison of means: t-test Example 1: Identify the study objective Research question: Among university students, is there a State the null & alternative hypothesis difference between the average weight for males versus females? Select the proper test statistic Null hypothesis (Ho): μ weight males = μ weight females Calculate the test statistic Alternative hypothesis (Ha): μ weight males ≠ μ weight females Take a statistical decision based on the p-value. Statistical test: t-test Reject or accept the null hypothesis Comparison of means: t-test Comparison of means: t-test T-Test (SPSS output) If this p-value is < 0.05 then reject null hypothesis and conclude that the variances are different (accept alternative) Group Statistics and hence check this p-value for the t-test. Std. Error gender N Mean Std. Deviation Mean weight male 804 75.92 12.843 .453 If this p-value is > 0.05 then accept null hypothesis and female 1135 56.47 8.923 .265 conclude that the variances are equal and hence check this p-value for the t-test.
    [Show full text]
  • Chapter Five Multivariate Statistics
    Chapter Five Multivariate Statistics Jorge Luis Romeu IIT Research Institute Rome, NY 13440 April 23, 1999 Introduction In this chapter we treat the multivariate analysis problem, which occurs when there is more than one piece of information from each subject, and present and discuss several materials analysis real data sets. We first discuss several statistical procedures for the bivariate case: contingency tables, covariance, correlation and linear regression. They occur when both variables are either qualitative or quantitative: Then, we discuss the case when one variable is qualitative and the other quantitative, via the one way ANOVA. We then overview the general multivariate regression problem. Finally, the non-parametric case for comparison of several groups is discussed. We emphasize the assessment of all model assumptions, prior to model acceptance and use and we present some methods of detection and correction of several types of assumption violation problems. The Case of Bivariate Data Up to now, we have dealt with data sets where each observation consists of a single measurement (e.g. each observation consists of a tensile strength measurement). These are called univariate observations and the related statistical problem is known as univariate analysis. In many cases, however, each observation yields more than one piece of information (e.g. tensile strength, material thickness, surface damage). These are called multivariate observations and the statistical problem is now called multivariate analysis. Multivariate analysis is of great importance and can help us enhance our data analysis in several ways. For, coming from the same subject, multivariate measurements are often associated with each other. If we are able to model this association then we can take advantage of the situation to obtain one from the other.
    [Show full text]
  • Package 'Biwavelet'
    Package ‘biwavelet’ May 26, 2021 Type Package Title Conduct Univariate and Bivariate Wavelet Analyses Version 0.20.21 Date 2021-05-24 Author Tarik C. Gouhier, Aslak Grinsted, Viliam Simko Maintainer Tarik C. Gouhier <[email protected]> Description This is a port of the WTC MATLAB package written by Aslak Grinsted and the wavelet program written by Christopher Torrence and Gibert P. Compo. This package can be used to perform univariate and bivariate (cross-wavelet, wavelet coherence, wavelet clustering) analyses. License GPL (>= 2) URL https://github.com/tgouhier/biwavelet BugReports https://github.com/tgouhier/biwavelet/issues LazyData yes LinkingTo Rcpp Imports fields, foreach, Rcpp (>= 0.12.2) Suggests testthat, knitr, rmarkdown, devtools RoxygenNote 7.1.1 NeedsCompilation yes Repository CRAN Date/Publication 2021-05-26 05:10:10 UTC R topics documented: biwavelet-package . .2 ar1.spectrum . .4 ar1_ma0_sim . .5 arrow ............................................6 arrow2 . .7 1 2 biwavelet-package check.data . .7 check.datum . .8 convolve2D . .9 convolve2D_typeopen . 10 enviro.data . 10 get_minroots . 11 MOTHERS . 12 phase.plot . 12 plot.biwavelet . 13 pwtc............................................. 17 rcpp_row_quantile . 19 rcpp_wt_bases_dog . 20 rcpp_wt_bases_morlet . 21 rcpp_wt_bases_paul . 22 smooth.wavelet . 23 wclust . 24 wdist . 25 wt.............................................. 26 wt.bases . 28 wt.bases.dog . 29 wt.bases.morlet . 29 wt.bases.paul . 30 wt.sig . 31 wtc.............................................. 32 wtc.sig . 34 wtc_sig_parallel . 36 xwt ............................................. 37 Index 40 biwavelet-package Conduct Univariate and Bivariate Wavelet Analyses Description This is a port of the WTC MATLAB package written by Aslak Grinsted and the wavelet program written by Christopher Torrence and Gibert P. Compo. This package can be used to perform uni- variate and bivariate (cross-wavelet, wavelet coherence, wavelet clustering) wavelet analyses.
    [Show full text]
  • Meta4diag: Bayesian Bivariate Meta-Analysis of Diagnostic Test Studies for Routine Practice
    meta4diag: Bayesian Bivariate Meta-analysis of Diagnostic Test Studies for Routine Practice J. Guo and A. Riebler Department of Mathematical Sciences, Norwegian University of Science and Technology, Trondheim, PO 7491, Norway. July 8, 2016 Abstract This paper introduces the R package meta4diag for implementing Bayesian bivariate meta-analyses of diagnostic test studies. Our package meta4diag is a purpose-built front end of the R package INLA. While INLA offers full Bayesian inference for the large set of latent Gaussian models using integrated nested Laplace approximations, meta4diag extracts the features needed for bivariate meta-analysis and presents them in an intuitive way. It allows the user a straightforward model- specification and offers user-specific prior distributions. Further, the newly proposed penalised complexity prior framework is supported, which builds on prior intuitions about the behaviours of the variance and correlation parameters. Accurate posterior marginal distributions for sensitivity and specificity as well as all hyperparameters, and covariates are directly obtained without Markov chain Monte Carlo sampling. Further, univariate estimates of interest, such as odds ratios, as well as the SROC curve and other common graphics are directly available for interpretation. An in- teractive graphical user interface provides the user with the full functionality of the package without requiring any R programming. The package is available through CRAN https://cran.r-project.org/web/packages/meta4diag/ and its usage will be illustrated using three real data examples. arXiv:1512.06220v2 [stat.AP] 7 Jul 2016 1 1 Introduction A meta-analysis summarises the results from multiple studies with the purpose of finding a general trend across the studies.
    [Show full text]
  • Correlation and Regression Analysis
    OIC ACCREDITATION CERTIFICATION PROGRAMME FOR OFFICIAL STATISTICS Correlation and Regression Analysis TEXTBOOK ORGANISATION OF ISLAMIC COOPERATION STATISTICAL ECONOMIC AND SOCIAL RESEARCH AND TRAINING CENTRE FOR ISLAMIC COUNTRIES OIC ACCREDITATION CERTIFICATION PROGRAMME FOR OFFICIAL STATISTICS Correlation and Regression Analysis TEXTBOOK {{Dr. Mohamed Ahmed Zaid}} ORGANISATION OF ISLAMIC COOPERATION STATISTICAL ECONOMIC AND SOCIAL RESEARCH AND TRAINING CENTRE FOR ISLAMIC COUNTRIES © 2015 The Statistical, Economic and Social Research and Training Centre for Islamic Countries (SESRIC) Kudüs Cad. No: 9, Diplomatik Site, 06450 Oran, Ankara – Turkey Telephone +90 – 312 – 468 6172 Internet www.sesric.org E-mail [email protected] The material presented in this publication is copyrighted. The authors give the permission to view, copy download, and print the material presented that these materials are not going to be reused, on whatsoever condition, for commercial purposes. For permission to reproduce or reprint any part of this publication, please send a request with complete information to the Publication Department of SESRIC. All queries on rights and licenses should be addressed to the Statistics Department, SESRIC, at the aforementioned address. DISCLAIMER: Any views or opinions presented in this document are solely those of the author(s) and do not reflect the views of SESRIC. ISBN: xxx-xxx-xxxx-xx-x Cover design by Publication Department, SESRIC. For additional information, contact Statistics Department, SESRIC. i CONTENTS Acronyms
    [Show full text]
  • Multiple Linear Regression
    Multiple Linear Regression Mark Tranmer Mark Elliot 1 Contents Section 1: Introduction............................................................................................................... 3 Exam16 ...................................................................................................................................................... Exam11...................................................................................................................................... 4 Predicted values......................................................................................................................... 5 Residuals.................................................................................................................................... 6 Scatterplot of exam performance at 16 against exam performance at 11.................................. 6 1.3 Theory for multiple linear regression................................................................................... 7 Section 2: Worked Example using SPSS.................................................................................. 10 Section 3: Further topics.......................................................................................................... 36 Stepwise................................................................................................................................... 46 Section 4: BHPS assignment................................................................................................... 46 Reading list.............................................................................................................................
    [Show full text]
  • WHY MULTIVARIATE ANALYSIS? • Used for Analysing Complicated
    WHY MULTIVARIATE ANALYSIS? • Used for analysing complicated data sets • When there are many Independent Variables (IVs) and/or many Dependent Variables (DVs) • When IVs and DVs are correlated with one Dr. Azmi Mohd Tamil another to varying degrees Jabatan Kesihatan Masyarakat Fakulti Perubatan UKM • When need to come up with Prediction Model Based on lecture notes by Dr Azidah Hashim • Parallels greater complexity of contemporary research USING MULTIVARIATE ANALYSIS CHOICE OF APPROPRIATE STATISTICAL METHOD BASED ON: • WHICH STATISTICAL PROCEDURE TO USE? • Nature of IVs and DVs • HOW TO PERFORM CHOSEN PROCEDURE? • Investigator’s Experience • Personal Preferences • HOW TO INFER FROM RESULTS OBTAINED? • Ease of Comfort with Methods Used • Literature Review References • ANY OTHER ALTERNATIVE APPROACH? • Consultation with Statistician BASIC PREMISE I: BASIC PREMISE I: TERM USED:- Relationship between X Y X Y Independent Variable (IV) Dependent Variable (DV) Predictor Outcome Explanatory Response Eg . Smoking Lung Cancer Risk Factor Effect Age Hypertension Covariates (Continuous) Maternal ANC Birthweight Factor (Categorical) Control INDEPENDENT DEPENDENT VARIABLE VARIABLE Confounders RISK FACTOR OUTCOME Nuisance 1 BASIC PREMISE II: TERMINOLOGY: Univariate Analysis DATA • Analysis in which there is a single DV Bivariate Analysis Qualitative Quantitative • Analysis of two variables (categorical) • Wish to simply study the relationship between -Dichotomous -Continuous the variables - Polynomial -Discrete Multivariate Analysis • Simultaneously analyse multiple DVs and IVs Rough Guide to Multivariate BASIC PREMISE 111:VARIATIONS OF Methods (1) THE SAME THEME Name Xs Y THE GENERAL LINEAR MODEL: Regression and Continuous Continuous Correlation (eg. age) (eg. BP) Y = a + b1x1 + b2x2 + b3x3 …. + bixi Analysis of Variance Categorical Continuous Used in the following procedures: (ANOVA) (eg.
    [Show full text]
  • Geographic Information Systems Correlation Modeling As A
    Boise State University ScholarWorks Anthropology Graduate Projects and Theses Department of Anthropology 6-1-2011 Geographic Information Systems Correlation Modeling as a Management Tool in the Study Effects of Environmental Variables’ Effects on Cultural Resources Brian Wallace Boise State University GEOGRAPHIC INFORMATION SYSTEMS CORRELATION MODELING AS A MANAGEMENT TOOL IN THE STUDY EFFECTS OF ENVIRONMENTAL VARIABLES’ EFFECTS ON CULTURAL RESOURCES by Brian N. Wallace A Project Submitted in Partial Fulfillment Of the Requirements for the Degree of Master of Applied Anthropology Boise State University June 2011 BOISE STATE UNIVERSITY GRADUATE COLLEGE DEFENSE COMMITTEE APPROVAL of the project submitted by Brian N. Wallace The project presented by Brian N. Wallace entitled Geographic Information Systems Correlation Modeling as a Management Tool in the Study Effects of Environmental Variables’ Effects on Cultural Resources is hereby approved: ______________________________________________________ Mark G. Plew Date Advisor ______________________________________________________ John Ziker Date Committee Member ______________________________________________________ Margaret A. Streeter Date Committee Member ______________________________________________________ John R. (Jack) Pelton Date Dean of the Graduate College DEDICATION To my parents Jerry and Cheryl Wallace ii ACKNOWLEDGMENTS The author would like to express his appreciation to those individuals and entities who contributed to the completion of this thesis. Dr. Mark Plew is owed a great debt of gratitude for his invaluable advice, guidance, and friendship throughout my graduate experience. It is because of his helpful suggestions both for this paper and his mentorship over the last few years that have allowed me succeed beyond my own expectations. To the Cultural Resource and GIS staff of the Boise District of the Bureau of Land Management (BLM) whom allowed for access to the necessary data to run analyses successfully and for helpful guidance along the way.
    [Show full text]
  • Bivariate Analysis: the Statistical Analysis of the Relationship Between Two Variables
    bivariate analysis: The statistical analysis of the relationship between two variables. cell frequency: The number of cases in a cell of a cross-tabulation (contingency table). chi-square (χ2) test for independence: A test of statistical significance used to assess the likelihood that an observed association between two variables could have occurred by chance. consistency checking: A data-cleaning procedure involving checking for unreasonable patterns of responses, such as a 12-year-old who voted in the last US presidential election. correlation coefficient: A statistical measure of the strength and direction of a linear relationship between two variables; it may vary from −1 to 0 to +1. data cleaning: The detection and correction of errors in a computer datafile that may have occurred during data collection, coding, and/or data entry. data matrix: The form of a computer datafile, with rows as cases and columns as variables; each cell represents the value of a particular variable (column) for a particular case (row). data processing: The preparation of data for analysis. descriptive statistics: Procedures for organizing and summarizing data. dummy variable: A variable or set of variable categories recoded to have values of 0 and 1. Dummy coding may be applied to nominal- or ordinal-scale variables for the purpose of regression or other numerical analysis. frequency distribution: A tabulation of the number of cases falling into each category of a variable. histogram: A graphic display in which the height of a vertical bar represents the frequency or percentage of cases in each category of an interval/ratio variable. imputation: A procedure for handling missing data in which missing values are assigned based on other information, such as the sample mean or known values of other variables.
    [Show full text]
  • A Statistical Study to Identify the Risk Factors of Heart Attack
    tri ome cs & Bi B f i o o l s t a a n t Fatima et al., J Biom Biostat 2019, 10:2 r i s u t i o c J s Journal of Biometrics & Biostatistics ISSN: 2155-6180 Research Article Article OpenOpen Access Access A Statistical Study to Identify the Risk Factors of Heart Attack Zoha Fatima, Itrat Batool Naqvi* and Sharoon Hanook Department of Statistics, Forman Christian College, Lahore, Pakistan Abstract A statistical study has been conducted to identify the risk factors of heart attack. The study design used in this research is an observational cross sectional. A semi structured questionnaire was designed and surveyed consisting of 25 questions which were filled from 246 patients from two hospitals ‘Gulab Devi’ and ‘Jinnah Hosptial’ Lahore, Pakistan. Respondents were asked questions regarding some of the possible reasons that may cause heart attack. Out of 246 patients, 123 were cases (people who had a heart attack) and remaining 123 were control (people who only had chest pain). We took 123 patients in each group because we needed comparison. Spss and R SOFTWARE were used to determine results of this research. By using univariate, bivariate and multivariate analysis it was observed that the significant factors from model are diabetes blood pressure, sweating, heart attack before, age, severity of pain, medication and pressure of the work. Keywords: Statistical study; Heart attack; Patients; Risk factors Determination of sample size p(1− pz) 2 Introduction n = where, e2 Loss of a loved one can never be forgotten. Forgetting the pain of losing someone is almost as difficult as forgetting the person.
    [Show full text]
  • Introducing the GLIMMIX Procedure for Generalized Linear Mixed Models Oliver Schabenberger, SAS Institute Inc., Cary, NC
    NESUG 18 Analysis Introducing the GLIMMIX Procedure for Generalized Linear Mixed Models Oliver Schabenberger, SAS Institute Inc., Cary, NC ABSTRACT This paper describes a new SAS/STAT procedure for fitting models to non-normal or normal data with correlations or nonconstant variability. The GLIMMIX procedure is an add-on for the SAS/STAT product in SAS 9.1 on the Windows platform. PROC GLIMMIX extends the SAS mixed model tools in a number of ways. For example, it • models data from non-Gaussian distributions • implements low-rank smoothing based on mixed models • provides new features for LS-means comparisons and display • enables you to use SAS programming statements to compute model effects, or to define link and variance functions • fits models to multivariate data in which observations do not all have the same distribution or link Applications of the GLIMMIX procedure include estimating trends in disease rates, modeling counts or proportions over time in a clinical trial, predicting probability of occurrence in time series and spatial data, and joint modeling of correlated binary and continuous data. This paper describes generalized linear mixed models and how to use the GLIMMIX procedure for estima- tion, inference, and prediction. INTRODUCTION The GLIMMIX procedure is a new procedure in SAS/STAT software. It is an add-on for the SAS/STAT product in SAS 9.1 on the Windows platform. It is currently downloadable for the SAS 9.1 release from Software Downloads at support.sas.com. PROC GLIMMIX performs estimation and statistical inference for generalized linear mixed models (GLMMs). A generalized linear mixed model is a statistical model that extends the class of generalized linear models (GLMs) by incorporating normally distributed random effects.
    [Show full text]
  • Statistical Analysis Plan (SAP) Intention-To-Treat Analysis
    EFFECTIVENESS OF A BRIEF INFORMATION ABOUT ADVANCE DIRECTIVES IN PRIMARY CARE: A RANDOMIZED CLINICAL TRIAL Document Date: September 14, 2017 Statistical Analysis Plan (SAP) Intention-to-treat analysis. No intermediate analysis is planned. A bilateral significance level of 5%, 95% CI and 80% power will be used. The SPSS V21 program will be used. Analysis of efficacy of unadjusted parallel clinical trials. Data will be analyzed according to the CONSORT guidance and comparisons between groups will be based on the intention-to-treat principle. First: analysis of the baseline comparability of the study groups in relation to the variables studied. Descriptive statistics of all variables collected. The Student's t test or the Mann- Whitney U test (according to the frequency distribution of the variables studied) will be used to compare continuous variables, and Chi-square or Fisher's exact test for comparison of categorical variables. In the case of finding a very different distribution of the independent baseline variables (1st visit) between the two groups, which would indicate an obvious selection bias despite randomization, the possibility of analyzing as a cohort study by means of propensity score techniques . • Main objective "To assess the effect of brief oral and written information on the Advance Directives (ADs) given by primary care physicians, compared to usual clinical practice, on the proportion of persons interested in or performing VAD in 3 months in a heterogeneous group of patients (different ages and ethnicities) who go to a previous appointment ": inferential analysis of effectiveness, analysis of proportions of two samples. A multivariate logistic regression analysis with the variables of interest or realization of the yes / no VAS as the dependent variable will be performed and the group assigned as an independent variable, adjusting for the potential confounding factors.
    [Show full text]