# -*- coding: utf-8 -*-
"""
Created on Wed Nov 10 15:52:56 2021

@author: Sim
"""
import statsmodels.stats.proportion as prop
import scipy.stats as st
import numpy as np

n=525; x=357
p0 = 0.55; phat=x/n

z0, pval = prop.proportions_ztest(x, n, value=p0, prop_var=p0, 
                                  alternative='larger')
print('검정값: ', z0, pval)

cip1 = st.norm.interval(.95, loc=phat, scale=np.sqrt(phat*(1-phat)/n))
print(cip1)

cip2 = prop.proportion_confint(x, n, alpha=0.05, method='normal')
print(cip2)