既然亲进来了,那应该就用过这个函数,所以就不再介绍,现在直接上问题,如果亲能帮忙解决,感激不尽!
gamlss(y~x+qrt,family=NBI,data=aids)
mod<-quote(gamlss(y ~ offset(this * x) + qrt, data = aids, family = NBI))
prof.term(mod, min=0.06, max=0.11, step=0.001)
以上是R中关于prof.term这个函数的一个例子,小弟明白此函数的作用,第一行是把NBI分布的位置参数表示成解释变量的一个线性函数,mu=a+bx+qrt...这样的。第二行就是用this指定要对x的系数 "b" 画profile deviance图,第三行是做图。
可是现在小弟的问题是我不是想针对x的系数,而是对常数项 "a" 画profile deviance图,这个this应该放在哪里,该怎么写呢?望高手指教,感激不尽! |