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

附件下载

所在主题：
文件名:  testt3.dta
资料下载链接地址: https://bbs.pinggu.org/a-1757347.html
附件大小:
尝试用极大似然估计做数值运算来估计三个未知参数，在设置相同初始值的情况下，使用nlminb和optim得出的估计值差异明显，而且使用不同的数据，optim的估计值不会发生变化，似乎仅与设定的初始值有关。帮助文件看起来非常晦涩，两者不知道如何取舍，向各位大神求教！先谢过了
  1. setwd("") #设定路径
  2. library(foreign)
  3. Data<-read.dta("testt3.dta")
  4. attach(Data) # Use variables independentaly
  5. LL<- function(x,data)
  6. {
  7. u<- x[1]
  8. w<- x[2]
  9. k<- x[3]

  11. inner.sum<-0
  12. outer.sum<-0
  13. extra.sum<-0
  14. for (a in 2:229) {
  15. for (b in 1:228) {
  16. if (data[b] >= data[a]) {
  17. next
  18. } else {
  19. inner.sum <- inner.sum + w * exp(w * (data[b] - data[a]))
  20. extra.sum <- extra.sum + exp(w * data[b])-exp(w * (data[b] - 3572))
  21. }

  23. }
  24. k0_inner.sum <- k * inner.sum
  25. log_inner.sum <- log(u + k0_inner.sum)
  26. outer.sum <- outer.sum + log_inner.sum
  27. k0_extra.sum <- k * extra.sum
  28. }
  29. c<- outer.sum-u*3572-k0_extra.sum
  30. return(-c)
  31. }
  32. Data.res<- nlminb(c(0.05,0.05,0.05) , LL , data = dday) #nlminb的极大似然估计
  33. Data.res$par
  34. Data.res<- optim(c(0.05,0.05,0.05) , LL , data = dday) #optim的极大似然估计
  35. Data.res$par
复制代码
附件是一组样本数据，提供给大神测试，多谢了！！！


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

扫码加我 拉你入群

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

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

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

GMT+8, 2025-12-31 19:46