请问怎么用R画这个对数似然函数图像？？？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

logmle<-function(theta){
          n<-length(x)
          loglik<-sum(log(1-cos(x-theta)))-n*log(2*pi)
          return(loglik)
      }
x<-c(3.91,4.85,2.28,4.06,3.70,4.04,5.46,3.53,2.28,1.96,
     2.53,3.88,2.22,3.47,4.82,2.46,2.99,2.54,0.52,2.50)
curve(logmle,-pi,pi)

为什么结果显示这样：
Error in curve(logmle, -pi, pi) : 'expr'没有评估为长度是'n'的对象
此外: Warning message:
In x - theta : 长的对象长度不是短的对象长度的整倍数

R语言新手求助各位大神[em29][em29]

扫码加我 拉你入群

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

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

分享0 收藏3 回帖

关键词：对数似然 似然函数 GML Log MLE R语言 R语言新手问题 r语言入门

因为没有向量化地编写函数

1. logmle<-function(theta, par = para){
   
2.     n <- length(para)
   
3.     loglik <- sapply(theta, function(x) sum(log(1 - cos(para - x))) - n * log(2 * pi))
   
4.     return(loglik)
   
5. }
   
6. para <- c(3.91, 4.85, 2.28, 4.06, 3.70, 4.04, 5.46, 3.53, 2.28, 1.96, 2.53, 3.88, 2.22, 3.47, 4.82, 2.46, 2.99, 2.54, 0.52, 2.50)
   
7. curve(logmle, -pi, pi)

复制代码

cheetahfly 发表于 2017-6-16 20:24
因为没有向量化地编写函数

哇问题解决了！谢谢您看来还有很多要学

cheetahfly 发表于 2017-6-16 20:24
因为没有向量化地编写函数

为什么这么做

等你等你等你 发表于 2017-6-16 19:02
logmle

并不是很懂为什么里面的方程变了

再提供种方法

1. logmle<-Vectorize(function(theta){
   
2.   n<-length(x)
   
3.   loglik<-sum(log(1-cos(x-theta)))-n*log(2*pi)
   
4.   return(loglik)
   
5. })
   
6. x<-c(3.91,4.85,2.28,4.06,3.70,4.04,5.46,3.53,2.28,1.96,
   
7.      2.53,3.88,2.22,3.47,4.82,2.46,2.99,2.54,0.52,2.50)
   
8. curve(logmle,-pi,pi)

复制代码

楼主，救救孩子，我们的作业是一样的

cheetahfly 发表于 2017-6-16 20:24
因为没有向量化地编写函数

这个题是不是还有几个问题，请问怎么解决，作业求助？
2、计算出瑟塔的矩估计
3、用第二问得到的矩估计作为初始值，找到瑟塔的MLE。并用2.7和-2.7作为初始值各计算一次。
4、用类似3的对称值作为初始值，计算200次，分析这200个结果。