Displaying 2 results from an estimated 2 matches for "kfac".
Did you mean:
fac
2007 May 03
3
factanal AIC?
...et
}
summary.multifactanal <- function(object,...){
do.call("rbind", lapply(object, summary.factanal))
}
print.multifactanal <- function(x,...){
ret <- summary.multifactanal(x)
print(ret, ...)
invisible(ret)
}
# simulate a true 4-factor model
n <- 1000
ktrue <- 4
kfac <- 5
true <- matrix(rnorm(n*ktrue), ncol=ktrue)
x <- matrix(rep(true, kfac)+rnorm(n*ktrue*kfac), ncol=ktrue*kfac)
dimnames(x) <- list(NULL, paste(rep(letters[1:ktrue], kfac), rep(1:kfac, rep(ktrue, kfac)), sep=""))
covmat <- cov.wt(x)
# run factanal for several numbers of f...
2007 May 03
3
factanal AIC?
...et
}
summary.multifactanal <- function(object,...){
do.call("rbind", lapply(object, summary.factanal))
}
print.multifactanal <- function(x,...){
ret <- summary.multifactanal(x)
print(ret, ...)
invisible(ret)
}
# simulate a true 4-factor model
n <- 1000
ktrue <- 4
kfac <- 5
true <- matrix(rnorm(n*ktrue), ncol=ktrue)
x <- matrix(rep(true, kfac)+rnorm(n*ktrue*kfac), ncol=ktrue*kfac)
dimnames(x) <- list(NULL, paste(rep(letters[1:ktrue], kfac), rep(1:kfac, rep(ktrue, kfac)), sep=""))
covmat <- cov.wt(x)
# run factanal for several numbers of f...