关于极大似然估计

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

实在头痛的不行了，用R软件做的极大似然估计，程序没有问题，数据没有问题，为什么就是运行不出来？坑爹啊程序和数据都附上
下面是我的程序，求大神指教我到底错哪里了

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

复制代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：极大似然估计 极大似然 似然估计 Likelihood length

怎么就没人帮我弄啊