statcond questions - F stat permutation

2 次查看（过去 30 天）

Xiaoyun Li 2012-10-10

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/50373-statcond-questions-f-stat-permutation

I need to analyse the significant activation electodes in the 129-channel EEG system in selected time epoch. I have two groups (H & L) with 12 subjects per each and 3 conditions (c, t, n), the time epoch includes 1200 time points. I hope to run the permutation analysis. Here is part of my script for the statistics.

H_c=zeros(129,12); H_t=zeros(129,12); H_n=zeros(129,12); L_c=zeros(129,12); L_t=zeros(129,12); L_n=zeros(129,12);

for i=1:size(Hi,2) for j=1:129 tmp1=squeeze(squeeze(mean(HC(s1:s2,j,i)))); tmp2=squeeze(squeeze(mean(HT(s1:s2,j,i)))); tmp3=squeeze(squeeze(mean(HN(s1:s2,j,i)))); H_c(j,i)=tmp1; H_t(j,i)=tmp2; H_n(j,i)=tmp3;

    end;
end;

for i=1:size(Lo,2) for j=1:129 tmp1=squeeze(squeeze(mean(LC(s1:s2,j,i)))); tmp2=squeeze(squeeze(mean(LT(s1:s2,j,i)))); tmp3=squeeze(squeeze(mean(LN(s1:s2,j,i)))); L_c(j,i)=tmp1; L_t(j,i)=tmp2; L_n(j,i)=tmp3; end; end;

A={H_c H_t H_n; L_c L_t L_n};

[F df p]=statcond(A,'paired','off','mode','perm','naccu',1000);

for i=1:3 for j=1:129 if p{i}(j)>0.05 F{i}(j)=0; end; end; end;

out = zeros(1200,129); for j=1:3 for k=1:129 out(1+(j-1)*100:j*100,k)=-F{j}(k); end; end;

outname = [dir '\epoch_' num2str(epoch(1,1)) '_' num2str(epoch(1,2)) '.asc'];

However, when I run the script,there is an error showed out.

??? Error using ==> cell.strmatch at 24 FLAG must be a string.

Error in ==> finputcheck>fieldtest at 201 if isempty(strmatch(lower(tmpval), lower(fieldval), 'exact'))

Error in ==> finputcheck at 122 res = fieldtest( fieldlist{index, NAME}, fieldlist{index, TYPE}, ...

Error in ==> statcond at 139 g = finputcheck( varargin, { 'naccu' 'integer' [1 Inf] 200;

Error in ==> comflaser_electrode_test at 125 [F df p]=statcond(A,'paired','off','mode','perm','naccu',1000);

I am not sure where is the problem and how to sout it out. It will be much grateful if anyone can guide me this. Many thanks.