求r大神指点MLE

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我在用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!
谢各位大神！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：MLE test code COD Est function return

It seems the inputs for the function to be optimized ( here it refers to "test") should be one dimension.