Python: Paired T-Test & Wilcoxon's Signed-Rank Test

[안재형]

from scipy import stats as st

import pandas as pd

# 데이터 FT.csv 읽음

FT.csv

201B

FT=pd.read_csv('~/data/FT.csv')

## 정규성 검정

st.shapiro(FT.Postwt-FT.Prewt)

# 결과

(0.9535844922065735, 0.5156101584434509)

## Paired T-test

st.ttest_1samp(FT.Postwt-FT.Prewt, popmean=0)

# 결과 겁나 간단하네요.

Ttest_1sampResult(statistic=4.184908135290033, pvalue=0.0007002531056005393)

## Wilcoxon's signed-rank test

st.wilcoxon(FT.Postwt - FT.Prewt)  또는

st.wilcoxon(FT.Postwt, FT.Prewt)

# 결과 - exact p-value

WilcoxonResult(statistic=11.0, pvalue=0.0008392333984375)

[cardiomoon]

주인장님, Mixed effect model 책은 언제 나오나요?

[안재형]

출판사에서 mixed models 책을 내면 안팔릴거라고 그냥 SAS책, R책에 넣고 있습니다. SAS 책부터 나와야 R책이 나올텐데요...

[강성찬]

열심히 사시는군요 ^^ 저는 40이 넘으니 아무것도 하기 싫어지는데요 ㅎㅎ

[안재형]

이런건 취미생활 정도죠^^