Allowing the code to calculate transmissivity with a provided P-value (41343b43) · Commits · AGS_public / pythermogis

src/statistics/calculate_pvalues.py

+5 −4

Original line number	Diff line number	Diff line
		@@ -2,7 +2,7 @@ import numpy as np
		from scipy import stats


		def generate_thickness_permeability_transmissivity_for_pvalue(thickness_mean, thickness_sd, ln_permeability_mean, ln_permeability_sd, Pvalue):
		def generate_thickness_permeability_transmissivity_for_pvalue(thickness_mean, thickness_sd, ln_permeability_mean, ln_permeability_sd, Pvalue, nSamples=10000):
		if Pvalue > 1.0:
		Pvalue /= 100

		@@ -12,8 +12,9 @@ def generate_thickness_permeability_transmissivity_for_pvalue(thickness_mean, th
		ln_permeability_dist = stats.norm(loc=ln_permeability_mean, scale=ln_permeability_sd)
		permeability_pvalue = np.exp(ln_permeability_dist.ppf(1 - Pvalue))

		nSamples = 10000
		# Sampling method for transmissivity
		transmissivity_samples = np.sort(np.exp(ln_permeability_dist.rvs(nSamples) + np.log(thickness_dist.rvs(nSamples))))
		transmissivity_pvalue = transmissivity_samples[int((1 - Pvalue) * nSamples)]
		transmissivity_pvalue_sampled = transmissivity_samples[int((1 - Pvalue) * nSamples)]

		return thickness_pvalue, permeability_pvalue, transmissivity_pvalue

		return thickness_pvalue, permeability_pvalue, transmissivity_pvalue_sampled