コードをここにコピペしてたけど
半角の<があるとそれ以降が改行のところまで消し飛んでた。
ってことは過去のも同じか。。。
保存用に残してたのに意味ナス。
2015-02-23
permutation testにおける多重比較補正
備忘録
リサンプリングによって得られた各指標のp値を算出。
最小のp値を抽出。
リサンプルを繰り返す。
最小p値の分布において実際に観測されたp値がどのぐらいの確率で出現するかを計算。
以下は対応のない2群比較にのみ対応。
これであってんのか自信なし。
いずれにしろ変数の数が増えるとやっぱ厳しい。
function Pcorrected=rbmt_ttst2(data,g1,g2,k)
% resampling-based multiple test
% input
% data: sub x var.
% g1: subjects of gourp 1, eg. 1:10
% g2: subjects of gourp 1, eg. 11:20
% k: number of permutation
s=size(data);
p1=zeros([s(2) 1]);
for ii=1:s(2)
[~, P]=ttest(data(g1,ii),data(g2,ii));
p1(ii)=P;
end;
Pcorrected=zeros([s(2) 1]);
for jj=1:k,
idx=randperm(s(1));
minP=1;
for ii=1:s(2),
[~, P]=ttest(data(idx(g1),ii),data(idx(g2),ii));
if P<minP, minP=P; end;
end;
Pcorrected=Pcorrected+(p1>minP);
end;
Pcorrected=Pcorrected/k;
なんかコピペミスしてたので,修正
リサンプリングによって得られた各指標のp値を算出。
最小のp値を抽出。
リサンプルを繰り返す。
最小p値の分布において実際に観測されたp値がどのぐらいの確率で出現するかを計算。
以下は対応のない2群比較にのみ対応。
これであってんのか自信なし。
いずれにしろ変数の数が増えるとやっぱ厳しい。
function Pcorrected=rbmt_ttst2(data,g1,g2,k)
% resampling-based multiple test
% input
% data: sub x var.
% g1: subjects of gourp 1, eg. 1:10
% g2: subjects of gourp 1, eg. 11:20
% k: number of permutation
s=size(data);
p1=zeros([s(2) 1]);
for ii=1:s(2)
[~, P]=ttest(data(g1,ii),data(g2,ii));
p1(ii)=P;
end;
Pcorrected=zeros([s(2) 1]);
for jj=1:k,
idx=randperm(s(1));
minP=1;
for ii=1:s(2),
[~, P]=ttest(data(idx(g1),ii),data(idx(g2),ii));
if P<minP, minP=P; end;
end;
Pcorrected=Pcorrected+(p1>minP);
end;
Pcorrected=Pcorrected/k;
なんかコピペミスしてたので,修正
2015-02-12
マルチカーネルSVM
Multimodal analysis of functional and structural disconnection in Alzheimer's disease using multiple kernel SVM.
Dyrba M, Grothe M, Kirste T, Teipel SJ.
Hum Brain Mapp. 2015 Feb 9. doi: 10.1002/hbm.22759.
これもマルチカーネル。
が,マルチサイトのデータがほとんど出てこない。
多分みんな苦戦してるんだろうけど,
そこをクリアしないと使えない。
Dyrba M, Grothe M, Kirste T, Teipel SJ.
Hum Brain Mapp. 2015 Feb 9. doi: 10.1002/hbm.22759.
これもマルチカーネル。
が,マルチサイトのデータがほとんど出てこない。
多分みんな苦戦してるんだろうけど,
そこをクリアしないと使えない。
2015-02-05
軸索定位をMRIで
Functional Consequences of Neurite Orientation Dispersion and Density in Humans across the Adult Lifespan.
J Neurosci. 2015 Jan 28;35(4):1753-62
Authors: Nazeri A, Chakravarty MM, Rotenberg DJ, Rajji TK, Rathi Y, Michailovich OV, Voineskos AN
軸索定位の散らばりと密度をMRIで測定できるらしい。
J Neurosci. 2015 Jan 28;35(4):1753-62
Authors: Nazeri A, Chakravarty MM, Rotenberg DJ, Rajji TK, Rathi Y, Michailovich OV, Voineskos AN
軸索定位の散らばりと密度をMRIで測定できるらしい。
登録:
投稿 (Atom)