各位前辈：
我用R拟合了一个VGAM模型(vector generalized additive model),
fit1=vgam(cbind(nfiq1,nfiq2)~s(force),data=data[data[,3]==1,],fam=cumulative(parallel=T))
我想得到该模型在原始数据上的拟合值的标准差，调用plot(fit1,se=T,main="30um")，我得到了该模型的拟合曲线(additive部分的拟合值)并且图上还画出了拟合曲线的置信区间。但是我调用predict命令，想要得到拟合值的标准差数据，却报错了，
predict(fit1,se=T)
错误于predict.vgam(object, ...) :
cannot handle this option (se.fit = TRUE) currently
如果不能输出，那为什么plot命令却画出图来了？或者还有其他办法吗?
说明一下，VGAM包中是不输出p值的，连参数估计值的标准差输出都比较麻烦，但是为了说明效应估计的显著性我们又需要这些标准差值。希望我把问题说清楚了，请大家帮帮忙！谢谢@！！！！