set.seed(1)
T=50;N=50;r=3
L=matrix(rnorm(N*r,0,1),N,r);F=matrix(rnorm(T*r,0,1),T,r)
E=matrix(rnorm(T*N,0,1),T,N)
X=F%*%t(L)+E
fit=PCA(X,3,"L")
t(fit$Lhat)%*%fit$Lhat/N
fit=PCA(X,3,"F")
t(fit$Fhat)%*%fit$Fhat/T
Run the code above in your browser using DataLab