SEM 1: Confirmatory Factor Analysis
|
|
- Susan Richard
- 5 years ago
- Views:
Transcription
1 SEM 1: Confirmatory Factor Analysis Week 3 - Measurement invariance and ordinal data Sacha Epskamp
2 General factor analysis framework: in which: y i = Λη i + ε i y N(0, Σ) η N(0, Ψ) ε N(0, Θ), y i is a p-length vector of item responses η i an m-length vector of latent variables ε i an p-length vector of residuals Λ a p m matrix of factor loadings Ψ an m m symmetric variance covariance matrix (assume always all latent variables are correlated) Θ is a p p symmetric variance covariance matrix, mostly diagonal (unless you explicitly expect violations of local independence)
3 The general framework: y i = Λη i + ε i y N(0, Σ) η N(0, Ψ) ε N(0, Θ), Allows you to derive the model-implied variance covariance matrix: Σ = ΛΨΛ + Θ
4 Two main rules: Identification Because the unit of the latent variable is unknown, we need to scale the latent variable by fixing its variance to 1 or fixing one (usually the first) factor loading to 1. We need at least as many observations (sample variances and covariances) as the number or parameters; we require non-negative degrees of freedom (DF) DF = a b a: number of observations: a = p(p + 1)/2 variances and covariances. b: number of parameters we need to estimate (do not count parameters we fixed for scaling) In general, we need 3 indicators for a single latent variable model, or 2 per factor for models with multiple (correlated) latent variables.
5 Lavaan # Install the package: install.packages("lavaan") # Load the package: library("lavaan") # Read data into R: Data <- read.csv("holzingerswineford1939.csv")
6 # cfa() model: Model <- ' visual =~ x1 + x2 + x3 textual =~ x4 + x5 + x6 speed =~ x7 + x8 + x9 '
7 # Fit in lavaan: fit <- cfa(model, Data) # Assess fit: fit ## lavaan ( ) converged normally after 35 iterations ## ## Number of observations 301 ## ## Estimator ML ## Minimum Function Test Statistic ## Degrees of freedom 24 ## P-value (Chi-square) 0.000
8 Fitting CFA models lavaan (R), Onyx and Jasp Testing for exact fit χ 2 test Assessing close fit RMSEA (below 0.5 to 0.8) SRMR (below 0.5) CFI, RNI, NFI, TLI, RFI, IFI (above 0.90 to 0.95) (A) GFI (above 0.90) Model comparison Likelihood ratio test Information criteria Modification indices
9 Lagrange Multiplier (LM) Tests mod <- modindices(fit) library("dplyr") mod %>% arrange(-mi) %>% head(10) ## lhs op rhs mi epc sepc.lv sepc.all sepc ## 1 visual =~ x ## 2 x7 ~~ x ## 3 visual =~ x ## 4 x8 ~~ x ## 5 textual =~ x ## 6 x2 ~~ x ## 7 textual =~ x ## 8 x2 ~~ x ## 9 x3 ~~ x ## 10 visual =~ x
10 # Adjusted model: Model <- ' visual =~ x1 + x2 + x3 + x9 textual =~ x4 + x5 + x6 speed =~ x7 + x8 + x9 x3 ~~ x5 '
11 # Fit in lavaan: fit <- cfa(model, Data) # Assess fit: fitmeasures(fit, c("rmsea","cfi","tli","rni","rfi","ifi","srmr","gfi")) ## rmsea cfi tli rni rfi ifi srmr gfi ##
12 Today Mean structure Multi-group CFA Measurement invariance Ordinal data
13 Including the mean structure: y i = τ + Λη i + ε i y N(µ, Σ) η N(α, Ψ) ε N(0, Θ), Allows you to derive the model-implied variance covariance matrix and means vector: Σ = ΛΨΛ + Θ µ = τ + Λα µ should resemble sample means ȳy as closely possible
14 Σ = ΛΨΛ + Θ µ = τ + Λα τ can cancel α out, hence we need to identify α = 0 Number of parameters: p(p + 1)/2 variances and covariances and p means! Number of parameters: p intercepts in τ p more observations, and p more parameters. This is why we normally ignore means!
15
16 Measurement Invariance f (y η, s) = f (y η) The distribution of y is independent on class membership s after knowing the latent trait η. Mellenbergh, G. J. (1989). Item bias and item response theory. International journal of educational research, 13(2),
17 Measurement invariance holds:
18 Measurement invariance is violated:
19 Measurement invariance is violated:
20 Multi-group model
21 Two models, for group 1: Σ 1 = Λ 1 Ψ 1 Λ 1 + Θ 1 µ 1 = τ 1 + Λ 1 α 1 With observed variance covariance matrix S 1 and observed means ȳȳȳ 1. For group 2: Σ 2 = Λ 2 Ψ 2 Λ 2 + Θ 2 µ 2 = τ 2 + Λ 2 α 2 With observed variance covariance matrix S 2 and observed means ȳȳȳ 2. Twice as many observations and parameters.
22 Steps to assess measurement invariance: Configural invariance: Is the configuration of the model the same? Weak Invariance: Are factor loadings the same? Strong Invariance: Are the intercepts the same? Strict Invariance: Are the residual variances the same?
23 Wicherts, J. M., & Dolan, C. V. (2010). Measurement invariance in confirmatory factor analysis: An illustration using IQ test performance of minorities. Educational Measurement: Issues and Practice, 29(3),
24 Wicherts, J. M., & Dolan, C. V. (2010). Measurement invariance in confirmatory factor analysis: An illustration using IQ test performance of minorities. Educational Measurement: Issues and Practice, 29(3),
25 Configural invariance: Does the same model fit in both groups?
26 Configural invariance For group 1: Σ 1 = Λ 1 Ψ 1 Λ 1 + Θ 1 µ 1 = τ 1 For group 2: Σ 2 = Λ 2 Ψ 2 Λ 2 + Θ 2 µ 2 = τ 2 (latent variable means constrained to be zero)
27 Model <- ' visual =~ x1 + x2 + x3 + x9 textual =~ x4 + x5 + x6 speed =~ x7 + x8 + x9 x3 ~~ x5 ' fit_configural <- cfa(model, Data, group = "school") fitmeasures(fit_configural, c("rmsea","cfi","tli","rni","rfi","ifi","srmr","gfi")) ## rmsea cfi tli rni rfi ifi srmr gfi ##
28 library("semplot") layout(t(1:2)) sempaths(fit_configural, "mod", "est", ask=false, layout = "tree", levels = c(1,2,4,5), edge.color = "black", reorder = FALSE, manifests = paste0("x",1:9), mar = c(2,2,4,2)) vsl txt spd vsl txt spd x1 x2 x3 x4 x5 x6 x7 x8 x9 x1 x2 x3 x4 x5 x6 x7 x8 x
29 Weak Invariance: Are factor loadings the same?
30 Weak Invariance Λ 1 = Λ 2 = Λ. For group 1: Σ 1 = ΛΨ 1 Λ + Θ 1 µ 1 = τ 1 For group 2: Σ 2 = ΛΨ 2 Λ + Θ 2 µ 2 = τ 2
31 Model <- ' visual =~ c(l1,l1)*x1 + c(l2,l2)*x2 + c(l3,l3)*x3 + c(l10,l10) * x9 textual =~ c(l4,l4)*x4 + c(l5,l5)*x5 + c(l6,l6)*x6 speed =~ c(l7,l7)*x7 + c(l8,l8)*x8 + c(l9,l9)*x9 x3 ~~ x5 ' fit_weak <- cfa(model, Data, group = "school") fit_weak ## lavaan ( ) converged normally after 41 iterations ## ## Number of observations per group ## Pasteur 156 ## Grant-White 145 ## ## Estimator ML ## Minimum Function Test Statistic ## Degrees of freedom 51 ## P-value (Chi-square) ##
32 Model <- ' visual =~ x1 + x2 + x3 + x9 textual =~ x4 + x5 + x6 speed =~ x7 + x8 + x9 x3 ~~ x5 ' fit_weak <- cfa(model, Data, group = "school", group.equal = "loadings") fit_weak ## lavaan ( ) converged normally after 41 iterations ## ## Number of observations per group ## Pasteur 156 ## Grant-White 145 ## ## Estimator ML ## Minimum Function Test Statistic ## Degrees of freedom 51 ## P-value (Chi-square) ##
33 layout(t(1:2)) sempaths(fit_weak, "mod", "est", ask=false, layout = "tree", levels = c(1,2,4,5), edge.color = "black", reorder = FALSE, manifests = paste0("x",1:9), mar = c(2,2,4,2)) vsl txt spd vsl txt spd x1 x2 x3 x4 x5 x6 x7 x8 x9 x1 x2 x3 x4 x5 x6 x7 x8 x
34 fitmeasures(fit_weak, c("rmsea","cfi","tli","rni","rfi","ifi","srmr","gfi")) ## rmsea cfi tli rni rfi ifi srmr gfi ## anova(fit_configural, fit_weak) ## Chi Square Difference Test ## ## Df AIC BIC Chisq Chisq diff Df diff Pr(>Chisq) ## fit_configural ## fit_weak (7 DF difference: 7 factor loadings)
35 Strong Invariance: Are the intercepts the same?
36 Strong Invariance τ 1 = τ 2 = τ. For group 1: Σ 1 = ΛΨ 1 Λ + Θ 1 µ 1 = τ For group 2: Σ 2 = ΛΨ 2 Λ + Θ 2 µ 2 = τ + Λα 2 α 2 is now identified!
37 Model <- ' visual =~ x1 + x2 + x3 + x9 textual =~ x4 + x5 + x6 speed =~ x7 + x8 + x9 x3 ~~ x5 ' fit_strong <- cfa(model, Data, group = "school", group.equal = c("loadings","intercepts")) fit_strong ## lavaan ( ) converged normally after 56 iterations ## ## Number of observations per group ## Pasteur 156 ## Grant-White 145 ## ## Estimator ML ## Minimum Function Test Statistic ## Degrees of freedom 57 ## P-value (Chi-square) ##
38 layout(t(1:2)) sempaths(fit_strong, "mod", "est", ask=false, layout = "tree", levels = c(1.25,2,4,5), edge.color = "black", reorder = FALSE, manifests = paste0("x",1:9), mar = c(2,2,4,2)) vsl txt spd vsl txt spd x1 x2 x3 x4 x5 x6 x7 x8 x9 x1 x2 x3 x4 x5 x6 x7 x8 x
39 fitmeasures(fit_strong, c("rmsea","cfi","tli","rni","rfi","ifi","srmr","gfi")) ## rmsea cfi tli rni rfi ifi srmr gfi ## anova(fit_configural, fit_weak, fit_strong) ## Chi Square Difference Test ## ## Df AIC BIC Chisq Chisq diff Df diff Pr(>Chisq) ## fit_configural ## fit_weak ## fit_strong e-07 ## --- ## Signif. codes: 0 '***' '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 (6 DF difference: 9 intercepts, but +3 latent variable means)
40 Where is the misfit? # Group 1: lavinspect(fit_strong, "mu")[[1]] - lavinspect(fit_strong, "sampstat")[[1]]$mean ## x1 x2 x3 x9 x4 x5 x6 x7 x8 ## # Group 2: lavinspect(fit_strong, "mu")[[2]] - lavinspect(fit_strong, "sampstat")[[2]]$mean ## x1 x2 x3 x9 x4 x5 x6 x7 x8 ##
41 Free two intercepts: Model <- ' visual =~ x1 + x2 + x3 + x9 textual =~ x4 + x5 + x6 speed =~ x7 + x8 + x9 x3 ~~ x5 x3 ~ c(t31,t32) * 1 x7 ~ c(t71,t72) * 1 ' fit_strong_mod <- cfa(model, Data, group = "school", group.equal = c("loadings","intercepts"))
42 fitmeasures(fit_strong_mod, c("rmsea","cfi","tli","rni","rfi","ifi","srmr","gfi")) ## rmsea cfi tli rni rfi ifi srmr gfi ## anova(fit_configural, fit_weak, fit_strong, fit_strong_mod) ## Chi Square Difference Test ## ## Df AIC BIC Chisq Chisq diff Df diff Pr(>Chisq) ## fit_configural ## fit_weak ## fit_strong_mod ## fit_strong e-08 ## --- ## Signif. codes: 0 '***' '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
43 Free two intercepts: Model <- ' visual =~ x1 + x2 + x3 + x9 textual =~ x4 + x5 + x6 speed =~ x7 + x8 + x9 x3 ~~ x5 ' fit_strong_mod <- cfa(model, Data, group = "school", group.equal = c("loadings","intercepts"), group.partial = c("x3 ~ 1","x7 ~ 1"))
44 fitmeasures(fit_strong_mod, c("rmsea","cfi","tli","rni","rfi","ifi","srmr","gfi")) ## rmsea cfi tli rni rfi ifi srmr gfi ## anova(fit_configural, fit_weak, fit_strong, fit_strong_mod) ## Chi Square Difference Test ## ## Df AIC BIC Chisq Chisq diff Df diff Pr(>Chisq) ## fit_configural ## fit_weak ## fit_strong_mod ## fit_strong e-08 ## --- ## Signif. codes: 0 '***' '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
45 Strict Invariance: Are the residual variances the same?
46 Strict Invariance Θ 1 = Θ 2 = Θ. For group 1: Σ 1 = ΛΨ 1 Λ + Θ µ 1 = τ For group 2: Σ 2 = ΛΨ 2 Λ + Θ µ 2 = τ + Λα 2
47 Model <- ' visual =~ x1 + x2 + x3 + x9 textual =~ x4 + x5 + x6 speed =~ x7 + x8 + x9 x3 ~~ x5 ' fit_strict <- cfa(model, Data, group = "school", group.equal = c("loadings","intercepts","residuals", "residual.covariances"), group.partial = c("x3 ~ 1","x7 ~ 1")) fit_strict ## lavaan ( ) converged normally after 61 iterations ## ## Number of observations per group ## Pasteur 156 ## Grant-White 145 ## ## Estimator ML ## Minimum Function Test Statistic ## Degrees of freedom 65
48 layout(t(1:2)) sempaths(fit_strong, "mod", "est", ask=false, layout = "tree", levels = c(1,2,4,5), edge.color = "black", reorder = FALSE, manifests = paste0("x",1:9), mar = c(2,2,4,2)) vsl txt spd vsl txt spd x1 x2 x3 x4 x5 x6 x7 x8 x9 x1 x2 x3 x4 x5 x6 x7 x8 x
49 fitmeasures(fit_strict, c("rmsea","cfi","tli","rni","rfi","ifi","srmr","gfi")) ## rmsea cfi tli rni rfi ifi srmr gfi ## anova(fit_configural, fit_weak, fit_strong_mod, fit_strict) ## Chi Square Difference Test ## ## Df AIC BIC Chisq Chisq diff Df diff Pr(>Chisq) ## fit_configural ## fit_weak ## fit_strong_mod ## fit_strict ## --- ## Signif. codes: 0 '***' '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Conclusion: Measurement variance holds adequately except for intercepts of x3 (Lozenges) and x7 (Speeded addition).
50 What about Ψ? If the goal is measurement, these do not matter (factor variances could easily differ per group) If the goal is explanatory, these (or structural relations between latents) could be constrained in a final step Group 1 Group 2 η
51 JASP DEMONSTRATION
52 In Onyx, you need to load separate datasets and copy the entire model. Very lengthy process that is much faster in Jasp or Lavaan..
53 Ordinal data If data is ordinal and consists of only a few levels of measurement data cannot be assumed normal Roughly less than five categories. Rhemtulla, M., Brosseau-Liard, P. É., & Savalei, V. (2012). When can categorical variables be treated as continuous? A comparison of robust continuous and categorical SEM estimation methods under suboptimal conditions. Psychological methods, 17(3), In this case threshold models should be used Then, it is assumed that underlying the response is a latent item that is normally distributed The covariance between this latent items and other such latent items or other continuous items can be estimated Polychoric correlation if both variables are ordinal Polyserial correlation if one item is ordinal and the other is continuous
54 I see myself as someone who is talkative
55
56
57 set.seed(1) # Setup: samplesize < cor <- 0.5 thresh1 <- c(-2,0,2) thresh2 <- c(-1,0.5,1.6) # Generate data: library("mvtnorm") cormat <- matrix(c(1,0.5,0.5,1),2,2) Data <- as.data.frame(rmvnorm(samplesize, sigma = cormat)) # Make catagorical: Data[,1] <- as.numeric(cut(data[,1],breaks = c(-inf,thresh1,inf))) Data[,2] <- as.numeric(cut(data[,2],breaks = c(-inf,thresh2,inf)))
58 # Pearson correlation: cor(data[,1], Data[,2]) ## [1] # Polychoric correlation: library("lavaan") DataOrdered <- Data DataOrdered[,1] <- ordered(data[,1]) DataOrdered[,2] <- ordered(data[,2]) lavcor(dataordered) ## V1 V2 ## V ## V
59 Polychoric correlations Lavaan will automatically treat variables that are made ordered factors via ordered() as ordinal variables and will include thresholds Alternatively, the operator can be used to define thresholds Polychoric and polyserial correlations relax the assumption of normality. However, they can sometimes go wrong! The crosstable should not have zero elements! When testing measurement invariance, now the thresholds need to be equated instead of intercepts
60 No thresholds: table(data) ## V2 ## V ## ## ## ## Zeroes.. So a bit dangerous!
61 No thresholds: Model <- ' a ~~ b ' names(data) <- c("a","b") fit <- cfa(model, Data) parameterestimates(fit) ## lhs op rhs est se z pvalue ci.lower ci.upper ## 1 a ~~ b ## 2 a ~~ a ## 3 b ~~ b
62 Thresholds: Model <- ' a ~~ b a t1 + t2 + t3 b t1 + t2 + t3 ' names(data) <- c("a","b") fit <- cfa(model, Data) parameterestimates(fit) ## lhs op rhs est se z pvalue ci.lower ci.upper ## 1 a ~~ b ## 2 a t ## 3 a t ## 4 a t ## 5 b t ## 6 b t ## 7 b t ## 8 a ~~ a NA NA ## 9 b ~~ b NA NA ## 10 a ~*~ a NA NA ## 11 b ~*~ b NA NA
63 Or use data with ordered columns: Model <- ' a ~~ b ' names(dataordered) <- c("a","b") fit <- cfa(model, DataOrdered) parameterestimates(fit) ## lhs op rhs est se z pvalue ci.lower ci.upper ## 1 a ~~ b ## 2 a t ## 3 a t ## 4 a t ## 5 b t ## 6 b t ## 7 b t ## 8 a ~~ a NA NA ## 9 b ~~ b NA NA ## 10 a ~*~ a NA NA ## 11 b ~*~ b NA NA ## 12 a ~ NA NA ## 13 b ~ NA NA
64 Sample Size How big is big enough? n : q ratio should be high Theory: to efficiently estimate lots of parameters, a larger sample is needed (5-10 per parameter) There s very little evidence that it matters (Jackson, 2003) This ratio is less important than absolute sample size n 200 people This is median SEM sample size (Shah & Goldstein, 2006) Appropriate for an average model with ML estimation Other recommendations: people minimum Use larger n if: Assumptions are violated (e.g., data are nonnormal) Model is complex (e.g., latent interactions, multilevel structure) Indicators have low reliability (factor loadings are low)
65 Big enough for what? Big enough that S is a precise estimate of Σ No estimation problems (model converges) Parameter estimates have small confidence intervals Power to detect model misspecification Chi-square test statistic has sufficient power Fit statistics are accurate
66 Power to detect non-zero parameters G-power cannot help you here: there are too many factors! To estimate power, you need to know (or estimate) the model, and all parameter values Simulation Method for Estimating Power 1. Specify a population model with all parameter values 2. Draw a large number of sample datasets of size n from this hypothetical population (e.g., 1000) 3. Fit the model to each dataset and record whether the parameter value you care about is significant 4. Count the proportion of significant parameter estimates out of 1000 datasets = power
67 Power to detect non-zero parameters G-power cannot help you here: there are too many factors! To estimate power, you need to know (or estimate) the model, and all parameter values Simulation Method for Estimating Power 1. Specify a population model with all parameter values 2. Draw a large number of sample datasets of size n from this hypothetical population (e.g., 1000) simulatedata in lavaan 3. Fit the model to each dataset and record whether the parameter value you care about is significant 4. Count the proportion of significant parameter estimates out of 1000 datasets = power
68 Power to Detect Misspecification Again, Simulation: 1. Specify a population model 2. Draw a large number of sample datasets of size n from this hypothetical population (e.g., 1000) 3. Fit a misspecified model to each dataset and record whether the chi-square test statistic is significant 4. Count the proportion of significant test statistics out of 1000 datasets = power
69 Power for Test of (Not-)Close Fit RMSEA estimates a population value Its sampling distribution has been worked out So we can put a confidence interval around it This confidence interval allows us to ask whether RMSEA is significantly different from a specified value If the population model fit is NOT CLOSE, what is power to reject H 0 by the test of close fit? If the population model fit is CLOSE, what is power to reject H 0 by the test of not-close fit? Method described in MacCallum et al. (1996) is implemented in online calculators: Power and minimum sample size for RMSEA: Power curves for RMSEA: See also findrmseasamplesize in semtools
70 Sample size required to reject RMSEA < 0.05 is the true RMSEA = 0.8 and DF = 20: library("semtools") findrmseasamplesize(rmsea0=.05, rmseaa=.08, df=20, power=0.80) ## [1] 434 Sample size required to reject RMSEA > 0.05 is the true RMSEA = 0.1 and DF = 20: findrmseasamplesize(rmsea0=.05, rmseaa=.01, df=20, power=0.80) ## [1] 474
71 Conclusion Means can be added to the CFA model In multiple-group CFA, a CFA model is fitted to several groups Measurement invariance can be assessed stepwise: Configural invariance: Is the configuration of the model the same? Weak Invariance: Are factor loadings the same? Strong Invariance: Are the intercepts the same? Strict Invariance: Are the residual variances the same? When data are ordinal, polychoric and polyserial correlations can be computed Sample size requirements are complicated, but power can be computed for RMSEA test of (non)close fit
SEM 1: Confirmatory Factor Analysis
SEM 1: Confirmatory Factor Analysis Week 3 - Measurement invariance and ordinal data Sacha Epskamp 17-04-2018 General factor analysis framework: in which: y i = Λη i + ε i y N(0, Σ) η N(0, Ψ) ε N(0, Θ),
More informationSEM 1: Confirmatory Factor Analysis
SEM 1: Confirmatory Factor Analysis Week 2 - Fitting CFA models Sacha Epskamp 10-04-2017 General factor analysis framework: in which: y i = Λη i + ε i y N(0, Σ) η N(0, Ψ) ε N(0, Θ), y i is a p-length vector
More informationThe lavaan tutorial. Yves Rosseel Department of Data Analysis Ghent University (Belgium) December 18, 2017
The lavaan tutorial Yves Rosseel Department of Data Analysis Ghent University (Belgium) December 18, 2017 Abstract If you are new to lavaan, this is the place to start. In this tutorial, we introduce the
More informationlavaan: an R package for structural equation modeling
lavaan: an R package for structural equation modeling Yves Rosseel Department of Data Analysis Belgium Modern Modeling Methods Conference May 22, 2012 University of Connecticut Yves Rosseel lavaan: an
More informationlavaan: an R package for structural equation modeling and more Version (BETA)
lavaan: an R package for structural equation modeling and more Version 0.4-9 (BETA) Yves Rosseel Department of Data Analysis Ghent University (Belgium) June 14, 2011 Abstract The lavaan package is developed
More informationUnified Visualizations of Structural Equation Models
Chapter 11 Unified Visualizations of Structural Equation Models Abstract Structural Equation Modeling (SEM) has a long history of representing models graphically as path diagrams. This chapter presents
More informationPerformance of Latent Growth Curve Models with Binary Variables
Performance of Latent Growth Curve Models with Binary Variables Jason T. Newsom & Nicholas A. Smith Department of Psychology Portland State University 1 Goal Examine estimation of latent growth curve models
More informationStudy Guide. Module 1. Key Terms
Study Guide Module 1 Key Terms general linear model dummy variable multiple regression model ANOVA model ANCOVA model confounding variable squared multiple correlation adjusted squared multiple correlation
More informationNew Features in qgraph Sacha Epskamp
New Features in qgraph Sacha Epskamp Contents New features in version.4. EBICglasso refitting............................................ Gaussian graphical model fit measures..................................
More informationMultiple Group CFA in AMOS (And Modification Indices and Nested Models)
Multiple Group CFA in AMOS (And Modification Indices and Nested Models) For this lab we will use the Self-Esteem data. An Excel file of the data is available at _www.biostat.umn.edu/~melanie/ph5482/data/index.html
More informationEstimation of Item Response Models
Estimation of Item Response Models Lecture #5 ICPSR Item Response Theory Workshop Lecture #5: 1of 39 The Big Picture of Estimation ESTIMATOR = Maximum Likelihood; Mplus Any questions? answers Lecture #5:
More informationBinary IFA-IRT Models in Mplus version 7.11
Binary IFA-IRT Models in Mplus version 7.11 Example data: 635 older adults (age 80-100) self-reporting on 7 items assessing the Instrumental Activities of Daily Living (IADL) as follows: 1. Housework (cleaning
More informationA Beginner's Guide to. Randall E. Schumacker. The University of Alabama. Richard G. Lomax. The Ohio State University. Routledge
A Beginner's Guide to Randall E. Schumacker The University of Alabama Richard G. Lomax The Ohio State University Routledge Taylor & Francis Group New York London About the Authors Preface xv xvii 1 Introduction
More informationCFA and More in R! Menglin Xu Department of Educational Studies (QREM)
CFA and More in R! Menglin Xu (xu.1384@osu.edu) Department of Educational Studies (QREM) 1 Outline Overview of lavaan and PISA data Data Screening in R, a brief overview Confirmatory factor analysis (CFA)
More informationWebSEM: Structural Equation Modeling Online
WebSEM: Structural Equation Modeling Online Zhiyong Zhang and Ke-Hai Yuan August 27, 2012 118 Haggar Hall, Department of Psychology, University of Notre Dame 1 Thanks The development of the path diagram
More informationEstimation of a hierarchical Exploratory Structural Equation Model (ESEM) using ESEMwithin-CFA
Estimation of a hierarchical Exploratory Structural Equation Model (ESEM) using ESEMwithin-CFA Alexandre J.S. Morin, Substantive Methodological Synergy Research Laboratory, Department of Psychology, Concordia
More informationHow to run a multiple indicator RI-CLPM with Mplus
How to run a multiple indicator RI-CLPM with Mplus By Ellen Hamaker September 20, 2018 The random intercept cross-lagged panel model (RI-CLPM) as proposed by Hamaker, Kuiper and Grasman (2015, Psychological
More informationCHAPTER 18 OUTPUT, SAVEDATA, AND PLOT COMMANDS
OUTPUT, SAVEDATA, And PLOT Commands CHAPTER 18 OUTPUT, SAVEDATA, AND PLOT COMMANDS THE OUTPUT COMMAND OUTPUT: In this chapter, the OUTPUT, SAVEDATA, and PLOT commands are discussed. The OUTPUT command
More informationCHAPTER 7 EXAMPLES: MIXTURE MODELING WITH CROSS- SECTIONAL DATA
Examples: Mixture Modeling With Cross-Sectional Data CHAPTER 7 EXAMPLES: MIXTURE MODELING WITH CROSS- SECTIONAL DATA Mixture modeling refers to modeling with categorical latent variables that represent
More informationMultiple Imputation with Mplus
Multiple Imputation with Mplus Tihomir Asparouhov and Bengt Muthén Version 2 September 29, 2010 1 1 Introduction Conducting multiple imputation (MI) can sometimes be quite intricate. In this note we provide
More informationThe linear mixed model: modeling hierarchical and longitudinal data
The linear mixed model: modeling hierarchical and longitudinal data Analysis of Experimental Data AED The linear mixed model: modeling hierarchical and longitudinal data 1 of 44 Contents 1 Modeling Hierarchical
More informationAnnotated multitree output
Annotated multitree output A simplified version of the two high-threshold (2HT) model, applied to two experimental conditions, is used as an example to illustrate the output provided by multitree (version
More informationSupplementary Notes on Multiple Imputation. Stephen du Toit and Gerhard Mels Scientific Software International
Supplementary Notes on Multiple Imputation. Stephen du Toit and Gerhard Mels Scientific Software International Part A: Comparison with FIML in the case of normal data. Stephen du Toit Multivariate data
More informationConducting a Path Analysis With SPSS/AMOS
Conducting a Path Analysis With SPSS/AMOS Download the PATH-INGRAM.sav data file from my SPSS data page and then bring it into SPSS. The data are those from the research that led to this publication: Ingram,
More informationModelling Proportions and Count Data
Modelling Proportions and Count Data Rick White May 4, 2016 Outline Analysis of Count Data Binary Data Analysis Categorical Data Analysis Generalized Linear Models Questions Types of Data Continuous data:
More informationBootstrapping Methods
Bootstrapping Methods example of a Monte Carlo method these are one Monte Carlo statistical method some Bayesian statistical methods are Monte Carlo we can also simulate models using Monte Carlo methods
More informationModelling Proportions and Count Data
Modelling Proportions and Count Data Rick White May 5, 2015 Outline Analysis of Count Data Binary Data Analysis Categorical Data Analysis Generalized Linear Models Questions Types of Data Continuous data:
More informationOrganizing data in R. Fitting Mixed-Effects Models Using the lme4 Package in R. R packages. Accessing documentation. The Dyestuff data set
Fitting Mixed-Effects Models Using the lme4 Package in R Deepayan Sarkar Fred Hutchinson Cancer Research Center 18 September 2008 Organizing data in R Standard rectangular data sets (columns are variables,
More informationWhy the Items versus Parcels Controversy Needn t Be One
Why the Items versus Parcels Controversy Needn t Be One Todd D. Little* University of Kansas Director, Quantitative Psychology Training Program Director, Center for Research Methods and Data Analysis Director,
More informationComparison of computational methods for high dimensional item factor analysis
Comparison of computational methods for high dimensional item factor analysis Tihomir Asparouhov and Bengt Muthén November 14, 2012 Abstract In this article we conduct a simulation study to compare several
More informationDescription Remarks and examples References Also see
Title stata.com intro 4 Substantive concepts Description Remarks and examples References Also see Description The structural equation modeling way of describing models is deceptively simple. It is deceptive
More informationConditional and Unconditional Regression with No Measurement Error
Conditional and with No Measurement Error /* reg2ways.sas */ %include 'readsenic.sas'; title2 ''; proc reg; title3 'Conditional Regression'; model infrisk = stay census; proc calis cov; /* Analyze the
More informationConfirmatory Factor Analysis on the Twin Data: Try One
Confirmatory Factor Analysis on the Twin Data: Try One /************ twinfac2.sas ********************/ TITLE2 'Confirmatory Factor Analysis'; %include 'twinread.sas'; proc calis corr; /* Analyze the correlation
More informationCHAPTER 12 ASDA ANALYSIS EXAMPLES REPLICATION-MPLUS 5.21 GENERAL NOTES ABOUT ANALYSIS EXAMPLES REPLICATION
CHAPTER 12 ASDA ANALYSIS EXAMPLES REPLICATION-MPLUS 5.21 GENERAL NOTES ABOUT ANALYSIS EXAMPLES REPLICATION These examples are intended to provide guidance on how to use the commands/procedures for analysis
More informationIntroduction to Mixed Models: Multivariate Regression
Introduction to Mixed Models: Multivariate Regression EPSY 905: Multivariate Analysis Spring 2016 Lecture #9 March 30, 2016 EPSY 905: Multivariate Regression via Path Analysis Today s Lecture Multivariate
More informationIntroduction to Factor Analysis for Marketing
Introduction to Factor Analysis for Marketing SKIM/Sawtooth Software Conference 2016, Rome Chris Chapman, Google. April 2016. Special thanks to Josh Lewandowski at Google for helpful feedback (errors are
More informationRegression Analysis and Linear Regression Models
Regression Analysis and Linear Regression Models University of Trento - FBK 2 March, 2015 (UNITN-FBK) Regression Analysis and Linear Regression Models 2 March, 2015 1 / 33 Relationship between numerical
More informationIntroduction to Mplus
Introduction to Mplus May 12, 2010 SPONSORED BY: Research Data Centre Population and Life Course Studies PLCS Interdisciplinary Development Initiative Piotr Wilk piotr.wilk@schulich.uwo.ca OVERVIEW Mplus
More informationSpatial Outlier Detection
Spatial Outlier Detection Chang-Tien Lu Department of Computer Science Northern Virginia Center Virginia Tech Joint work with Dechang Chen, Yufeng Kou, Jiang Zhao 1 Spatial Outlier A spatial data point
More informationPSY 9556B (Feb 5) Latent Growth Modeling
PSY 9556B (Feb 5) Latent Growth Modeling Fixed and random word confusion Simplest LGM knowing how to calculate dfs How many time points needed? Power, sample size Nonlinear growth quadratic Nonlinear growth
More informationANNOUNCING THE RELEASE OF LISREL VERSION BACKGROUND 2 COMBINING LISREL AND PRELIS FUNCTIONALITY 2 FIML FOR ORDINAL AND CONTINUOUS VARIABLES 3
ANNOUNCING THE RELEASE OF LISREL VERSION 9.1 2 BACKGROUND 2 COMBINING LISREL AND PRELIS FUNCTIONALITY 2 FIML FOR ORDINAL AND CONTINUOUS VARIABLES 3 THREE-LEVEL MULTILEVEL GENERALIZED LINEAR MODELS 3 FOUR
More informationSTATS PAD USER MANUAL
STATS PAD USER MANUAL For Version 2.0 Manual Version 2.0 1 Table of Contents Basic Navigation! 3 Settings! 7 Entering Data! 7 Sharing Data! 8 Managing Files! 10 Running Tests! 11 Interpreting Output! 11
More informationStatistical Analysis of List Experiments
Statistical Analysis of List Experiments Kosuke Imai Princeton University Joint work with Graeme Blair October 29, 2010 Blair and Imai (Princeton) List Experiments NJIT (Mathematics) 1 / 26 Motivation
More informationsemplot: Unified visualizations of Structural Equation Models
semplot: Unified visualizations of Structural Equation Models Sacha Epskamp University of Amsterdam Department of Psychological Methods Psychoco 204 3-02-204 semplot R package dedicated to visualizing
More informationPoisson Regression and Model Checking
Poisson Regression and Model Checking Readings GH Chapter 6-8 September 27, 2017 HIV & Risk Behaviour Study The variables couples and women_alone code the intervention: control - no counselling (both 0)
More informationRonald H. Heck 1 EDEP 606 (F2015): Multivariate Methods rev. November 16, 2015 The University of Hawai i at Mānoa
Ronald H. Heck 1 In this handout, we will address a number of issues regarding missing data. It is often the case that the weakest point of a study is the quality of the data that can be brought to bear
More informationMulticollinearity and Validation CIVL 7012/8012
Multicollinearity and Validation CIVL 7012/8012 2 In Today s Class Recap Multicollinearity Model Validation MULTICOLLINEARITY 1. Perfect Multicollinearity 2. Consequences of Perfect Multicollinearity 3.
More informationCHAPTER 1 INTRODUCTION
Introduction CHAPTER 1 INTRODUCTION Mplus is a statistical modeling program that provides researchers with a flexible tool to analyze their data. Mplus offers researchers a wide choice of models, estimators,
More informationCSSS 510: Lab 2. Introduction to Maximum Likelihood Estimation
CSSS 510: Lab 2 Introduction to Maximum Likelihood Estimation 2018-10-12 0. Agenda 1. Housekeeping: simcf, tile 2. Questions about Homework 1 or lecture 3. Simulating heteroskedastic normal data 4. Fitting
More informationHandbook of Statistical Modeling for the Social and Behavioral Sciences
Handbook of Statistical Modeling for the Social and Behavioral Sciences Edited by Gerhard Arminger Bergische Universität Wuppertal Wuppertal, Germany Clifford С. Clogg Late of Pennsylvania State University
More informationIntroduction to SAS proc calis
Introduction to SAS proc calis /* path1.sas */ %include 'SenicRead.sas'; title2 ''; /************************************************************************ * * * Cases are hospitals * * * * stay Average
More informationPath Analysis using lm and lavaan
Path Analysis using lm and lavaan Grant B. Morgan Baylor University September 10, 2014 First of all, this post is going to mirror a page on the Institute for Digital Research and Education (IDRE) site
More informationMixed Effects Models. Biljana Jonoska Stojkova Applied Statistics and Data Science Group (ASDa) Department of Statistics, UBC.
Mixed Effects Models Biljana Jonoska Stojkova Applied Statistics and Data Science Group (ASDa) Department of Statistics, UBC March 6, 2018 Resources for statistical assistance Department of Statistics
More informationChapters 5-6: Statistical Inference Methods
Chapters 5-6: Statistical Inference Methods Chapter 5: Estimation (of population parameters) Ex. Based on GSS data, we re 95% confident that the population mean of the variable LONELY (no. of days in past
More informationPackage influence.sem
Type Package Package influence.sem April 14, 2018 Title Case Influence in Structural Equation Models Version 2.2 Date 2018-04-14 Author Massimiliano Pastore & Gianmarco Altoe' Depends lavaan Suggests tcltk
More informationData Mining. ❷Chapter 2 Basic Statistics. Asso.Prof.Dr. Xiao-dong Zhu. Business School, University of Shanghai for Science & Technology
❷Chapter 2 Basic Statistics Business School, University of Shanghai for Science & Technology 2016-2017 2nd Semester, Spring2017 Contents of chapter 1 1 recording data using computers 2 3 4 5 6 some famous
More informationStandard Errors in OLS Luke Sonnet
Standard Errors in OLS Luke Sonnet Contents Variance-Covariance of ˆβ 1 Standard Estimation (Spherical Errors) 2 Robust Estimation (Heteroskedasticity Constistent Errors) 4 Cluster Robust Estimation 7
More informationChapter 6: Linear Model Selection and Regularization
Chapter 6: Linear Model Selection and Regularization As p (the number of predictors) comes close to or exceeds n (the sample size) standard linear regression is faced with problems. The variance of the
More informationPackage nonnest2. January 23, 2018
Title Tests of Non-Nested Models Version 0.5-1 Date 2018-01-22 Package nonnest2 January 23, 2018 Testing non-nested models via theory supplied by Vuong (1989) . Includes tests of model
More informationPart I. Hierarchical clustering. Hierarchical Clustering. Hierarchical clustering. Produces a set of nested clusters organized as a
Week 9 Based in part on slides from textbook, slides of Susan Holmes Part I December 2, 2012 Hierarchical Clustering 1 / 1 Produces a set of nested clusters organized as a Hierarchical hierarchical clustering
More informationPerforming Cluster Bootstrapped Regressions in R
Performing Cluster Bootstrapped Regressions in R Francis L. Huang / October 6, 2016 Supplementary material for: Using Cluster Bootstrapping to Analyze Nested Data with a Few Clusters in Educational and
More informationTitle. Description. stata.com. intro 12 Convergence problems and how to solve them
Title statacom intro 12 Convergence problems and how to solve them Description or, or, Description Remarks and examples Also see It can be devilishly difficult for software to obtain results for SEMs Here
More informationEconomics of Cybercrime The Influence of Perceived Cybercrime Risk on Online Service Adoption of European Internet Users
MÜNSTER Economics of Cybercrime The Influence of Perceived Cybercrime Risk on Online Service Adoption of European Internet Users Markus Riek Rainer Böhme Tyler Moore June 23 2014 Agenda MÜNSTER Economics
More informationFlorida State University Libraries
Florida State University Libraries Electronic Theses, Treatises and Dissertations The Graduate School 2013 Use of Item Parceling in Structural Equation Modeling with Missing Data Fatih Orcan Follow this
More informationWeek 10: Heteroskedasticity II
Week 10: Heteroskedasticity II Marcelo Coca Perraillon University of Colorado Anschutz Medical Campus Health Services Research Methods I HSMP 7607 2017 c 2017 PERRAILLON ARR 1 Outline Dealing with heteroskedasticy
More informationRepeated Measures Part 4: Blood Flow data
Repeated Measures Part 4: Blood Flow data /* bloodflow.sas */ options linesize=79 pagesize=100 noovp formdlim='_'; title 'Two within-subjecs factors: Blood flow data (NWK p. 1181)'; proc format; value
More informationThe Mplus modelling framework
The Mplus modelling framework Continuous variables Categorical variables 1 Mplus syntax structure TITLE: a title for the analysis (not part of the syntax) DATA: (required) information about the data set
More informationAdvanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras
Advanced Operations Research Prof. G. Srinivasan Department of Management Studies Indian Institute of Technology, Madras Lecture - 35 Quadratic Programming In this lecture, we continue our discussion on
More informationRegression on SAT Scores of 374 High Schools and K-means on Clustering Schools
Regression on SAT Scores of 374 High Schools and K-means on Clustering Schools Abstract In this project, we study 374 public high schools in New York City. The project seeks to use regression techniques
More informationPackage merror. November 3, 2015
Version 2.0.2 Date 2015-10-20 Package merror November 3, 2015 Author Title Accuracy and Precision of Measurements N>=3 methods are used to measure each of n items. The data are used
More informationCHAPTER 11 EXAMPLES: MISSING DATA MODELING AND BAYESIAN ANALYSIS
Examples: Missing Data Modeling And Bayesian Analysis CHAPTER 11 EXAMPLES: MISSING DATA MODELING AND BAYESIAN ANALYSIS Mplus provides estimation of models with missing data using both frequentist and Bayesian
More informationSTAT 113: Lab 9. Colin Reimer Dawson. Last revised November 10, 2015
STAT 113: Lab 9 Colin Reimer Dawson Last revised November 10, 2015 We will do some of the following together. The exercises with a (*) should be done and turned in as part of HW9. Before we start, let
More informationTHIS IS NOT REPRESNTATIVE OF CURRENT CLASS MATERIAL. STOR 455 Midterm 1 September 28, 2010
THIS IS NOT REPRESNTATIVE OF CURRENT CLASS MATERIAL STOR 455 Midterm September 8, INSTRUCTIONS: BOTH THE EXAM AND THE BUBBLE SHEET WILL BE COLLECTED. YOU MUST PRINT YOUR NAME AND SIGN THE HONOR PLEDGE
More informationUnit 5: Estimating with Confidence
Unit 5: Estimating with Confidence Section 8.3 The Practice of Statistics, 4 th edition For AP* STARNES, YATES, MOORE Unit 5 Estimating with Confidence 8.1 8.2 8.3 Confidence Intervals: The Basics Estimating
More informationGeneralized least squares (GLS) estimates of the level-2 coefficients,
Contents 1 Conceptual and Statistical Background for Two-Level Models...7 1.1 The general two-level model... 7 1.1.1 Level-1 model... 8 1.1.2 Level-2 model... 8 1.2 Parameter estimation... 9 1.3 Empirical
More informationPredictive Checking. Readings GH Chapter 6-8. February 8, 2017
Predictive Checking Readings GH Chapter 6-8 February 8, 2017 Model Choice and Model Checking 2 Questions: 1. Is my Model good enough? (no alternative models in mind) 2. Which Model is best? (comparison
More informationPanel Data 4: Fixed Effects vs Random Effects Models
Panel Data 4: Fixed Effects vs Random Effects Models Richard Williams, University of Notre Dame, http://www3.nd.edu/~rwilliam/ Last revised April 4, 2017 These notes borrow very heavily, sometimes verbatim,
More informationDS Machine Learning and Data Mining I. Alina Oprea Associate Professor, CCIS Northeastern University
DS 4400 Machine Learning and Data Mining I Alina Oprea Associate Professor, CCIS Northeastern University September 20 2018 Review Solution for multiple linear regression can be computed in closed form
More informationTHE ANALYSIS OF CONTINUOUS DATA FROM MULTIPLE GROUPS
THE ANALYSIS OF CONTINUOUS DATA FROM MULTIPLE GROUPS 1. Introduction In practice, many multivariate data sets are observations from several groups. Examples of these groups are genders, languages, political
More informationLatent Curve Models. A Structural Equation Perspective WILEY- INTERSCIENŒ KENNETH A. BOLLEN
Latent Curve Models A Structural Equation Perspective KENNETH A. BOLLEN University of North Carolina Department of Sociology Chapel Hill, North Carolina PATRICK J. CURRAN University of North Carolina Department
More informationMultivariate Analysis Multivariate Calibration part 2
Multivariate Analysis Multivariate Calibration part 2 Prof. Dr. Anselmo E de Oliveira anselmo.quimica.ufg.br anselmo.disciplinas@gmail.com Linear Latent Variables An essential concept in multivariate data
More informationMotivating Example. Missing Data Theory. An Introduction to Multiple Imputation and its Application. Background
An Introduction to Multiple Imputation and its Application Craig K. Enders University of California - Los Angeles Department of Psychology cenders@psych.ucla.edu Background Work supported by Institute
More informationSimulating power in practice
Simulating power in practice Author: Nicholas G Reich This material is part of the statsteachr project Made available under the Creative Commons Attribution-ShareAlike 3.0 Unported License: http://creativecommons.org/licenses/by-sa/3.0/deed.en
More informationFrequently Asked Questions Updated 2006 (TRIM version 3.51) PREPARING DATA & RUNNING TRIM
Frequently Asked Questions Updated 2006 (TRIM version 3.51) PREPARING DATA & RUNNING TRIM * Which directories are used for input files and output files? See menu-item "Options" and page 22 in the manual.
More information5.5 Regression Estimation
5.5 Regression Estimation Assume a SRS of n pairs (x, y ),..., (x n, y n ) is selected from a population of N pairs of (x, y) data. The goal of regression estimation is to take advantage of a linear relationship
More informationResources for statistical assistance. Quantitative covariates and regression analysis. Methods for predicting continuous outcomes.
Resources for statistical assistance Quantitative covariates and regression analysis Carolyn Taylor Applied Statistics and Data Science Group (ASDa) Department of Statistics, UBC January 24, 2017 Department
More informationFMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu
FMA901F: Machine Learning Lecture 3: Linear Models for Regression Cristian Sminchisescu Machine Learning: Frequentist vs. Bayesian In the frequentist setting, we seek a fixed parameter (vector), with value(s)
More informationEXST 7014, Lab 1: Review of R Programming Basics and Simple Linear Regression
EXST 7014, Lab 1: Review of R Programming Basics and Simple Linear Regression OBJECTIVES 1. Prepare a scatter plot of the dependent variable on the independent variable 2. Do a simple linear regression
More informationChapter 7: Dual Modeling in the Presence of Constant Variance
Chapter 7: Dual Modeling in the Presence of Constant Variance 7.A Introduction An underlying premise of regression analysis is that a given response variable changes systematically and smoothly due to
More informationSTA431 Handout 9 Double Measurement Regression on the BMI Data
STA431 Handout 9 Double Measurement Regression on the BMI Data /********************** bmi5.sas **************************/ options linesize=79 pagesize = 500 noovp formdlim='-'; title 'BMI and Health:
More informationNotes on Simulations in SAS Studio
Notes on Simulations in SAS Studio If you are not careful about simulations in SAS Studio, you can run into problems. In particular, SAS Studio has a limited amount of memory that you can use to write
More informationStat 500 lab notes c Philip M. Dixon, Week 10: Autocorrelated errors
Week 10: Autocorrelated errors This week, I have done one possible analysis and provided lots of output for you to consider. Case study: predicting body fat Body fat is an important health measure, but
More informationSection 3.2: Multiple Linear Regression II. Jared S. Murray The University of Texas at Austin McCombs School of Business
Section 3.2: Multiple Linear Regression II Jared S. Murray The University of Texas at Austin McCombs School of Business 1 Multiple Linear Regression: Inference and Understanding We can answer new questions
More informationOptimal designs for comparing curves
Optimal designs for comparing curves Holger Dette, Ruhr-Universität Bochum Maria Konstantinou, Ruhr-Universität Bochum Kirsten Schorning, Ruhr-Universität Bochum FP7 HEALTH 2013-602552 Outline 1 Motivation
More informationEstimating DCMs Using Mplus. Chapter 9 Example Data
Estimating DCMs Using Mplus 1 NCME 2012: Diagnostic Measurement Workshop Chapter 9 Example Data Example assessment 7 items Measuring 3 attributes Q matrix Item Attribute 1 Attribute 2 Attribute 3 1 1 0
More informationA noninformative Bayesian approach to small area estimation
A noninformative Bayesian approach to small area estimation Glen Meeden School of Statistics University of Minnesota Minneapolis, MN 55455 glen@stat.umn.edu September 2001 Revised May 2002 Research supported
More informationData Analysis and Solver Plugins for KSpread USER S MANUAL. Tomasz Maliszewski
Data Analysis and Solver Plugins for KSpread USER S MANUAL Tomasz Maliszewski tmaliszewski@wp.pl Table of Content CHAPTER 1: INTRODUCTION... 3 1.1. ABOUT DATA ANALYSIS PLUGIN... 3 1.3. ABOUT SOLVER PLUGIN...
More informationZ-TEST / Z-STATISTIC: used to test hypotheses about. µ when the population standard deviation is unknown
Z-TEST / Z-STATISTIC: used to test hypotheses about µ when the population standard deviation is known and population distribution is normal or sample size is large T-TEST / T-STATISTIC: used to test hypotheses
More informationPhysics 736. Experimental Methods in Nuclear-, Particle-, and Astrophysics. - Statistical Methods -
Physics 736 Experimental Methods in Nuclear-, Particle-, and Astrophysics - Statistical Methods - Karsten Heeger heeger@wisc.edu Course Schedule and Reading course website http://neutrino.physics.wisc.edu/teaching/phys736/
More informationMODELING FORCED-CHOICE DATA USING MPLUS 1
MODELING FORCED-CHOICE DATA USING MPLUS 1 Fitting a Thurstonian IRT model to forced-choice data using Mplus Anna Brown University of Cambridge Alberto Maydeu-Olivares University of Barcelona Author Note
More information