PstdII <-
function(par,sizes,phat,length){
f<-matrix(nrow=length,ncol=1)
for(i in 1:length){
z<-(par[i]-phat)/sqrt((phat*(1-phat))/sizes[i])
f[i,]<-z}
std<-f*sizes
qcc(std,type="p",sizes,center=0,limits=c(-3,3),title="P standardized phase II")
}