在用R做极大似然估计是出错了,哪位大侠帮忙看看错误在那里?谢谢
源代码:
> LL<-function(params,data){
+ t1<-function(params,data){
+ t1<-params[1]*params[2](data)^(params[1]-1)*(1-data^params[1])^(params[2])
+ }
+ II<-sum(log(t1))
+ return(-II)
+ }
> dd.res<-nlminb(c(0.1,0.1),LL,data=bb1,lower=c(0.0001,-Inf,0.0001,Inf,-Inf,0.0001),upper=c(15,Inf,Inf,Inf,Inf))
错误于log(t1) : 数学函数中用了非数值变元
再往后就不能继续了。