glxfgsh 发表于 2014-8-3 16:25
没有,如果你解决了,麻烦也跟我说一声,谢谢
发现ywh19860616提供的包只能做cre,做不了固定效应的,会报错:
错误于cbind(deparse.level, ...) : args have differing numbers of rows
此外: 警告信息:
In as.matrix.csr(x, nrow = nrow, ncol = ncol, eps = eps) :
ncol*nrow indivisable by length(x)
后来看了代码发现是有一处cbind应该改为as.matrix(对比了Roger Koenker的2004年源程序),这样就可以做固定效应了,另外det函数那段,会出之前det(),matrix.csr的错误,不知道怎么改合适,我把它注释掉了,那是个判断语句,只要保证数据没问题应该就好了~