-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathT_test_1.m
60 lines (41 loc) · 1.86 KB
/
T_test_1.m
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
function p = T_test_1(person_ans_heartrate_AllDay, person_IAAFT_MSEaverage_AllDay)
% 対標本t検定を実行
[h, p, ci, stats]=ttest(person_ans_heartrate_AllDay, person_IAAFT_MSEaverage_AllDay);
% 結果の表示
fprintf('t統計量: %.4f\n', stats.tstat);
fprintf('p値: %.4f\n', p);
% 有意水準を設定
alpha = 0.05;
% p値を用いて帰無仮説の棄却/採択を判断
if p < alpha
fprintf('帰無仮説を棄却: サンプル間に統計的に有意な違いがあります\n');
else
fprintf('帰無仮説を採択: サンプル間に統計的に有意な違いがありません\n');
end
% 有意水準を設定
alpha = 0.05;
% ループを使用してデータごとに対標本t検定を実行
for i = 1:10000 % 前は12000
% disp(['MSE: ',num2str(person_ans_heartrate_AllDay_1(i))]);
% disp(['IAAFTafMSE: ',num2str(person_IAAFT_MSEaverage_AllDay_1(i))]);
%
% sample1 = person_ans_heartrate_AllDay_1(i); % サンプル1から1つのデータを選択
% sample2 = person_IAAFT_MSEaverage_AllDay_1(i); % サンプル2から対応するデータを選択
x=person_ans_heartrate_AllDay(:,1);
y=person_IAAFT_MSEaverage_AllDay(:,1);
%disp(x)
[h, p(i), ci, stats] = ttest2(x, y);
% 対標本t検定を実行
% [h, p, ci, stats] = ttest2(sample1, sample2);
%結果の表示
%
% fprintf('データセット %d に対するt統計量: %.4f\n', i,stats.tstat);
% fprintf('データセット %d に対するt統計量: %.4f, p値: %.4f\n', i, stats.tstat, p);
%
% % p値を用いて帰無仮説の棄却/採択を判断
% if p < alpha
% fprintf('データセット %d: 帰無仮説を棄却, 統計的に有意な違いがあります\n', i);
% else
% fprintf('データセット %d: 帰無仮説を採択, 統計的に有意な違いがありません\n', i);
% end
end