全部版块 我的主页
论坛 数据科学与人工智能 数据分析与数据科学 R语言论坛
5127 10
2015-04-07
给出一组数值,用对数极大似然函数进行参数估计,不用专门的函数,先大概给出参数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]

求大神帮忙解答


二维码

扫码加我 拉你入群

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

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

全部回复
2015-4-7 23:19:09
不错的程序!
二维码

扫码加我 拉你入群

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

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

2015-4-8 00:28:40
把函数f里面的sum去掉可以吗?
复制代码
二维码

扫码加我 拉你入群

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

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

2015-4-8 09:11:57
sum()之后就变成长度为1的数了
二维码

扫码加我 拉你入群

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

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

2015-4-8 12:19:55
xuehe 发表于 2015-4-7 23:19
不错的程序!
可是运行不了
二维码

扫码加我 拉你入群

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

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

2015-4-8 12:40:03
毓闻 发表于 2015-4-8 00:28
把函数f里面的sum去掉可以吗?
去掉了是可以运行了。可是sum可以去掉吗,极大似然函数不是每个的密度函数连乘,然后取对数的话就是加和。out<-sum(dnorm(x=data,mean=mu,sd=sigma,log=T))  里的x=data,不是x=1,前面写错了。
二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群