sst1 <- smartSST(mu_Barm=c(G1=30, G0=20), sigsq_Barm=c(G1=16,G0=16),
Barm=1, sens=seq(.6, 1, by=.1), spec=seq(.6, 1, by=.1),
nsubject=252)
sst2 <- smartSST(mu_Barm=c(G1=20, G0=30), sigsq_Barm=c(G1=16,G0=16),
Barm=2, sens=seq(.6, 1, by=.1), spec=seq(.6, 1, by=.1),
nsubject=252)
psst12 <- powerSST(sst1, sst2)
print(psst12) ## plot(psst12)
Run the code above in your browser using DataLab