please review the code for randomization test

1 次查看（过去 30 天）

显示更早的评论

BYUNGSUK KIM 2019-12-13

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/496394-please-review-the-code-for-randomization-test

Randomization test-1.xlsx

在 MATLAB Online 中打开

Hi.

I created code for Randomization t-test (2-samples)

The final p value appears as 1 for all row data.

Please review the code below and the attached data for correctness.

The data I used is attached to Excel.

Thank you.

[code]

eA= y- yhatA; 
eB= y- yhatB; 
diff=eA.^2-eB.^2; 
meandiff=mean(diff); 
n=length(diff); 
niter = 199 
sum= 0; 
for k = 1: niter 
     randomsign=2*round(rand(1,n))-1; 
     signeddiff =randomsign. *diff; 
     meansigneddiff=mean(signeddiff); 
     sum=sum+... 
     (abs(meansigneddiff)> 
     =abs(meandiff)); 
end 
pvalue=(sum+1)/(niter+1)