2014-10-28 15 views
9

ben bunu yapar Python bir test için arıyorum prop.test için:Python Oran testi R

> survivors <- matrix(c(1781,1443,135,47), ncol=2) 
> colnames(survivors) <- c('survived','died') 
> rownames(survivors) <- c('no seat belt','seat belt') 
> survivors 
      survived died 
no seat belt  1781 135 
seat belt  1443 47 
> prop.test(survivors) 

    2-sample test for equality of proportions with continuity correction 

data: survivors 
X-squared = 24.3328, df = 1, p-value = 8.105e-07 
alternative hypothesis: two.sided 
95 percent confidence interval: 
-0.05400606 -0.02382527 
sample estimates: 
    prop 1 prop 2 
0.9295407 0.9684564 

ben p-value hesaplamada çoğunlukla ilgileniyorum.

örnek Anladım düşünüyorum formu here

+0

Sorumu geliştirmek için ne yapmalıyım? – Akavall

+0

Sanırım biraz "bir bakıma benim için lütfen" sorusu gibi görünüyor. Statsmodels paketi + belgeleri de bakmak için ilginç bir yer olabilir: http://statsmodels.sourceforge.net/devel/stats.html – cel

+0

@cel, belki de bu şekilde, kendi çalışmamı yaptım ve sanırım sonunda anladım. Eğer işlevin ismini çok farklı görüyorsanız ve "orantı testi" için googling ettikten sonra bazı problemlerim vardı. – Akavall

cevap

11

alınır:

In [11]: from scipy import stats 

In [12]: import numpy as np 

In [13]: survivors = np.array([[1781,135], [1443, 47]]) 

In [14]: stats.chi2_contingency(survivors) 
Out[14]: 
(24.332761232771361,  # x-squared 
8.1048817984512269e-07, # p-value 
1, 
array([[ 1813.61832061, 102.38167939], 
     [ 1410.38167939, 79.61832061]])) 
İlgili konular