GAMLSS模型中的BCPE分布

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

各位大神好！
作为R语言初学者，在使用GAMLSS模型中， 出现了一些问题，还望大神解惑。
对于一个模型BCPE(X=ageλ，df（u）=9，df（σ）=4，ν=1，τ=2），为了找到某个λ值时，其GD最小。若λ=0.05，
代码为下：dp<-mydata
dp$xpower<-dp$x^0.05
m1<gamlss(y~cs(xpower,df=9),sigma.fo=~cs(xpower,df=5),nu.fo=~xpower,tau.fo=~xpower,
family=BCPE,data = dp)
此代码是否正确？还是说应该是m1<gamlss(y~cs(xpower,df=9),sigma.fo=~cs(xpower,df=5),nu.fo=~(xpower,df=3),tau.fo=~(xpower,df=2),family=BCPE,data = dp)
当ν=1，τ=2时，表示该模型为正态分布，是否应该代码应为m1<gamlss(y~cs(xpower,df=9),sigma.fo=~cs(xpower,df=5),family=NO,data = dp)
所以对于v,τ有固定值时，其自由度应该怎么算呢？
望大神指点！谢谢！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：CPE GAM R语言初学者 初学者 R语言 国际热点

楼主解决了吗？我用AIC和sbc算出来的最优自由度，如果过度拟合了，要怎么选取？
fn <- function(p) AIC(gamlss(cd4 ~ cs(age, df = p[1]), data = CD4, + trace = FALSE), k = 2) ；
opAIC <- optim(par = c(3), fn, method = "L-BFGS-B", lower = c(1), upper = c(15))
opAIC$par
这个是例子中的代码，能得到最优自由度，但是画图显示可能过度拟合