EBSeq/R/DenNHist.R

   1 DenNHist <-
   2 function(QList,Alpha,Beta,name,AList="F",GroupName)
   3 {
   4     if(!is.list(QList)) QList=list(QList)
   5         for (i in 1:length(QList)){
   6                 if (AList=="F") alpha.use=Alpha
   7                         if(AList=="T")  alpha.use=Alpha[i]
   8         hist(QList[[i]][QList[[i]]<.98&QList[[i]]>0],prob=T,col="blue",breaks=100,main=paste(GroupName[i],name,sep=" "),xlim=c(0,1),xlab=paste("Q alpha=",round(alpha.use,2)," beta=",round(Beta[i],2),sep=""))
   9         tmpSize=length(QList[[i]][QList[[i]]<.98])
  10         tmpseq=seq(0.001,1,length=1000)
  11         #tmpdensity=dbeta(tmpseq,AlphaResult,BetaResult[i])
  12         #points(tmpseq,tmpdensity, type="l",col="green")
  13         #ll=dbeta(tmpseq,Alpha,Beta[i])
  14         ll=tmpseq
  15                  lines(ll,dbeta(ll,alpha.use,Beta[i]),col="green",lwd=2)
  16         legend("topright",c("Data","Fitted density"),col=c("blue","green"),lwd=2)
  17 }
  18
  19         }
  20