有没有人能救救我啊，为什么我按照federici那篇文章做到这一步全是0呀，是不是他在前边已经选择了主成分个数，所以才有时候是一维的矩阵呀mat ci=J(1,7,0)
forvalues i=1/55{
mat cih = (zi[`i',1...]-overall_average)*eigenvectors
mat ci = ci\cih
}
mat ci=ci[2...,1...]
mat list ci