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

附件下载

所在主题：
文件名:  Jennison_Turnbull_Sequential.pdf
资料下载链接地址: https://bbs.pinggu.org/a-1244265.html
附件大小:
Group Sequenitial Methods with Appliaction to clinical trial 作者是：Jennision 书中的常数：R如何让求？
书中在page54有说利用（3.8）和（3.10），但是不知道如何写积分上下线，尝试着用R写了code，但是结果跑出来是：错误于uniroot(R, c(0.01, 10000), tol = 1e-04) : 位于极点边的f()值之正负号不相反
觉得自己积分线还是有问题，（3.8）式期望值有正有负，该如何处理？希望高手指点。
  1. library(mvtnorm)
  2. R <- function (r)
  3. {
  4. alpha=0.05
  5. beta=0.2
  6. z2=1.977
  7. z1=sqrt(2)*z2
  8. ez1=(qnorm(1-alpha)+qnorm(1-beta))*sqrt(1/2*r)
  9. ez2=(qnorm(1-alpha)+qnorm(1-beta))*sqrt(r)
  10. lower1=c(ez1-z1,ez2-z2)
  11. upper1=c(Inf,Inf)
  12. lower2=c(-Inf,-Inf)
  13. upper2=c(z1-ez1,z2-ez2)
  14. arho=sqrt(1/2)
  15. corr=matrix(c(1,arho,arho,1),2)
  16. alpha=0.05
  17. pmvnorm(lower=lower2,upper=upper2,mean=c(0,0),corr)[1]+
  18. pmvnorm(lower=lower1,upper=upper1,mean=c(0,0),corr)[1]
  19. -2+beta
  20. }
  21. rr=uniroot(R,c(0.01,10000), tol = 0.0001)$root
复制代码


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

扫码加我 拉你入群

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

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

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

GMT+8, 2026-1-1 09:02