求最大似然估计的问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

L<-read.csv("C:/Users/CTLab/Desktop/Left.csv",head=F)

> LL<-function(theta,L){
+ mu<-theta[1]
+ sigma<-theta[2]
+ logL<-sum(log(pnorm(L[,1],mu,sigma)-pnorm(L[,2],mu,sigma)))
+ return(-logL)
+ }
> res<-optim(c(1,1),LL，x=L)
结果·········
错误于fn(par, ...) : 参数没有用(x = list(V1 = ····
这个哪里错了啊？不明白啊，求教

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：最大似然估计 似然估计 最大似然 function return sigma

> L<-read.csv("C:/Users/CTLab/Desktop/Left.csv",head=F)
> R<-read.csv("C:/Users/CTLab/Desktop/Right.csv",head=F)
>  LL<-function(theta,L,R){
+ mu<-theta[1]
+ sigma<-theta[2]
+ logL<-sum(log(pnorm(R[,1],mu,sigma)-pnorm(L[,1],mu,sigma)))
+ return(-logL)
+ }
> res<-optim(c(1,1),LL,L=L,R=R)
错误于optim(c(1, 1), LL, L = L, R = R) : 无法用初始参数来评估函数
这个又是怎么回事啊

应该用dnorm()（不是pnorm()）

res<-optim(c(1,1),LL，x=L)
把x=L删了，
在函数定义LL<-function(theta,L)里把形参L删了。

soccy 发表于 2015-4-6 17:12
应该用dnorm()（不是pnorm()）

似然函数就是正太分布函数，dnorm不是密度函数么？

nuomin 发表于 2015-4-6 17:35
res

按照你说的做了还是不对，结果是
错误于optim(c(1, 1), LL) : 无法用初始参数来评估函数

直接用maxLik包吧，很多琐碎的工作都省掉了

soccy 发表于 2015-4-6 17:53
直接用maxLik包吧，很多琐碎的工作都省掉了

我还是想知道问题在哪，而且用maxLik省事儿很多？这个也就输入似然函数，然后初值就行啦

optim()有很多选项需要你自己调整，不是你想象的那么简单

fdcch 发表于 2015-4-6 17:44
按照你说的做了还是不对，结果是
错误于optim(c(1, 1), LL) : 无法用初始参数来评估函数

换一个初始参数，一个比较好的初始值是最小二乘估计的1.5倍。