R软件极大似然非数值参数

8232

收藏 2012-10-03

R软件初学，求大神挑错。
> LL<-function(params,data){
+ t1<-function(params,data){
+ t1<-params[1]*(1+(params[2]*data)/params[1])^(1+1/params[2])
+ }
+ II<-sum(log(t1))
+ return(-II)
+ }
> data1.res<-nlminb(c(0.0001,0.0001),LL,data=data1[,1],
+ lower=c(0.000001,0.000001),upper=c(20,20))
错误于log(t1) : 数学函数中用了非数值参数
。。。。

扫码加我拉你入群

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

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

全部回复

南冰

2012-10-3 17:11:29

t1里面有未知参数或者非数值参数吧

扫码加我拉你入群

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

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

鱼小汐

2012-10-3 22:00:52

南冰发表于 2012-10-3 17:11
t1里面有未知参数或者非数值参数吧

本身是想用极大似然求参数的不可能避免位置参数呀。。

扫码加我拉你入群

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

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

wngbaq

2012-10-18 11:49:37

试一下这样行不行

LL<-function(params,data){

t1<-params[1]*(1+(params[2]*data)/params[1])^(1+1/params[2])

II<-sum(log(t1))
return(-II)
}
data1.res<-nlminb(c(0.0001,0.0001),LL,data=data1[,1],
lower=c(0.000001,0.000001),upper=c(20,20))

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群