求R大神来指点下MLE

2122

收藏 2013-12-29

我在用R做一个MLE, 写的code是> test<-function(theta,da){
+ c1<-theta[1]
+ c2<-theta[2]
+ c3<-theta[3]
+ k<-nrow(da)
+ logl<- -.5*k*log(2*pi)-0.5*sum((c1*da[,1]+c2*da[,2]+c3*da[,3])**2)
+ return(-logl)
+ }
> optim(theta <- c(0,1,2), test,da, hessian=TRUE)
但是显示的错误是
错误于nrow(da) : 缺少参数"da",也没有缺省值

注：da是一个矩阵，类似于

1	0	-1
0	1	-1
1	0	1
1	-1	0
1	-1	0
-1	1	0
1	-1	0
1	-1	0
1	0	-1
1	0	-1
1	0	-1
1	0	1
0	1	-1
0	1	-1
0	-1	1
0	-1	1
0	-1	1
0	1	-1

求大神指点，究竟是哪里写错了?
如果要估计的参数从3个扩展到300个，有没有更好的写法，求给一个loop!
谢各位大神！

扫码加我拉你入群

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

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

全部回复

ntsean

2013-12-29 19:50:21

试试 optim 里面用 da = da

扫码加我拉你入群

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

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

showayvivian

2013-12-29 20:23:59

ntsean 发表于 2013-12-29 19:50
试试 optim 里面用 da = da

还是不行啊，错误于fn(par, ...) :
参数没有用(da = list(V1 = c(1, 0, 1, 1, 1, -1, 1, 1, 1, 1, 1, 1, 0, 0, 0, 0, 0, 0), V2 = c(0, 1, 0, -1, -1, 1, -1, -1, 0, 0, 0, 0, 1, 1, -1, -1, -1, 1), V3 = c(-1, -1, 1, 0, 0, 0, 0, 0, -1, -1, -1, 1, -1, -1, 1, 1, 1, -1)))

扫码加我拉你入群

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

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

qoiqpwqr

2013-12-29 21:21:24

> optim(c(0,1,2), test, da = da, hessian=TRUE)
$par
[1]  7.611328e-05 -4.540658e-05 -4.769417e-05

$value
[1] 16.54089

$counts
function gradient
   114    NA

$convergence
[1] 0

$message
NULL

$hessian
   [,1] [,2] [,3]
[1,] 11 -5 -2
[2,] -5 12 -7
[3,] -2 -7 13

扫码加我拉你入群

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

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

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

扫码加我 拉你入群

分享

扫码加好友，拉您进群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群

扫码加我拉你入群