# -*- coding: utf-8 -*-
"""
Created on Wed Jan  5 11:44:17 2022

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

obs = np.array([32, 2, 28, 14])
pp = np.array([0.342, 0.271, 0.274, 0.113])       
expd = sum(obs) * pp

chi0, pval = st.chisquare(obs, expd)
print(chi0, pval)