problem with the anova calculation

7 visualizzazioni (ultimi 30 giorni)
Michael
Michael il 6 Dic 2022
Commentato: Jeff Miller il 10 Dic 2022
I have a problem with the calculation of the ANOVA. What is wrong with my code?
C_t = [30; 30; 30; 30; 30; 30; 30; 120; 120; 120; 120; 120; 120; 120; 120];
pick = [0; 0.026; 0.0550; 0.120 ; 0.320 ; 0.62 ; 1.2; 0 ;0.03; 0.066; 0.110 ; 0.220; 0.4; 0.740; 1.73];
tbl = table(C_t,pick,'VariableNames',{'C_t','pick'});
tbl.C_t = categorical(tbl.C_t);
mdl = fitlm(tbl,'pick')
tbl = anova(mdl,'summary'

Risposta accettata

Jeff Miller
Jeff Miller il 6 Dic 2022
Just change one line to this
mdl = fitlm(tbl,'pick~C_t')
  4 Commenti
Michael
Michael il 10 Dic 2022
Ok, but how can i calculate leak of fit and purre error from this data?
Jeff Miller
Jeff Miller il 10 Dic 2022
You can't. This is not a software problem, it is a limitation of your design. With only two values of X (30, 120), a straight line goes through the two (X, avg Y) pairs perfectly, so your data provide no information about whether there would be any deviation from a straight line fit ("lack of fit") if you had a third X value.

Accedi per commentare.

Più risposte (1)

James Water Bird
James Water Bird il 6 Dic 2022
C_t = [30; 30; 30; 30; 30; 30; 30; 120; 120; 120; 120; 120; 120; 120; 120];
pick = [0; 0.026; 0.0550; 0.120 ; 0.320 ; 0.62 ; 1.2; 0 ;0.03; 0.066; 0.110 ; 0.220; 0.4; 0.740; 1.73];
% Use the table() function to create a table from the data
tbl = table(C_t,pick,'VariableNames',{'C_t','pick'});
% Convert the 'C_t' column to a categorical variable
tbl.C_t = categorical(tbl.C_t);
% Use the fitlm() function to fit a linear regression model to the data
mdl = fitlm(tbl,'pick')
% Use the anova() function to perform an analysis of variance on the model
tbl = anova(mdl,'summary')
  1 Commento
Michael
Michael il 6 Dic 2022
I understand that you have added descriptions, but this does not solve my problem.

Accedi per commentare.

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by