偏斜正态分布的极大似然估计的数值模拟程序

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

SampleNorm=function(n,mu,sigma,lambda){
    U=rnorm(n,mean=mu,sd=sigma)
    Gu=pnorm(lambda*U,mean=mu,sd=sigma)
    p=runif(n);X1=U
    for(i in 1:n){if(p[i]>=Gu[i]){X1[i]=-X1[i]}}
    hist(X1)
    return(X1)   
}
SkewNorm=function(parameter){
  n=length(bb)
  mu=parameter[1]
  sigma=parameter[2]
  lambda=parameter[3]
  p1=dnorm(bb,mu,sigma)
  e1=sum(log(p1))
  p2=pnorm(lambda*bb,mu,sigma)
  e2=sum(log(p2))
  Lfunction=n*log(2)+e1+e2
  return(Lfunction)
}
b<-c(1:300)
b1<-as.matrix(b,ncol=3)
for(i in 1:100){
   bb=SampleNorm(30,0,1,1)
   d=SkewNorm(c(0,1,1))
   b1[i,]=nlm(SkewNorm,c(0,1,1))$estimate
   c=na.omit(b1)
}
colMeans(c)

目的是，每个样本容量n，模拟100次求平均，使估计误差相对较小。但是运行时，程序第三部分就出错了，麻烦大神看看到底该怎么修改。

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：极大似然估计 正态分布 数值模拟 似然估计 极大似然 正态分布 程序

请问解决了吗？