RDP gives P-value greater than 1

[WilliamZekaiWang]

RDP gives P-value greater than 1

[ArtPoon]

This seems to be where p-values are being calculated for the RDP method:

OpenRDP/openrdp/rdp.py

Lines 165 to 175 in bb02000

	# Calculate p_value
	val = 0
	log_n_fact = np.sum(np.log(np.arange(1, n+1))) # Convert to log space to prevent integer overflow
	for i in range(m, n):
	log_i_fact = np.sum(np.log(np.arange(1, i+1)))
	log_ni_fact = np.sum(np.log(np.arange(1, n-i+1)))
	with np.errstate(divide='ignore'): # Ignore floating point error
	val += np.math.exp((log_n_fact - (log_i_fact + log_ni_fact)) + np.log(p**n) + np.log((1-p)**(n-i)))
	uncorr_pvalue = (len_trp / n) * val
	corr_p_value = G * uncorr_pvalue

[ArtPoon]

On line 172, we should be taking the power out of $\log(p^n) = n\log(p)$ when calculating the binomial probability
Also need to look at how G is being calculated as a possible issue