# package들을 불러들이고
from scipy import stats as st
import pandas as pd
# alcohol.csv에 입력한 데이터 (첨부됨)를 읽고
alcohol=pd.read_csv('~/data/alcohol.csv')
# 정규성 검정 (Shapiro-Wilk Test)
st.shapiro(alcohol)
# 이게 결과인데 제 R 책 4.4절 W, p-value 값과 동일합니다.
(0.9234137535095215, 0.3863089680671692)
# ttest_1sampe()이라는 함수로 mu=8.1인지 분석을 해봤습니다.
st.ttest_1samp(alcohol, popmean=8.1)
# 결과가 심하게 단순하게 t-값과 p-value만 나오는데 좀 더 많이 나오게 하는 방법은 아직 파악 못했습니다.
Ttest_1sampResult
(statistic=array([0.65299814]), pvalue=array([0.53008293]))
첨부파일첨부된 파일이 1개 있습니다.
다음검색