[R] Blocking and Nesting ANOVA
Bob O'Hara
rn|@boh @end|ng |rom gm@||@com
Tue Jul 2 15:46:25 CEST 2019
No, it isn't correct I'm afraid. From what you describe, Population
isn't nested within Sex or Maturity. Population 1 can contain both
male and female trees, as well as trees of different ages.
A simpler model would be
lm (PAC ~ Season + Population + Sex + Maturity)
So that Sex and Maturity have the same effect in each population. If
you think their effects vary between populations, then you would want
lm (PAC ~ Season + Population*Sex + Population*Maturity)
which is the same as lm (PAC ~ Season + Population*(Sex + Maturity)).
Bob
On Tue, 2 Jul 2019 at 15:27, Eleftheria Dalmaris <edalmaris using gmail.com> wrote:
>
> Dear all,
>
> I have a blocking and nested design.
> I have Season (3 levels; three different collection periods). That's my
> Block.
> I have three Populations that I sampled from (Plant material)
> From each population, I sampled 30 trees.
> Each Tree is either Female or Male or Unidentified (Sex factor) and
> Each tree is either Juvenile or Mature or Old (Maturity factor)
>
> Maturity factor and Sex factor is not balanced within the population.
>
> So, Population is nested within Sex and Maturity
>
> I measure a chemical from the leaves of those trees with the name of PAC.
>
> Is the way that I write the code correct?
> lm (PAC ~ Season + Population+ Sex/Population + Maturity/Population)
>
> Thanks a lot,
>
> Eleftheria
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help using r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
--
Bob O'Hara
Institutt for matematiske fag
NTNU
7491 Trondheim
Norway
Mobile: +47 915 54 416
Journal of Negative Results - EEB: www.jnr-eeb.org
More information about the R-help
mailing list