参数的联合置信域作图问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

x<-c(0.5013  1.1522  1.3654  2.2170  2.2448  2.4681  2.5502  2.7717  4.0041  4.1201  6.5891  7.3616  8.5736 10.1656 17.3799 23.0542 24.5661 29.9821 39.0952)
怎么写代码画出那个图啊，用matlab也行，就是不会写代码

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：作图问题 联合

1. x = c(0.5013, 1.1522, 1.3654, 2.2170, 2.2448, 2.4681, 2.5502,
   
2.       2.7717, 4.0041, 4.1201, 6.5891, 7.3616, 8.5736, 10.1656,
   
3.       17.3799, 23.0542, 24.5661, 29.9821, 39.0952)
   
4. ngrids = 200L
   
5. beta = seq(0.6860, 2.0791, length = ngrids)
   
6. alpha1 = 21.4064/2/colSums(outer(x,-beta,FUN="^"))
   
7. alpha2 = 54.3055/21.4064*alpha1
   
8. plot(beta, alpha1, xlim=c(0.5,2.5),ylim=c(0.5,5),type="l",
   
9.      xlab=expression(beta), ylab=expression(alpha))
   
10. lines(beta, alpha2)
    
11. segments(beta[1],alpha1[1],beta[1],alpha2[1])
    
12. segments(beta[ngrids],alpha1[ngrids],beta[ngrids],alpha2[ngrids])

复制代码