kleine p-Werte (rank sum test) < Matlab < Mathe-Software < Mathe < Vorhilfe
|
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
Hallo.
Die Ausgangslage ist folgende: Ich habe in einer Versuchsreihe Probanden pseudorandomisiert mehrere Konditionen gezeigt (zB einmal einen roten Punkt und einmal einen blauen) und dann die Gehirnaktivität dazu über die Zeit gemessen. Demnach habe ich für jede Kondition und jeden Zeitschritt x Datenpunkte- mehrere Datenpunkte, weil eine große Fläche gleichzeitig gemessen wird.
Um blau vs rot zu vergleichen, nutze ich einen zweiseitigen Wilcoxon rank sum test in matlab- eine Normalverteilung liegt oft nicht vor.
Dabei werden mir aber selbst, wenn die Mittelwerte der jeweilgen Populationen sehr nah beieinander liegen, höchst signifikante Unterschiede ausgegeben mit Werten für p von zB 1*10e-27. Das passiert auch noch in den ersten Zeitschritten, in denen die Probanden noch gar nichts sehen. Zwar ist die Methode einigermaßen stör- bzw rauschanfällig, aber signifikante Unterscheide- zumal in diesen Dimensionen- sollten eigentlich nicht auftreten.
Die Frage ist also, wo der Fehler liegt (im code, falscher Test gar, ...).
Ich füge zum besseren Verständnis einen plot der beiden Mittelwerte an. Sterne oberhalb der Kurve signalisieren signifikante Unterschiede laut rank sum test, errorbars zeigen die einfache Standardabweichung.
[Dateianhang nicht öffentlich]
hier noch die leicht gekürzten Zeile(n) mit dem Test:
>> %size(datamatrix) = pixels - frames - conditions
>> f = 20; % frame (time)
>> conditionA = 1; % red
>> conditionB = 2; % blue
>> A=nanmean(squeeze(datamatrix(:,f,conditionA)),1);
>> B=nanmean(squeeze(datamatrix(:,f,conditionB)),1);
>> [h,p,ci,stats]=ttest(A,B,[],[],1);
Ich bedanke mich hier schon mal für eure Hilfe!
Dateianhänge:
Anhang Nr. 1 (Typ: jpg) [nicht öffentlich]
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 19:20 Di 07.01.2014 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
