EBSeq/R/QQP.R

   1 QQP <-
   2 function(QList,AlphaResult,BetaResult,name,AList="F",GroupName){
   3
   4                     for (i in 1:length(BetaResult)){
   5                                 tmpSize=length(QList[[i]][QList[[i]]<1 & !is.na(QList[[i]])])
   6                         if (AList=="F") rdpts=rbeta(tmpSize,AlphaResult,BetaResult[i])
   7                                 else rdpts=rbeta(tmpSize,AlphaResult[i],BetaResult[i])
   8         qqplot(QList[[i]][QList[[i]]<1], rdpts,xlab="estimated q's", ylab="simulated q's from fitted beta",main=paste(name,GroupName[i],sep=" "),xlim=c(0,1),ylim=c(0,1))
   9         fit=lm(sort(rdpts)~sort(QList[[i]][QList[[i]]<1  & !is.na(QList[[i]])]))
  10         abline(fit,col="red")
  11
  12                         }
  13 }
  14