plot 似然函数 出错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

par.est.theta <- function(sample){
n=length(sample)
log.lik<-function(theta){
  return(-n*log(pi)-sum(log(1+(theta-sample)^2)));
}
diff.f<-function(theta){
return(-2*sum((theta-sample)/(1+(theta-sample)^2)))
}
hess.f<-function(theta){
return(matrix(-2*sum((1-(theta-sample)^2)/(1+(theta-sample)^2)^2),nrow=1));
}
theta.value=seq(-10,10,by=1);
plot(theta.value,log.lik(theta.value),xlab="parameter",
        ylab="log-likelihood", type="l");
theta.est=nlminb(start=-11,-log.lik,-diff.f,-hess.f,lower=-10,upper=10);
return(theta.est);
}
a=c(1.77,-.23,2.76,3.8,3.47,56.75,-1.34,4.24,-2.44,3.29,3.711,-2.4,4.53,-.07,-1.05,13.87,-2.53,-1.75,0.27,43.21)
par.est.theta（a)

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：plot 似然函数 Theta Est ETA

上面太复杂，简化一下问题：

a=c(1.77,-.23,2.76,3.8,3.47,56.75,-1.34,4.24,-2.44,3.29,3.711,-2.4,4.53,-.07,-1.05,13.87,-2.53,-1.75,0.27,43.21)
sample=a
n=length(sample)
log.lik<-function(theta){
  return(-n*log(pi)-sum(log(1+(theta-sample)^2)));
}
theta.value=seq(-10,10,by=1);
plot(theta.value,log.lik(theta.value),xlab="parameter",
        ylab="log-likelihood", type="l");

提示：
Error in xy.coords(x, y, xlabel, ylabel, log) :
  'x' and 'y' lengths differ
In addition: Warning message:
In theta - sample :
  longer object length is not a multiple of shorter object length

就是一个简单的plot，设置x为theta.value，y为log.lik函数值。不明白为什么出错

你的log.lik的返回值只有一个数，和theta.value的长度不一样，当然会报错了。

解决办法：

1. 用循环
2. 在plot前面加上
log.lik <- Vectorize(log.lik)

qoiqpwqr 发表于 2012-2-13 06:11
解决办法：

1. 用循环

谢谢谢~
二楼 问题解决了。
可是一楼的plot还是存在问题……
Error in plot.window(...) : need finite 'xlim' values
In addition: Warning messages:
1: In min(x) : no non-missing arguments to min; returning Inf
2: In max(x) : no non-missing arguments to max; returning -Inf
3: In min(x) : no non-missing arguments to min; returning Inf
4: In max(x) : no non-missing arguments to max; returning -Inf

qoiqpwqr 发表于 2012-2-13 06:11
解决办法：

1. 用循环

没问题了，谢谢版主~~

友情帮顶

Another way of define log.lik, others remain the same:

1. log.lik <- function(theta) {return(-n*log(pi)-rowSums(1+outer(theta, sample, "-")^2))}

复制代码