用R做极大似然估计的问题~

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

      要用R语言做极大似然估计，我也知道可以用optim（）、maxLik（）等函数来做，可是如果我要给我的估计参数一个
范围的话应该怎么写？求解答~不胜感激~

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：极大似然估计 极大似然 似然估计 MAXLIK Optim 不胜感激

constrOptim

qoiqpwqr 发表于 2013-5-19 06:30
constrOptim

谢谢~可以详细点给个例子么？

niji_13 发表于 2013-5-19 15:41
谢谢~可以详细点给个例子么？

?constrOptim

的帮助文件里面有好几个例子

Here is a maximum log likelihood estimation of linear model.

set.seed(2013)
# random generate simulation data with GPD distribution with 3 parameters
alpha<- 1; beta<- 2;

n<- 100;
err<- rnorm(n, sd=0.5)
x <-  rnorm(n)


# model y defined as
y<- alpha+beta*x + err

data<- data.frame(x, y)



# define -2*loglikelihood function
LL<-function(params){
  alpha<- params[1]; beta<- params[2]; sigma<- params[3]
  
  linkfun<- alpha+beta*data$x
  
  error<-  data$y-linkfun
  
  ;
  
  f<- (1/(sigma*sqrt(2*pi))) * exp(-0.5*(error/sigma)**2)
  
  ll<- sum(log(f))
  return(-2*ll)
}

#optimized with initial condistion
optim(c(0,0,5),LL)

楼上的例子很好