docs: def _qqplot_bar(M=1000000, alphaLevel = 0.05, distr = 'log10'):
  	"""calculate theoretical expectations for qqplot"""
  	mRange=10**(sp.arange(sp.log10(0.5),sp.log10(M-0.5)+0.1,0.1));#should be exp or 10**?
  	numPts=len(mRange);
  	betaalphaLevel=sp.zeros(numPts);#down in the plot
  	betaOneMinusalphaLevel=sp.zeros(numPts);#up in the plot
  	betaInvHalf=sp.zeros(numPts);
  	for n in range(numPts):
  	   m=mRange[n]; #numPLessThanThresh=m;
  	   betaInvHalf[n]=st.beta.ppf(0.5,m,M-m);
  	   betaalphaLevel[n]=st.beta.ppf(alphaLevel,m,M-m);
  	   betaOneMinusalphaLevel[n]=st.beta.ppf(1-alphaLevel,m,M-m);
  	betaDown=betaInvHalf-betaalphaLevel;
  	betaUp=betaOneMinusalphaLevel-betaInvHalf;
  
  	theoreticalPvals=mRange/M;
  	return betaUp, betaDown, theoreticalPvals