搜索
人大经济论坛 › 附件下载

附件下载

所在主题：
文件名:  新建文本文档.txt
资料下载链接地址: https://bbs.pinggu.org/a-1340650.html
附件大小:
26.65 KB   举报本内容
实在头痛的不行了，用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)
复制代码


    熟悉论坛请点击新手指南
下载说明
1、论坛支持迅雷和网际快车等p2p多线程软件下载，请在上面选择下载通道单击右健下载即可。
2、论坛会定期自动批量更新下载地址,所以请不要浪费时间盗链论坛资源,盗链地址会很快失效。
3、本站为非盈利性质的学术交流网站,鼓励和保护原创作品，拒绝未经版权人许可的上传行为。本站如接到版权人发出的合格侵权通知，将积极的采取必要措施；同时，本站也将在技术手段和能力范围内，履行版权保护的注意义务。
(如有侵权，欢迎举报)
二维码

扫码加我 拉你入群

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

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

京ICP备16021002号-2 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问：王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2026-1-9 02:36