用R语言做极大似然估计的时候遇到一些问题

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

R语言数据.xls (89.5 KB)

2013-5-26 15:30:53 上传

数据

d<-read.csv(file=file.choose(),header=T)
>  s1<-d$ln人口密度
>  s2<-d$ln面积
>  x1<-d$GDP
>  x2<-d$距离
>  x3<-d$是否中心城区
> params<-c(0,0,0,0,0,0,0,0,0)
> #定义log-likelihood函数
> LL<-function(params,d)
+ {#参数"params"是一个向量.
+ #参数"d"，是整理之后的数据.
+ f<-(-x4)*log(1+exp(-params[1]*(params[2]*s1+params[3]*s2+params[4]*x1+params[5]*x2+params[6]*x3-params[7]*x1+params[8]*x2+params[9]*x3)))
+ ll<-sum(log(f))
+ ll<-sum(log(f))
+ #log-likelihood函数
+ return(ll)
+ #函数是最小化一个函数的值，但我们但我们是要最大化log-
+ #likeilhood函数，所以需要在“ll”前加个“-”号。
+ }
> ##最大似然值的计算
> res<-optim(params,LL,method="BFGS",hessian=TRUE,control=list(fnscale=-1))
错误于optim(params, LL, method = "BFGS", hessian = TRUE, control = list(fnscale = -1)) :
  'vmmin'的初始值不能为无穷大
此外: 警告信息：
1: In log(f) : 产生了NaNs
2: In log(f) : 产生了NaNs
> ##参数估计
> b<-res$par
错误: 找不到对象'res'
> hhh<-res$hessian
错误: 找不到对象'res'
> ##t值的计算
> tval<-b/sqrt(-diag(solve(hhh)))
错误: 找不到对象'b'
> ##初步似然值
> L0<-colSums（(-x4)*log(1+exp(-params[1]*(params[2]*s1+params[3]*s2+params[4]*x1+params[5]*x2+params[6]*x3-params[7]*x1+params[8]*x2+params[9]*x3)））
> ##最终似然值
> LL<-res$value
错误: 找不到对象'res'
> ##拟合度计算
> ##结果输出
> ##p^2的值
> cat("roh = ",(L0-LL)/L0)
错误于cat("roh = ", (L0 - LL)/L0) : 找不到对象'L0'
> ##修正的p^2值
> cat("rohbar=",(L0-(LL-length(b)))/L0)
错误于cat("rohbar=", (L0 - (LL - length(b)))/L0) : 找不到对象'L0'
> print(res)
错误于print(res) : 找不到对象'res'
> print(tval)
错误于print(tval) : 找不到对象'tval'




请问那个optim函数出错的原因是什么？还有此外: 警告信息：
1: In log(f) : 产生了NaNs
2: In log(f) : 产生了NaNs
这是什么意思

下面，有我改过的代码，现在把数据奉上

扫码加我 拉你入群

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

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

分享0 收藏1 回帖

关键词：极大似然估计 似然估计 极大似然 R语言 Likelihood 语言

发程序代码最好用“代码”功能，那样比较清楚。

你的LL函数里面的x4是什么？

把数据放上来看看。

qoiqpwqr 发表于 2013-5-25 23:00
你的LL函数里面的x4是什么？

把数据放上来看看。

1. >  d<-read.csv(file=file.choose(),header=T)
   
2. >  s1<-d$ln人口密度
   
3. >  s2<-d$ln面积
   
4. >  x1<-d$GDP
   
5. >  x2<-d$距离
   
6. >  x3<-d$是否中心城区
   
7. >  x4<-d$区外总出行
   
8. >  x5<-d$区外不出行
   
9. 
   
10. >params<-c(0,0,0,0,0,0,0,0,0)
    
11. > #定义log-likelihood函数
    
12. > LL<-function(params,d)
    
13. + {#参数"params"是一个向量.
    
14. + #参数"d"，是整理之后的数据.
    
15. + f<-(x4*ln(1+exp(-params[1]*(params[2]*s1+params[3]*s2+params[4]*x1+params[5]*x2+params[6]*x3-params[7]*x1-params[8]*x2-params[9]*x3))+x5*log(1+exp(params[1]*(params[2]*s1+params[3]*s2+params[4]*x1+params[]*x2+params[6]*x3-params[7]*x1-params[8]*x2-params[9]*x3)))
    
16. +ll<-sum(-f)
    
17. + #log-likelihood函数
    
18. + return(ll)
    
19. + }
    
20. >##最大似然值的计算
    
21. >res<-optim(params,LL,method="BFGS",hessian=TRUE,control=list(fnscale=-1))
    
22. >##参数估计
    
23. >b<-res$par
    
24. >hhh<-res$hessian
    
25. >##t值的计算
    
26. >tval<-b/sqrt(-diag(solve(hhh)))
    
27. >##初步似然值
    
28. >L0<-colSums（(-x4)*log(1+exp(-params[1]*(params[2]*s1+params[3]*s2+params[4]*x1+params[5]*x2+params[6]*x3-params[7]*x1-params[8]*x2-params[9]*x3)））
    
29. >##最终似然值
    
30. >LL<-res$value
    
31. >##拟合度计算
    
32. >##结果输出
    
33. >##p^2的值
    
34. >cat("roh = ",(L0-LL)/L0)
    
35. >##修正的p^2值
    
36. >cat("rohbar=",(L0-(LL-length(b)))/L0)
    
37. >print(res)
    
38. >print(tval)

复制代码

上面是我的代码

qoiqpwqr 发表于 2013-5-25 23:00
你的LL函数里面的x4是什么？

把数据放上来看看。

数据已经发了，能麻烦帮我看看么

区外不出行是哪个变量？

qoiqpwqr 发表于 2013-5-26 20:12
区外不出行是哪个变量？

是x5 啊

pjq521 发表于 2013-5-26 20:38
是x5 啊

额，我已经弄出来了，谢谢啊

pjq521 发表于 2013-5-26 15:25
上面是我的代码

我也碰到同样的问题了，你是怎么解决的？

楼上的都是怎么解决的，我也碰到了这样的问题。