2 function(QList,Alpha,Beta,name,AList="F",GroupName)
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])
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)