import scipy.stats as stats import math # Specify the number of occurences (x), the sample size (n), and the proportion claimed in the null-hypothesis (p) x = 10 n = 40 p = 0.2 # Calculate the sample proportion p_hat = x/n # Calculate and print the test statistic print((p_hat-p)/(math.sqrt((p*(1-p))/(n))))