Group Sequenitial Methods with Appliaction to clinical trial

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

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

复制代码

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：Clinical Methods Method Action clinic clinical

Jennison_Turnbull_Sequential.pdf (3.6 MB)

2013-1-1 22:52:35 上传

Jennision

谢谢