Konfidenzintervalle einer Anova (Kontrast eines Faktors)

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

Localhorst

Forum-Fortgeschrittener


	Beiträge: 86

	Anmeldedatum: 18.02.08

	Wohnort: ---

	Version: ---

Verfasst am: 27.08.2015, 07:52 Titel: Konfidenzintervalle einer Anova (Kontrast eines Faktors)

Hallo Forum,

ich habe eine Anova und es interessieren mich die Konfidenzintervalle des Kontrastes des Factors "factor_area_numeric" gegen alle anderen Faktoren.

Ich habe die Konfidenzintervalle bzgl. der Differenz der Means bereits berechnet, aber das kommt glaube ich eher einem t-test nahe, da ich dann natürlich nicht die Faktoren einberechne, die ich in der Anova eigentlich integriert habe.

Ich habe aus meinem Code mal ein lauffähiges Script extrahiert. Evtl. kann mir einer von euch helfen. Wie in den Kommentaren zu sehen, habe ich mich schon an verschiedenen Anleitungen orientiert, da ich hier wirklich an die Grenzen meiner Mathematik stoße. Habt also bitte etwas Geduld mit mir Smile

Grüße
Localhorst

Code:

i=1;
combined_zscores=[0.306146153846154;0.198415384615385;0.256323076923077;0.509276923076923;0.234315384615385;0.291061538461539;0.442869230769231;0.550761538461539;0.347815384615385;0.149115384615385;0.243084615384615;0.303692307692308;0.431038461538462;0.504846153846154;0.496653846153846;0.666969230769231;0.679646153846154;0.346553846153846;0.236461538461538;0.347730769230769;0.295923076923077;0.272715384615385;0.0740384615384615;0.150015384615385;0.198415384615385;0.224330769230769;0.286615384615385;0.149115384615385;0.284753846153846;0.184323076923077;0.0447538461538462;0.158123076923077;0.133823076923077;0.129000000000000;0.274061538461538;0.190576923076923;0.225392307692308;0.190353846153846;0.236461538461538;0.302630769230769];
factor_area_numeric=[1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2];
factor_specimen=[3 3 3 3 3 2 2 2 5 5 5 4 4 4 4 1 1 6 6 6 6 3 3 3 3 2 5 5 4 6 6 6 6 6 6 6 6 6 6 6];
factor_age=[49 49 49 49 49 31 31 31 24 24 24 55 55 55 55 57 57 39 39 39 39 49 49 49 49 31 24 24 55 39 39 39 39 39 39 39 39 39 39 39];
factor_ethn=[2 2 2 2 2 3 3 3 1 1 1 3 3 3 3 3 3 1 1 1 1 2 2 2 2 3 1 1 3 1 1 1 1 1 1 1 1 1 1 1];

factors={factor_area_numeric factor_specimen factor_age factor_ethn};
varnames = {'Area';'Specimen';'Age';'ethn'};

% Anova => factor area == factor of interest
% factors specimen, age, ethn included as confounding factors
[p,tab,stats] = anovan(combined_zscores(:,i),factors,'varnames',varnames,'continuous',3,'display','off');

% calculation of eta2 effectsize (accourding to
% http://jalt.org/test/PDF/Brown28.pdf)
ss_total=tab{7,2};
ss_between_group_area=tab{2,2};
Reference_Anovan_eta2(i)=ss_between_group_area/ss_total;

%Difference between mean (accourding to http://onlinestatbook.com/2/estimation/difference_means.html)

% A confidence interval on the difference between means is computed using the following formula:
% Lower Limit = M1 - M2 -(tCL)(S m1-m2)
% Upper Limit = M1 - M2 +(tCL)(S m1-m2)
%
% where M1 - M2 is the difference between sample means,
% tCL is the t for the desired level of confidence, and S m1-m2
% is the estimated standard error of the difference between sample means.

mse=(var(combined_zscores(factor_area_numeric==1,i))+ var(combined_zscores(factor_area_numeric==2,i)))/2;
sm1m2=sqrt((2*mse)/size(combined_zscores,1));
% t for 95% with 48 degree of freedom confidence interval 2.011
% (http://onlinestatbook.com/2/calculators/inverse_t_dist.html)
Reference_Anovan_CI_l(i)=(mean(combined_zscores(factor_area_numeric==1,i))- mean(combined_zscores(factor_area_numeric==2,i)))-2.011*sm1m2;
Reference_Anovan_CI_h(i)=(mean(combined_zscores(factor_area_numeric==1,i))- mean(combined_zscores(factor_area_numeric==2,i)))+2.011*sm1m2;
Reference_Anovan_diff_mean(i)=mean(combined_zscores(factor_area_numeric==1,i))-mean(combined_zscores(factor_area_numeric==2,i));

Funktion ohne Link?


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2024 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.