Prev: Problem using Kalman Filter estimation
Next: please change my E-mail Id
From: Leopold on 6 Aug 2010 09:19
Hallo
I am new to Matlab and this forum:
i'd like to convert the following code of a two-way Scheirer-Ray-Hare-Test (non-parametric test similiar to ANOVA) into the three-way version:

[function [p,H,p2] = schreirer_hare(data,grp)

[R,tie] = tiedrank(data);

% compute the anova SS's on the ranks:
[p1,p2] = anovan(R,grp,'display','off','model','interaction');
nvar = size(p2,1)-2;
for l=1:nvar
df(l) = p2{l+1,3};
SS(l) = p2{l+1,2};
end
% total MS
MS = sum(SS)/sum(df);
% compute ratios
for l=1:nvar-1
H(l) = SS(l)/MS;
end

% these are compared to a chi-2-distribution with df deg of freedom:
for l=1:nvar-1
p(l) = 1-chi2cdf(H(l),df(l));
end
]

Any ideas ?
Thanks in advance
 | 
Pages: 1
Prev: Problem using Kalman Filter estimation
Next: please change my E-mail Id