statcond questions - F stat permutation

2 次查看（过去 30 天）
Xiaoyun Li
Xiaoyun Li 2012-10-10
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

请先登录，再回答此问题。

回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 EEG/MEG/ECoG 的更多信息

标签

Community Treasure Hunt

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

Start Hunting!

Translated by