Description Usage Arguments Details Author(s) References Examples
Bootstrap the LRT, or any other statistic (or vector of statistics) you can extract from a fitted lavaan object.
1 2 3 4 5 6 7 8 9 10  bootstrapLavaan(object, R = 1000L, type = "ordinary", verbose = FALSE,
FUN = "coef", warn = 1L, return.boot = FALSE,
parallel = c("no", "multicore", "snow"),
ncpus = 1L, cl = NULL, iseed = NULL, h0.rmsea = NULL, ...)
bootstrapLRT(h0 = NULL, h1 = NULL, R = 1000L, type="bollen.stine",
verbose = FALSE, return.LRT = FALSE, double.bootstrap = "no",
double.bootstrap.R = 500L, double.bootstrap.alpha = 0.05,
warn = 1L, parallel = c("no", "multicore", "snow"),
ncpus = 1L, cl = NULL, iseed = NULL)

object 
An object of class 
h0 
An object of class 
h1 
An object of class 
R 
Integer. The number of bootstrap draws. 
type 
If 
FUN 
A function which when applied to the 
... 
Other named arguments for 
verbose 
If 
warn 
Sets the handling of warning messages. See 
return.boot 
Not used for now. 
return.LRT 
If 
parallel 
The type of parallel operation to be used (if any). If
missing, the default is 
ncpus 
Integer: number of processes to be used in parallel operation: typically one would chose this to the number of available CPUs. 
cl 
An optional parallel or snow cluster for use if

iseed 
An integer to set the seed. Or NULL if no reproducible seeds are
needed. To make this work, make sure the first
RNGkind() element is 
h0.rmsea 
Only used if 
double.bootstrap 
If 
double.bootstrap.R 
Integer. The number of bootstrap draws to be use for the double bootstrap. 
double.bootstrap.alpha 
The significance level to compute the adjusted alpha based on the plugin pvalues. 
The FUN function can return either a scalar or a numeric vector.
This function can be an existing function (for example coef
) or
can be a custom defined function. For example:
1 2 3 4 5 
If parallel="snow"
, it is imperative that the require(lavaan)
is included in the custom function.
Yves Rosseel and Leonard Vanbrabant. Ed Merkle contributed Yuan's bootstrap. Improvements to Yuan's bootstrap were contributed by Hao Wu and Chuchu Cheng.
Bollen, K. and Stine, R. (1992) Bootstrapping Goodness of Fit Measures in Structural Equation Models. Sociological Methods and Research, 21, 205–229.
Yuan, K.H., Hayashi, K., & Yanagihara, H. (2007). A class of population covariance matrices in the bootstrap approach to covariance structure analysis. Multivariate Behavioral Research, 42, 261–281.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18  # fit the Holzinger and Swineford (1939) example
HS.model < ' visual =~ x1 + x2 + x3
textual =~ x4 + x5 + x6
speed =~ x7 + x8 + x9 '
fit < cfa(HS.model, data=HolzingerSwineford1939, se="none")
# get the test statistic for the original sample
T.orig < fitMeasures(fit, "chisq")
# bootstrap to get bootstrap test statistics
# we only generate 10 bootstrap sample in this example; in practice
# you may wish to use a much higher number
T.boot < bootstrapLavaan(fit, R=10, type="bollen.stine",
FUN=fitMeasures, fit.measures="chisq")
# compute a bootstrap based pvalue
pvalue.boot < length(which(T.boot > T.orig))/length(T.boot)

Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.