###################################################
### chunk number 1: 
###################################################
#line 59 "vignettes/GeneticsDesign/inst/doc/GPC.Rnw"
library(GeneticsDesign)


###################################################
### chunk number 2: 
###################################################
#line 64 "vignettes/GeneticsDesign/inst/doc/GPC.Rnw"
set1<-seq(from=0.1, to=0.5, by=0.1)
set2<-c(1.25, 1.5, 1.75, 2.0)
len1<-length(set1)
len2<-length(set2)
mat<-matrix(0, nrow=len1, ncol=len2)
rownames(mat)<-paste("MAF=", set1, sep="")
colnames(mat)<-paste("RRAA=", set2, sep="")
for(i in 1:len1)
{ a<-set1[i]
  for(j in 1:len2)
  { b<-set2[j]
    res<-GPC.default(pA=a,pD=0.1,RRAa=(1+b)/2, RRAA=b, Dprime=1,pB=a, quiet=T)
    mat[i,j]<-res$power
  }
}
print(round(mat,3))