Maar het is dus correct dat ik een one-sample t-test neem met alle data van hetzelfde gen als basis t-stat en niet een two-sample t-test waarbij de data van ALL en AML zijn gescheiden?quote:Op zondag 21 februari 2016 23:10 schreef Anoonumos het volgende:
De nulhypothese is dat gen 1 en gen 2 niet differentially expressed zijn.
Dus onder de nulhypothese maakt het niet uit als je de groepen door elkaar gooit. Dat doe je dus met die permutaties.
Maar blijkbaar is dat niet gerechtvaardigd want t (de echte waarneming) is veel groter dan de meeste t* (uitkomsten van de permutaties), met andere woorden de p-value is klein, dus je concludeert dat ze wel differentially expressed zijn.
Ja boeie, ik snap je.quote:Op zondag 21 februari 2016 23:15 schreef Anoonumos het volgende:
[..]
Oke deze zin klopt biologisch niet geloof ik
Of een gen differentially expressed is op sample 1 en 2
Nee moet two-sample zijn lijkt mequote:Op zondag 21 februari 2016 23:17 schreef Rezania het volgende:
[..]
Maar het is dus correct dat ik een one-sample t-test neem met alle data van hetzelfde gen als basis t-stat en niet een two-sample t-test waarbij de data van ALL en AML zijn gescheiden?
Oké, aangepast. Nu heb ik dit:quote:Op zondag 21 februari 2016 23:22 schreef Anoonumos het volgende:
[..]
Nee moet two-sample zijn lijkt me
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 | % calculate the original t-statistic for CCND3 matlab_tstat = t.tstat % select the number of permutation nperm = 1000 for i = 1 : nperm % generate a new label vector (for now all labels are zero -> ALL) l = zeros(length(labels),1); % randomly generate 11 indices between 1 and 38 r = randperm(38,11); % assign ones to the random 11 indices (random AML samples) l(r) = 1; % find the indices of the ALL and ALL samples randALL = find(l == 0); randAML = find(l == 1); % extract the random samples from the original data rALL = data(1042, randALL); rAML = data(1042, randAML); % calculate the t-statistic between rALL and rAML for iteration i [~,~,~,g] = ttest2(rALL,rAML) rand_tstat(i) = g.tstat end % plot the histogram of all the t-statistics you calculated figure, hist(rand_tstat) % draw a line at the original t-statistic value line([matlab_tstat matlab_tstat], [0 200], ... 'LineWidth',2,'Color','r', 'LineStyle', '--') % calculate the permutation-based p-value perm_pvalue = (find(rand_tstat >= t.tstat))/nperm |
1 2 3 | perm_pvalue = Empty matrix: 1-by-0 |
Zou geen verschil moeten maken qua resultaten, maar kan het voor de zekerheid wel bijzetten.quote:Op zondag 21 februari 2016 23:42 schreef Anoonumos het volgende:
Je moet denk ik wel absolute waardes van t bekijken. Tweezijdige test dit dat.
En ik ben al m'n memory leaks kwijt.quote:Op dinsdag 23 februari 2016 14:57 schreef Rezania het volgende:
Gewoon geen één error gekregen met Matlab vandaag.
quote:Op woensdag 24 februari 2016 04:16 schreef wimjongil het volgende:
[..]
En ik ben al m'n memory leaks kwijt.
C ja. In het eerste jaar van informatica krijg je daarin de meest fundamentele vakken. Omdat je overal op moet letten leer je er wel netjes door programmeren.quote:
Dat heeft voordelen. Wij moeten al die shit maar zelf zien uit te vogelen.quote:Op woensdag 24 februari 2016 17:43 schreef wimjongil het volgende:
[..]
C ja. In het eerste jaar van informatica krijg je daarin de meest fundamentele vakken. Omdat je overal op moet letten leer je er wel netjes door programmeren.
Ik kies juist groepen waar ik eindelijk weer op het lab mag staan.quote:Op woensdag 24 februari 2016 22:41 schreef Anoonumos het volgende:
Ik koos de enige groep waar je geen experimenten hoefde te doen.
+1quote:Op woensdag 24 februari 2016 22:58 schreef Rezania het volgende:
http://www.observantonlin(...)m-nu-student-te-zijn
Heel veel klaagzang over prestatiedruk, maar wat me opvalt is dat er alleen maar studenten worden genoemd die van die massastudies volgen.
Je lijkt een enorme Lijenaar, waarom kom je dan naar Delft?quote:Op woensdag 24 februari 2016 22:44 schreef Rezania het volgende:
[..]
Ik kies juist groepen waar ik eindelijk weer op het lab mag staan.
Omdat Leiden meh is. Hoe kom je erbij dat ik een Lijenaar ben?quote:Op woensdag 24 februari 2016 23:35 schreef Bravebart het volgende:
[..]
Je lijkt een enorme Lijenaar, waarom kom je dan naar Delft?
Qua dingen waar je het hier over hebt; genetisch modificeren, genen, plasmiden enzo. Heel Leids imhoquote:Op woensdag 24 februari 2016 23:35 schreef Rezania het volgende:
[..]
Omdat Leiden meh is. Hoe kom je erbij dat ik een Lijenaar ben?
IMB doet bijna niks anders. Het ligt eraan wat je ermee doet, niet waar je het doet.quote:Op woensdag 24 februari 2016 23:43 schreef Bravebart het volgende:
[..]
Qua dingen waar je het hier over hebt; genetisch modificeren, genen, plasmiden enzo. Heel Leids imho
Klopt, IMB is ook een stukje Leiden dat verdwaald is in Delftquote:Op woensdag 24 februari 2016 23:43 schreef Rezania het volgende:
[..]
IMB doet bijna niks anders. Het ligt eraan wat je ermee doet, niet waar je het doet.
Forum Opties | |
---|---|
Forumhop: | |
Hop naar: |