scipy计算φ相关系数

1.问题描述

为了研究青年大学生性别与对某项心理测验的关系，选取了170名青年进行心理测验，得到如下表所示的数据，计算性别与测验态度的相关系数。

2.源码实现

import numpy as np
from scipy.stats import chi2

C = np.array([[22, 88], [18, 42]])

a = C[0][0]
b = C[0][1]
c = C[1][0]
d = C[1][1]

N = a + b + c + d

u = N * (a*d - b*c)**2 / ((a+b)*(a+c)*(b+d)*(c+d))
r = np.sqrt(u / N)

print(chi2.ppf(0.95, df=1))
print(u)
print(r)

3.运行及其结果

$ python3 example.py 
3.84145882069
2.15769230769
0.11266014243

4.结果解析

因为，所以求得的列联系数不具有显著意义。