def Nominator(d, c_min, c_max, c1Len, c2Len, readLen): nomin1 = (erf( (c_min + c_max + d + 1 - mean) / (2**0.5 * float(stdDev))) + erf( (mean - d - c_max - readLen - 1) / (2**0.5 * float(stdDev)))) * (pi / 2)**0.5 nomin2 = -(erf( (c_min + d + readLen + 1 - mean) / (2**0.5 * float(stdDev))) + erf( (mean - d - 2 * readLen + 1) / (2**0.5 * float(stdDev)))) * (pi / 2)**0.5 nomin = nomin1 + nomin2 return nomin
def Denominator(d,c1Len,c2Len,readLen): #term 1,2 and 3 denodes what part of the function we are integrating term1 for first (ascending), etc... term2=(c_min-readLen+1)/2.0*(erf((c_max+d+readLen-mean)/((2**0.5)*stdDev))- erf((c_min+d+readLen-mean)/((2**0.5)*stdDev)) ) first=-((pi/2)**0.5)*(d+2*readLen-mean-1)*( erf((c_min+d+readLen-mean)/(2**0.5*float(stdDev))) - erf((d+2*readLen-1-mean)/(2**0.5*float(stdDev))) ) second=stdDev*( 2.718**(-( (d+2*readLen-1-mean)**2)/(float(2*stdDev**2))) - 2.718**(-( (c_min+d+readLen-mean)**2)/(float(2*stdDev**2)))) term1=first+second first=((pi/2)**0.5)*(c_min+c_max+d-mean+1)*( erf((c_min+c_max+d-mean)/(2**0.5*float(stdDev))) - erf((c_max+readLen+d-mean)/(2**0.5*float(stdDev))) ) second=stdDev*( 2.718**(-( (c_min+c_max+d-mean)**2)/(float(2*stdDev**2))) - 2.718**(-( (c_max+readLen+d-mean)**2)/(float(2*stdDev**2)))) term3=first+second denom=term1+term2+term3 return denom
def Denom(d, readLen, mean, stdDev): first = -((pi / 2)**0.5) * (d + 2 * readLen - mean - 1) * (1 + erf( (mean - d - 2 * readLen + 1) / (2**0.5 * float(stdDev)))) second = stdDev * 2.718**(-((mean - d - 2 * readLen + 1)**2) / (float(2 * stdDev**2))) denom = first + second return denom
def Denominator(d, c1Len, c2Len, readLen): #term 1,2 and 3 denodes what part of the function we are integrating term1 for first (ascending), etc... term2 = (c_min - readLen + 1) / 2.0 * (erf( (c_max + d + readLen - mean) / ((2**0.5) * stdDev)) - erf( (c_min + d + readLen - mean) / ((2**0.5) * stdDev))) first = -((pi / 2)**0.5) * (d + 2 * readLen - mean - 1) * (erf( (c_min + d + readLen - mean) / (2**0.5 * float(stdDev))) - erf( (d + 2 * readLen - 1 - mean) / (2**0.5 * float(stdDev)))) second = stdDev * (2.718**(-((d + 2 * readLen - 1 - mean)**2) / (float(2 * stdDev**2))) - 2.718**(-((c_min + d + readLen - mean)**2) / (float(2 * stdDev**2)))) term1 = first + second first = ((pi / 2)**0.5) * (c_min + c_max + d - mean + 1) * (erf( (c_min + c_max + d - mean) / (2**0.5 * float(stdDev))) - erf( (c_max + readLen + d - mean) / (2**0.5 * float(stdDev)))) second = stdDev * (2.718**(-((c_min + c_max + d - mean)**2) / (float(2 * stdDev**2))) - 2.718**(-((c_max + readLen + d - mean)**2) / (float(2 * stdDev**2)))) term3 = first + second denom = term1 + term2 + term3 return denom
def Nominator(d,c_min,c_max,c1Len,c2Len,readLen): nomin1=( erf((c_min+c_max+d+1-mean)/(2**0.5*float(stdDev))) +erf((mean-d-c_max-readLen-1)/(2**0.5*float(stdDev))))*(pi/2)**0.5 nomin2=-(erf((c_min+d+readLen+1-mean)/(2**0.5*float(stdDev)))+erf((mean-d-2*readLen+1)/(2**0.5*float(stdDev))))*(pi/2)**0.5 nomin=nomin1+nomin2 return nomin
def Denom(d,readLen,mean,stdDev): first=-((pi/2)**0.5)*(d+2*readLen-mean-1)*(1+erf((mean-d-2*readLen+1)/(2**0.5*float(stdDev)))) second=stdDev*2.718**(-((mean-d-2*readLen+1)**2)/(float(2*stdDev**2))) denom=first+second return denom
def Nom(z,mean,stdDev): nom=-(1+erf((mean-d-2*readLen+1)/(2**0.5*float(stdDev))))*(pi/2)**0.5 return nom
def Nom(z, mean, stdDev): nom = -(1 + erf((mean - d - 2 * readLen + 1) / (2**0.5 * float(stdDev)))) * (pi / 2)**0.5 return nom