Hi,
The plant and Treatment are aliased (if you know the value of plant you know
the value of treatment) so you shouldn't add them together in the same
model. That being said when you enter treatment first and then plant you are
adding the extra information from plant as the second effect and in the
second model when you try to enter treatment as the second effect there's
nothing more to fit. If you look closely at the two model you will notice
that the total variance and the total degree of freedom are the same. Also
when you use *, what you are fitting is plant + treatment +
plant:treatment but you don't see the interaction in the model because it is
useless in the model.
Hope it helps,
Charles-?douard
Treatment Qn1 Qn2 Qn3 Qc1 Qc3 Qc2 Mn3 Mn2 Mn1 Mc2 Mc3 Mc1
nonchilled 7 7 7 0 0 0 7 7 7 0 0 0
chilled 0 0 0 7 7 7 0 0 0 7 7 7
-----Message d'origine-----
De : R-help <r-help-bounces at r-project.org> De la part de konrad via R-help
Envoy? : 17 mars 2022 10:17
? : r-help at r-project.org
Objet : [R] Possible bug
Dear all,
I recently conducted aov calculations of datasets::CO2.
When I'm using the following code:
'summary(aov(CO2$uptake~CO2$Treatment*CO2$Plant))'. Everything is fine.
I get a result for Plant and Treatment.
However, if using: 'summary(aov(CO2$uptake~CO2$Plant*CO2$Treatment))'
'Treatment' is ignored.
I'm using R 3.6 on Linux Mint. Any idea why this happens?
Regards,
Konrad
______________________________________________
<mailto:R-help at r-project.org> R-help at r-project.org mailing list -- To
UNSUBSCRIBE and more, see <https://stat.ethz.ch/mailman/listinfo/r-help>
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
<http://www.R-project.org/posting-guide.html>
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.