how to transform t_value to p_value ?

Question

jinwoo lee 2021-2-3

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/734997-how-to-transform-t_value-to-p_value

回答： Hiro Yoshino 2021-2-3

在 MATLAB Online 中打开

Hi, all

I would like to know how to calculate p_value through calculated t_value.

Referring to the above formula, I calculated t_value and since it is a calculation for 118 channels, the dimension of t_value is 118x 118.

Now, I need to calculate p_value based on t_value.

% corr_RH = [112 x 118 x 118] bp_RH = [ 112 x 250 x 118] ( trials, time , ch)
for rh_idx = 1: l
    corr_RH(rh_idx,:,:) = corrcoef(squeeze(bp_RH(rh_idx,:,:)));
end
for rf_idx = 1:r
    corr_RF(rf_idx,:,:) = corrcoef(squeeze(bp_RF(rf_idx,:,:)));
end
% mean corrleation coefficient
mean_corr_RH = squeeze(mean(corr_RH, 1));
mean_corr_RF = squeeze(mean(corr_RF, 1));
% find t-statistic
% S_p : pooled standard deviation case_2 
% freedom = l + r - 2
% paper equation (6) ref.
freedom = l + r - 2;
for rh_idx = 1: l
    s_temp_rh(rh_idx, :,:) = (squeeze(corr_RH(rh_idx,:,:)) - mean_corr_RH).^2;    
end
for rf_idx = 1: l
    s_temp_rf(rf_idx, :,:) = (squeeze(corr_RF(rf_idx,:,:)) - mean_corr_RH).^2;    
end
% attenion () 
s_rh = sum(s_temp_rh,1) ./ (l-1) ;
s_rf = sum(s_temp_rf,1) ./ (r-1) ;
s_p = sqrt( ((l-1).*squeeze(s_rh) + (r-1).*squeeze(s_rf)) ./ freedom );
t_value = ( mean_corr_RH - mean_corr_RF ) ./ s_p * sqrt( 1/l + 1/r ); 
% t_value --> p_value transform 
% P_thr = 0.05
P_thr = 0.05;