I have data that has been subjected to one of two treatments (solution added or not) and is in a random block design. I have 8 blocks (numbered 1-8) but in the aov, df=1. R didn't recognise there were 8, so I lettered them (B1, B2 etc). Now there is no F value or p value in the two-way aov, although now the df=7.
Originally:
soilaov <- aov(Si.arc.leaf~Block*Treatment, data=soilbd)
summary(soilaov)
Output
`Df` `Sum Sq` `Mean Sq` `F value` `Pr(>F)`
Block `1` `5.174e-05` `5.174e-05` `3.941` `0.0705 .`
Treatment `1` `2.526e-05` `2.526e-05` `1.924` `0.1907`
Block:Treatment `1` `7.310e-06` `7.310e-06` `0.557` `0.4699`
Residuals `12` `1.576e-04` `1.313e-05`
Now:
soilaov <- aov(Si.arc.leaf~Block*Treatment, data=soilbd)
summary(soilaov)
Output
`Df` `Sum Sq` `Mean Sq`
Block `7` `1.647e-04` `2.352e-05`
Treatment `1` `2.526e-05` `2.526e-05`
Block:Treatment `7` `5.193e-05` `7.419e-06`