Displaying 1 result from an estimated 1 matches for "ishap".
Did you mean:
isham
2011 Aug 17
0
vglm regression with weibull distribution
...x)){
????? n<-length(x)
????? a<-p[1]
????? b<-p[2]
????? t1<-(mu-b*gamma(1+1/a))^2
????? t2<-(sig-((b^2)*gamma(1+2/a)-mu^2))^2
????? rval<-t1+t2
????? rval
? }
? optim(c(1,1),f)$par
}
With the founded initial values I have I have the message. But When try great initial value for ishape and iscale I have this
regWeib0<-vglm(eqn0,family=weibull(lshape='logoff',lscale='loge',eshape=list(offset=-2),ishape=22,iscale=6),subset(br, fdh<1))
Messages d'avis :
1: In checkwz(wz, M = M, trace = trace, wzeps = control$wzepsilon) :
? 1525 elements replaced by 1.81...