# -*- coding: utf-8 -*-
"""
Created on Sun Nov 14 11:56:40 2021

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

obs = np.array([26, 49, 25])
pp = np.array([0.25, 0.5, .25])
exp = sum(obs) * pp

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