library(OpenMx)
data(demoOneFactor)
manifests <- names(demoOneFactor)
latents <- c("G")
factorModel <- mxModel("One Factor",
type="RAM",
manifestVars = manifests,
latentVars = latents,
mxPath(from=latents, to=manifests, values=0.8),
mxPath(from=manifests, arrows=2,values=1),
mxPath(from=latents, arrows=2,
free=FALSE, values=1.0),
mxPath(from="one", to=manifests),
mxData(demoOneFactor, type="raw"))
factorModelFit <- mxRun(factorModel)
indModel <- factorModelFit
indModel$A$values['x1','G'] <- 0.3
indModel$A$free['x1','G'] <- FALSE
indModel <- mxRun(indModel)
got <- mxPower(factorModelFit, indModel, power=NA, n=100, method='ncp')
detail <- attr(got, 'detail')
omxCheckEquals(sort(names(detail)),
c("method", "n", "power", "sig.level", "statistic"))
omxCheckEquals(detail$method, 'ncp')
omxCheckEquals(detail$n, 100)
omxCheckCloseEnough(detail$power, .9, .01)
omxCheckEquals(detail$sig.level, .05)
omxCheckEquals(detail$statistic, 'LRT')
got <- mxPower(factorModelFit, indModel, power=.8, method='ncp')
detail <- attr(got, 'detail')
omxCheckEquals(sort(names(detail)),
c("method", "n", "power", "sig.level", "statistic"))
omxCheckEquals(detail$method, 'ncp')
omxCheckEquals(detail$n, 80)
omxCheckEquals(detail$power, .8)
omxCheckEquals(detail$sig.level, .05)
omxCheckEquals(detail$statistic, 'LRT')
got <- mxPower(factorModelFit, indModel, probes=30, power=.8, n=40)
detail <- attr(got, 'detail')
omxCheckEquals(sort(names(detail)),
c("method", "n", "parameter", "parameterDiff", "power", "probes", "sig.level", "statistic"))
omxCheckEquals(detail$method, 'empirical')
omxCheckEquals(detail$n, 40)
omxCheckEquals(detail$parameter, 'One Factor.A[1,6]')
omxCheckCloseEnough(detail$parameterDiff, 0.154, .03)
omxCheckEquals(detail$power, .8)
omxCheckEquals(detail$probes, 30)
omxCheckEquals(detail$sig.level, .05)
omxCheckEquals(detail$statistic, 'LRT')
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.