statcond questions - F stat permutation

1 view (last 30 days)
Xiaoyun Li
Xiaoyun Li on 10 Oct 2012
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.
Evelyn

Answers (0)

Categories

Find more on EEG/MEG/ECoG in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!