2019-01-01から1ヶ月間の記事一覧

対応のあるt検定

10点テスト学習した結果、 サンプル数が一緒。対応があるt検定を行う import scipy as sp b= [9,8,10,7,5,9,10,10,8,10,10,6,8,9,10,9,10,9] a =[9,9,10,7,6,10,10,9,8,10,7,8,10,10,10,10,10,10] print(len(a),len(b)) print(sp.stats.ttest_rel( b,a)) Tte…

信頼区間

このサンプル数から95%の確率で母平均が含まれるだろう区間を求めよ。 それが信頼区間。 import scipy as sp x = [45,39,42,57,28,33,40,52] se = stats.sem(x) sp.stats.t.interval(alpha=0.95,loc=(np.mean(x)),df=(len(x)-1),scale=se) (34.1066670792881…

対応の無いt検定

対応の無いt検定 import scipy as sp b = [7,8,10,5,8,7,9,5,6,9,10,6,7,8,7,9,10,10] a = [9,9,6,10,9,8,10,7,9,10,6,8,9,9,10,7,8,8,10,9] t_value = (sp.mean(a)-sp.mean(b))/sp.sqrt((sp.var(b,ddof=1)/len(b) + sp.var(a,ddof=1)/len(a))) t_value 1.46…