[R] [OFF] Nested or not nested, this is the question.

Ronaldo Reis Jr. chrysopa at insecta.ufv.br
Wed Apr 9 21:07:57 CEST 2003


Hi,

sorry by this off.

I'm still try to understand nested design.

I have the follow example (fiction):

I have 12 plots in 4 sizes in 3 replicates (4*3 = 12)
In each plot I put 2 species (A and B) to reproduce.
After a period I make samples in each board and count the number of 
individuals total (tot) and individuals A and B (nsp). Others individuals 
excepts A and B are in total of individuals.

This make a dataset with the 24 lines and not 12. Its smell pseudoreplication 
in a nested design, OK?

I need to know:

the species are different in proportion?

the size affect the species's proportion?

existe interaction between size and species?

I make the analysis.

> m.lme <- lme(nsp/tot~size*specie,random=~1|size/specie)
> anova(m.lme)
            numDF denDF  F-value p-value
(Intercept)     1    16 374.7121  <.0001
size            1     2  37.8683  0.0254
specie          1     2  18.2036  0.0508
size:specie     1     2   9.3203  0.0926
> 

This is the correct mean to make this analysis?

or

> m.lme <- lme(nsp/tot~size*specie,random=~1|plot/specie)
> anova(m.lme)
            numDF denDF  F-value p-value
(Intercept)     1    10 579.8853  <.0001
size            1    10  58.6030  <.0001
specie          1    10  59.5235  <.0001
size:specie     1    10  30.4760   3e-04
> 

or neither?

I know about the distribution (binomial in this case), but I try to understand 
the nested design.

Thanks for any help.

The dataset is:

   plot size specie nsp tot
1     1   10      A   2  20
2     1   10      B   6  20
3     5   10      A   3  20
4     5   10      B   5  20
5     9   10      A   1  20
6     9   10      B   4  20
7     2   20      A   5  20
8     2   20      B   8  20
9     6   20      A   6  20
10    6   20      B   9  20
11   10   20      A   4  20
12   10   20      B   6  20
13    3   30      A   8  20
14    3   30      B   9  20
15    7   30      A   9  20
16    7   30      B  10  20
17   11   30      A   7  20
18   11   30      B   8  20
19    4   40      A  10  20
20    4   40      B   9  20
21    8   40      A   9  20
22    8   40      B  10  20
23   12   40      A   9  20
24   12   40      B   9  20
-- 
Great minds run in great circles.
--
|   // | \\   [*****************************][*******************]
|| ( õ   õ )  [Ronaldo Reis Júnior          ][PentiumIII-600     ]
|      V      [UFV/DBA-Entomologia          ][HD: 30 + 10 Gb     ]
||  /     \   [36571-000 Viçosa - MG        ][RAM: 128 Mb        ]
|  /(.''`.)\  [Fone: 31-3899-2532           ][Video: SiS620-8Mb  ]
||/(: :'  :)\ [chrysopa at insecta.ufv.br      ][Modem: Pctel-onboar]
|/ (`. `'` ) \[ICQ#: 5692561                ][Kernel: 2.4.18     ]
||  ( `-  )   [*****************************][*******************]
||| _/   \_Powered by GNU/Debian W/Sarge D+ || Lxuser#: 205366



More information about the R-help mailing list