mumat13 <- cbind(G1=c(30,35), G0=c(20,28))
varmat13 <- cbind(G1=c(100,100),G0=c(100,100))
dtr13 <- smartDTR(mu_Barm=mumat13, sigsq_Barm=varmat13,
Barm=c(1,3), nsubject=252, pG_A1=0.8)
mumat24 <- cbind(G1=c(25,32), G0=c(18,23))
varmat24 <- cbind(G1=c(100,100),G0=c(100,100))
dtr24 <- smartDTR(mu_Barm=mumat24, sigsq_Barm=varmat24,
Barm=c(2,4), nsubject=252, pG_A1=0.8, pG_A2=0.8)
pdtr13vs24 <- powerDTR(dtr13, dtr24)
print(pdtr13vs24) ## plot(pdtr13vs24)
Run the code above in your browser using DataLab