Nooit wat ik zoek in ieder geval. Over R is veel meer te vinden, in ieder geval in mijn straatje. Komt natuurlijk deels ook doordat matlab wat breder is en R meer toegespitst is op statistiek, maar daar gebruik ik matlab ook voornamelijk voor.quote:Op woensdag 17 februari 2016 02:01 schreef Bravebart het volgende:
[..]
Vent, er zijn 60.000 topics over MatLab op stack
quote:Op woensdag 17 februari 2016 12:08 schreef Anoonumos het volgende:
En dan was de bedenker van de t-verdeling nog wel een bier brouwende halve scheikundige.
Jij hebt het over het feit dat AMD niet drivers maakt voor de Linux maar in plaats hiervan de code geeft en de OS gemeenschap ermee laat knutselen? Ik kan me voorstellen dat dat allemaal wat minder soepel loopt gezien de korte cyclus van grafische kaarten. Ik heb begrepen dat 'oudere' AMD-kaarten prima werken maar dat je bij Linux vaak 2 generaties of zo achterloopt.quote:Op woensdag 17 februari 2016 01:31 schreef wimjongil het volgende:
[..]
Als je ooit iets anders dan Windows een kans wilt geven: blijf weg bij AMD.
Eens.quote:Die matlab documentatie is echt verschrikkelijk. Daarnaast is er niks op stackoverflow ofzo te vinden.
Dat is vaak veel te specifiek voor die ene gebruiker zijn ene toepassing. Soms heb je er wat aan, vaak heb je er niets aan.quote:[quote] Op woensdag 17 februari 2016 02:01 schreef Bravebart het volgende:
[..]
Vent, er zijn 60.000 topics over MatLab op stack
Omdat de eagles niemand van dienst zijn.quote:Op donderdag 18 februari 2016 21:52 schreef Anoonumos het volgende:
Waarom vlogen de eagles Frodo niet gewoon naar Mordor?
1 2 3 4 5 6 7 8 9 10 11 12 | % load the Golub dataset (make sure you are in the right directory) load('golub.mat') % set the index of gene Gdf5 (it is the 2058th row) gene_idx = 2058; % get the expression data of Gdf5 a = data(gene_idx,:); % obtain the labels of each class ALL = find(labels == 0); AML = find(labels == 1); % calculate the mean of CRH in the ALL and AML samples b = mean(a(ALL)); c = mean(a(AML)); |
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 = ??? % select the number of permutation nperm = ??? 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(gene_idx, randALL); rAML = data(gene_idx, randAML); % calculate the t-statistic between rALL and rAML for iteration i ??? = ??? rand_tstat(i) = ??? 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 = ??? |
1 | [~,~,~,tstat]=ttest2(rALL,rAML) |
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 % select the number of permutation nperm = 100 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(gene_idx, randALL); rAML = data(gene_idx, randAML); % calculate the t-statistic between rALL and rAML for iteration i [~,~,~,tstat] = ttest2(rALL,rAML) rand_tstat(i) = ??? 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 = ??? |
Ja maar, die varieert per i lijkt me. We krijgen het niet voor mekaar om die berekening van de tstat te koppelen aan rand_stat.quote:Op vrijdag 19 februari 2016 15:44 schreef Anoonumos het volgende:
Bij rand_tstat moet je de waarde van de t-statistic invullen
Oh, er moet nog iets voor? Dat wist ik niet. Ik deed gewoon rand_tstat=tstat.quote:Op vrijdag 19 februari 2016 16:11 schreef Anoonumos het volgende:
[~,~,~,tstat] = ttest2(rALL,rAML)
rand_tstat = tstat.tstat ( ) ?
[~,~,~,a] = ttest2(rALL,rAML)
rand_tstat = a.tstat
Staat dus nergens in de 'handleiding' die we hebben gekregen en ik zou echt niet weten hoe ik dit zou moeten Googlen.quote:Op vrijdag 19 februari 2016 16:53 schreef Bram_van_Loon het volgende:
Een moooie illustratie van het echte probleem voor studenten die met Matlab werken. Het gevonden krijgen hoe je in godsnaam die functie gebruikt.
Ingewikkeld is het niet maar zoek er maar eens naar.
Voor de duidelijkheid, ik ben het volstrekt met je eens en ik schat in dat het niet aan jou ligt aangezien ik dit zelf als student ook vaak genoeg heb meegemaakt.quote:Op vrijdag 19 februari 2016 16:55 schreef Rezania het volgende:
[..]
Staat dus nergens in de 'handleiding' die we hebben gekregen en ik zou echt niet weten hoe ik dit zou moeten Googlen.
Je moet hem dan ook integralen geven waar zijn rekenmachine op stuk loopt. Ik heb er nog wel een paar voor je. Moet je ze zelf natuurlijk wel kunnen oplossen, anders sta je in je hemd.quote:Op zaterdag 20 februari 2016 16:39 schreef Rezania het volgende:
Ik zat gister bijles te geven, wiskunde B 5 vwo, integralen. Hij moest 8/x˛ integreren van 1 naar 8. Pakt hij zijn rekenmachine, tikt wat in, en het antwoord kwam er zo uitrollen. Wtf. Gelijk rekenmachine afgepakt en gezegd dat hij vanaf nu alles handmatig moet doen.
Zeker met complexe getallen?quote:Op zaterdag 20 februari 2016 18:02 schreef Riparius het volgende:
[..]
Je moet hem dan ook integralen geven waar zijn rekenmachine op stuk loopt. Ik heb er nog wel een paar voor je. Moet je ze zelf natuurlijk wel kunnen oplossen, anders sta je in je hemd.
Het gaat wel, alleen het heeft allemaal valkuilen die nergens staan vermeld.quote:Op zaterdag 20 februari 2016 18:05 schreef b2py het volgende:
Ik moet volgend kwartiel met Matlab, aan de slag. Ben nu al verdrietig over het idee er aan..
Partiële integratie, integratie door substitutie, en integratie door breuksplitsing zijn keuzeonderwerpen voor wiskunde B... Hij hoeft die dus niet per se op te kunnen lossen...quote:Op zaterdag 20 februari 2016 18:02 schreef Riparius het volgende:
[..]
Je moet hem dan ook integralen geven waar zijn rekenmachine op stuk loopt. Ik heb er nog wel een paar voor je. Moet je ze zelf natuurlijk wel kunnen oplossen, anders sta je in je hemd.
Nee hoor. Hier heb je er eentje die Casio rekenmachines vast laat lopen.quote:
Forum Opties | |
---|---|
Forumhop: | |
Hop naar: |