-
Notifications
You must be signed in to change notification settings - Fork 1
/
Copy patht_test.py
72 lines (64 loc) · 3.56 KB
/
t_test.py
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
import numpy as np
from scipy import stats
# Sample data for stroke patients
stroke_scores = [56.66666667, 41.66666667, 38.33333333, 81.66666667, 53.33333333, 58.33333333
, 46.66666667, 43.33333333, 71.66666667, 40., 43.33333333, 56.66666667
, 66.66666667, 55., 45., 55., 83.33333333, 66.66666667, 75.
, 61.66666667, 58.33333333, 65., 43.33333333, 56.66666667, 68.33333333
, 58.33333333, 85., 68.33333333, 60., 53.33333333, 41.66666667
, 55., 81.66666667, 58.33333333, 65., 53.33333333]
# Sample data for healthy individuals
healthy_scores = [68.33333333, 45., 45., 45., 68.33333333, 50.
, 43.33333333, 56.66666667, 65., 75., 43.33333333, 43.33333333
, 56.66666667, 61.66666667, 50., 65., 55., 60.
, 75., 41.66666667, 75., 60., 53.33333333, 60.
, 40., 55., 85., 68.33333333, 58.33333333, 83.33333333
, 58.33333333, 70., 73.33333333, 71.66666667, 81.66666667, 73.33333333
, 65., 88.33333333, 60., 95., 81.66666667, 80.
, 60., 71.66666667, 88.33333333, 88.33333333, 91.66666667, 88.33333333
, 71.66666667, 75., 91.66666667, 78.33333333, 73.33333333, 58.33333333
, 73.33333333, 80., 70., 70., 78.33333333, 86.66666667
, 68.33333333, 86.66666667, 88.33333333, 91.66666667, 78.33333333, 78.33333333
, 91.66666667, 91.66666667, 70., 71.66666667, 88.33333333, 80.
, 85., 78.33333333, 88.33333333, 76.66666667, 80., 85.
, 68.33333333, 75., 70., 56.66666667, 55., 85.
, 70., 85., 66.66666667, 91.66666667, 83.33333333, 70.
, 80., 63.33333333, 56.66666667, 78.33333333, 61.66666667, 70.
, 55., 70., 85., 66.66666667, 91.66666667, 83.33333333
, 61.66666667, 96.66666667, 83.33333333, 81.66666667]
# Perform t-test
t_value, p_value = stats.ttest_ind(stroke_scores, healthy_scores)
# Print results
print("t-value: ", t_value)
print("p-value: ", p_value)
# Set significance level
alpha = 0.05
if p_value < alpha:
print("Reject the null hypothesis. There is a significant difference between the scores of stroke patients and healthy individuals.")
else:
print("Fail to reject the null hypothesis. There is no significant difference between the scores of stroke patients and healthy individuals.")
#[56.66666667, 41.66666667, 38.33333333, 81.66666667, 53.33333333, 58.33333333
#, 46.66666667, 43.33333333, 71.66666667, 40., 43.33333333, 56.66666667
#, 66.66666667, 55., 45., 55., 83.33333333, 66.66666667, 75.
#, 61.66666667, 58.33333333, 65., 43.33333333, 56.66666667, 68.33333333
#, 58.33333333, 85., 68.33333333, 60., 53.33333333, 41.66666667
#, 55., 81.66666667, 58.33333333, 65., 53.33333333]
#
#[68.33333333, 45., 45., 45., 68.33333333, 50.
#, 43.33333333, 56.66666667, 65., 75., 43.33333333, 43.33333333
#, 56.66666667, 61.66666667, 50., 65., 55., 60.
#, 75., 41.66666667, 75., 60., 53.33333333, 60.
#, 40., 55., 85., 68.33333333, 58.33333333, 83.33333333
#, 58.33333333, 70., 73.33333333, 71.66666667, 81.66666667, 73.33333333
#, 65., 88.33333333, 60., 95., 81.66666667, 80.
#, 60., 71.66666667, 88.33333333, 88.33333333, 91.66666667, 88.33333333
#, 71.66666667, 75., 91.66666667, 78.33333333, 73.33333333, 58.33333333
#, 73.33333333, 80., 70., 70., 78.33333333, 86.66666667
#, 68.33333333, 86.66666667, 88.33333333, 91.66666667, 78.33333333, 78.33333333
#, 91.66666667, 91.66666667, 70., 71.66666667, 88.33333333, 80.
#, 85., 78.33333333, 88.33333333, 76.66666667, 80., 85.
#, 68.33333333, 75., 70., 56.66666667, 55., 85.
#, 70., 85., 66.66666667, 91.66666667, 83.33333333, 70.
#, 80., 63.33333333, 56.66666667, 78.33333333, 61.66666667, 70.
#, 55., 70., 85., 66.66666667, 91.66666667, 83.33333333
#, 61.66666667, 96.66666667, 83.33333333, 81.66666667]