请问各位前辈们：      本人是个儿科医生，初次使用R语言和gamlss，目前想做儿童的身长生长曲线，在使用Quantile sheets建立模型qs1之后，需要使用resid(qs1)计算残差，我编的程序如下：

ppp<-findPower(lttermgirl$length,lttermgirl$PMA)
qs1<-quantSheets(length, PMA, data = lttermgirl,
                 cent=c(0.4, 2, 10, 25, 50, 75, 90, 98, 99.6),
                 x.lambda=1, p.lambda=10, logit=TRUE, power=ppp)
res1<-resid(qs1)


结果显示：
> ppp<-findPower(lttermgirl$length,lttermgirl$PMA)
*** Checking for transformation for x ***
*** power parameters  0.3489284 ***  
> qs1<-quantSheets(length, PMA, data = lttermgirl,
+                  cent=c(0.4, 2, 10, 25, 50, 75, 90, 98, 99.6),
+                  x.lambda=1, p.lambda=10, logit=TRUE, power=ppp)
% of cases below  0.4 centile is  0.07620015
% of cases below  2 centile is  1.485903
% of cases below  10 centile is  11.68402
% of cases below  25 centile is  25.86995
% of cases below  50 centile is  50.4318
% of cases below  75 centile is  74.18085
% of cases below  90 centile is  88.03658
% of cases below  98 centile is  98.0061
% of cases below  99.6 centile is  99.5682
> res1 <- resid(qs1)
Error in predMat[i, ] : 下标出界

麻烦问一下前辈们这种情况怎么处理呢？
太感谢了！！！