非参数回归lme调用

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

求助。。。

library(nlme)
rate.risk.df<-read.csv("D:/D/毕设/rate.risk.2.csv",,header = TRUE)
n=10344
No=rep(1,n)
# 生成节点向量 kk

K<-35
step0<-1/(K+2)
kk<-numeric(K)
kk[1]<-2*step0
for(i in 2:(K)){
  kk<- kk[i-1]+step0
}

# 生成矩阵Z
Z<-matrix( nrow = n, ncol = K)
for(j in 1:K) Z[,j]<-time(i)            # 每列 都为 time
for(i in 1:n) Z[i,]<-Z[i,]-kk     # 每行 都 减去 kk
Z[Z<0]<-0 #将所有小于0的元素变为0
# 生成数据框
PenalizedSp<-data.frame(No=No,time=rate.risk.df$time,rate=rate.risk.df$rate,Z=Z)


# 以下为随机效应分析的数据结构
PenalizedSp.group<- groupedData(rate ~ time|No, data = PenalizedSp)
#线性样条拟合
fm1<-lme(fixed= rate~time, data =PenalizedSp.group,random= pdIdent(~-1+Z))

错误于logLik.lmeStructInt(lmeSt, lmePars) : 外接函数调用时不能有NA/NaN/Inf(arg1)

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：非参数回归 LME 非参数 Library BRARY library

请问楼主解决了吗

我的也是，运行以后的参数结果中提示我有缺失值，很无奈，找不到原因