谁帮着修改下程序（R软件实现样条拟合）

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

library(foreign)
> library(nlme)
> file<-read.spss("file.sav",to.data.frame=TRUE)
> attach(file)
> x <- file$LaggedValue
> y <- file$StockIndex
> default.knots <- function(x,num.knots)
+ {
+ if (missing(num.knots))
+ num.knots <- max(5,min(floor(length(unique(x))/4),35))
+ return(quantile(unique(x),seq(0,1,length=
+ (num.knots+2))[-c(1,(num.knots+2))]))
+ }
> knots<-default.knots(x)
> n <- length(x)
> X <- cbind(rep(1,n),x)
> Z <- outer(x,knots,"-")
> Z <- Z*(Z>0)
> fit <- lme(y~-1+X,random=pdIdent(~-1+Z))
错误于getGroups.data.frame(dataMix, groups) :
  Invalid formula for groups
> beta.hat <- fit$coef$fixed
错误: 找不到对象'fit'
> u.hat <- unlist(fit$coef$random)
错误于unlist(fit$coef$random) : 找不到对象'fit'
> f.hat <- X%*%beta.hat + Z%*%u.hat
错误: 找不到对象'beta.hat'
> sig.eps.hat <- fit$sigma
错误: 找不到对象'fit'
> sig.u.hat <- sig.eps.hat*exp(unlist(fit$modelStruct))

运行到fit那一栏就不行了，lme（）这里面有问题吧，请各位高手帮帮忙

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：r软件 Library Default Foreign attach 程序 软件 拟合