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
