if (FALSE) {
## INITIALISATION
library(xtable)
library(nFactors)
nFactors <- 3
unique <- 0.2
loadings <- 0.5
nsubjects <- 180
repsim <- 10
var <- 36
pmjc <- 12
reppar <- 10
zwick <- generateStructure(var=var, mjc=nFactors, pmjc=pmjc,
loadings=loadings,
unique=unique)
## SIMULATIONS
mzwick <- structureSim(fload=as.matrix(zwick), reppar=reppar,
repsim=repsim, details=TRUE,
N=nsubjects, quantile=0.5)
## TEST OF structureSim METHODS
is(mzwick)
summary(mzwick, index=1:5, eigenSelect=1:10, digits=3)
print(mzwick, index=1:10)
plot(x=mzwick, index=c(1:10), cex.axis=0.7, col="red")
boxplot(x=mzwick, nFactors=3, vLine="blue", col="red")
}
Run the code above in your browser using DataLab