请高手指点1：我自己做了数据，用广义加性模型拟合了，但不会分析！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

代码：
rm(list=ls())
data<-read.table("d:/广义加性模型试验.csv",head=TRUE,sep=",")
head(data)
colnames(data)<-c("y","x1","x2","x3","x4","x5")
data<-as.data.frame(data)
data
head(data)
library(mgcv)
library(nlme)
whole.data<-gam(y~s(x1,bs="ts")+s(x2,bs="ts")+s(x3,bs="ts")
                +s(x4,bs="ts")+s(x5,bs="ts"),
                data=data,fimily="poisson")
par(mfrow=c(3,2),mar=c(4,4,0.5,0.5))
plot(whole.data,scale=0,pages=0,select=1,
     xlab="x1",ylab="y",
     residuals=T,shade=T,lwd=2,pch=1,cex=0.5)
plot(whole.data,scale=0,pages=0,select=2,
     xlab="x2",ylab="y",
     residuals=T,shade=T,lwd=2,pch=1,cex=0.5)
plot(whole.data,scale=0,pages=0,select=3,
     xlab="x3",ylab="y",
     residuals=T,shade=T,lwd=2,pch=1,cex=0.5)
plot(whole.data,scale=0,pages=0,select=4,
     xlab="x4",ylab="y",
     residuals=T,shade=T,lwd=2,pch=1,cex=0.5)
plot(whole.data,scale=0,pages=0,select=5,
     xlab="x5",ylab="y",
     residuals=T,shade=T,lwd=2,pch=1,cex=0.5)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：广义加性模型 模型拟合 高手指点 Data list 模型 广义加性模型