set.seed(123456)
T=20;p=10;q=10;k=2;r=2
R=matrix(runif(p*k,min=-1,max=1),p,k)
C=matrix(runif(q*r,min=-1,max=1),q,r)
X=Y=E=array(0,c(T,p,q))
for(i in 1:T){
Y[i,,]=R%*%matrix(rnorm(k*r),k,r)%*%t(C)
E[i,,]=matrix(rnorm(p*q),p,q)
}
X=Y+E
fn=MKER(X,9)
fn$khat;
fn$rhat
Run the code above in your browser using DataLab