用r画3D图求解极大似然函数估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

给出一组数值，用对数极大似然函数进行参数估计，不用专门的函数，先大概给出参数mu和sigma范围，然后希望通过画出一个3D可视图进行观测

data=seq(-5,5,0.1)
mu<-seq(-5,5,0.1)
sigma<-(1,10,0.1)

f<-function(mu,sigma)
{
out<-sum(dnorm(x=1,mean=mu,sd=sigma,log=T))
}
z<-outer(mu,sigma,f)
persp(mu,sigma,z)

运行z时总是会提示Error in outer(mu, sigma, lognorm) :
  dims [product 5151] do not match the length of object [1]

求大神帮忙解答

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：极大似然 似然函数 3D图 Sigma Data function sigma Error

不错的程序！

把函数f里面的sum去掉可以吗？

1. f<-function(mu,sigma)
   
2. {
   
3. out<-dnorm(x=1,mean=mu,sd=sigma,log=T)
   
4. }

复制代码

sum()之后就变成长度为1的数了

xuehe 发表于 2015-4-7 23:19
不错的程序！

可是运行不了

毓闻 发表于 2015-4-8 00:28
把函数f里面的sum去掉可以吗？

去掉了是可以运行了。可是sum可以去掉吗，极大似然函数不是每个的密度函数连乘，然后取对数的话就是加和。out<-sum(dnorm(x=data,mean=mu,sd=sigma,log=T))  里的x=data，不是x=1，前面写错了。

nuomin 发表于 2015-4-8 09:11
sum()之后就变成长度为1的数了

out<-sum(dnorm(x=data,mean=mu,sd=sigma,log=T))  不是算出一个数吗？
错误提示里的 object [1] 指的是什么呢

因为sum这个函数的存在，每次只生成一个结果，x和y的数目不匹配，试试apply这个函数吧

fwsagdj 发表于 2015-4-8 20:18
因为sum这个函数的存在，每次只生成一个结果，x和y的数目不匹配，试试apply这个函数吧

x,y是指outer里面的前两个参数吗？

zhyalo 发表于 2015-4-8 12:40
去掉了是可以运行了。可是sum可以去掉吗，极大似然函数不是每个的密度函数连乘，然后取对数的话就是加和。 ...

明白你的意思了，试试这个

1. data<-seq(-5,5,0.1)
   
2. mu<-seq(-5,5,0.1)
   
3. sigma<-seq(1,10,0.1)
   
4. 
   
5. f<-function(mu,sigma)
   
6. {
   
7.   rowSums(sapply(data, function(x) dnorm(x, mean = mu, sd = sigma, log = TRUE)))
   
8. }
   
9. z<-outer(mu,sigma,f)
   
10. persp(mu,sigma,z)

复制代码