悬赏 30 个论坛币 已解决
刚接触R,做GLM分布为poisson时没有问题,但设为gamma时总报错,求问大神们为什么..谢谢
自己用excel随机做了一个test用的csv,R运行如下:
> tal.data<-read.csv("test.csv",header=T);tal.data
c1 c2 c3 c4 c5 c6 c7 o1
1 1 1 1 0 0 0 0 61.5
2 0 0 0 0 1 1 1 26.1
........
25 0 0 0 0 0 1 0 73.9
尝试了三种语句都报错..(上面是我输入的语句,红色是报错)
第一种:
> glm.gamma<-glm(o1~c1+c2+c3+c4+c5+c6+c7,family=gamma,data=tal.data)
错误于family() : 0个参数给'gamma',但它只需要1个
第二种:
> glm.gamma<-glm(o1~c1+c2+c3+c4+c5+c6+c7,family=gamma(link=inverse),data=tal.data)
错误于glm(o1 ~ c1 + c2 + c3 + c4 + c5 + c6 + c7, family = gamma(link = inverse), :
找不到对象'inverse'
第三种:
> glm.gamma<-glm(o1~c1+c2+c3+c4+c5+c6+c7,family=gamma(link=log),data=tal.data)
错误于gamma(link = log) : 提供的参数名'link'同'x'不匹配
求问1、这三种错误的原因,
以及2、怎么才能用R将gamma分布的GLM做出来呀..
论坛币不多,谢谢了!