Displaying 1 result from an estimated 1 matches for "reigv".
Did you mean:
reign
2003 May 08
3
Avoiding loops to spare time and memory
...%*%S%*%(1/sqrt(diag(S))*diag(items)),only.values=T)$values
for (i in 2:sample)
{
X=matrix(rnorm(cases*items),nrow=cases,byrow=F)
S=crossprod(X-rep(1,cases) %*% t(colMeans(X)))
EV=rbind(EV,eigen((1/sqrt(diag(S))*diag(items))%*%S%*%(1/sqrt(diag(S))*diag(items)),only.values=T)$values)
}
REigV=(cbind(1:items,colMeans(EV)))
REigV[,2]=as.numeric(formatC(REigV[,2],format="f",digits=7,flag="
",width=10))
colnames(REigV)=c(' ','Eigenvalue')
rownames(REigV)=rep('',items)
return(REigV)
}
# ---------------------------------------------
Thanks...