Hi all,

I am an honours student and very new to R. Any help as to obtaining Std
Deviations for each interaction term in my LMER model would be greatly
appreciated.

This is a subset of my original data, and is in a table as follows:

id

genus

species1

species

pop

treat

bed

range

week

w

1

1

1

Rytidosperma penicillatum

31

1

11

2

6

5.5

3

1

1

Rytidosperma penicillatum

31

1

11

2

6

6.5

4

1

2

Rytidosperma nudiflorum

29

1

11

1

6

2.3

5

1

2

Rytidosperma nudiflorum

29

1

11

1

6

2.5

6

1

1

Rytidosperma penicillatum

26

1

11

2

6

5.5

7

1

1

Rytidosperma penicillatum

26

1

11

2

6

6

8

2

3

Senecio extensus

14

1

11

1

6

1.5

9

2

4

Senecio gunii

13

1

11

2

6

1.7

10

2

4

Senecio gunii

11

1

11

2

6

1

11

2

3

Senecio extensus

14

1

11

1

6

3

857

1

1

Rytidosperma penicillatum

26

2

3

2

6

2.5

863

13

16

Argyrotegium fordianum

89

2

3

2

6

1.5

864

13

15

Argyrotegium nitidulum

76

2

3

1

6

0.6

868

13

16

Argyrotegium fordianum

86

2

3

2

6

1.8

869

13

15

Argyrotegium nitidulum

76

2

3

1

6

0.6

Where w=width

            > str(wid)

'data.frame':            5741 obs. of  10 variables:

 $ id      : Factor w/ 1528 levels "1","3","4","5",..: 1 2 3 4 5 6 7 8 9 10
...

 $ genus   : Factor w/ 8 levels "1","2","6","7",..: 1 1 1 1 1 1 2 2 2 2 ...

 $ species1: Factor w/ 16 levels "1","2","3","4",..: 1 1 2 2 1 1 3 4 4 3 ...

 $ species : Factor w/ 25 levels "Argyrotegium fordianum",..: 18 18 17 17
18 18 21 22 22 21 ...

 $ pop     : Factor w/ 42 levels "1","10","11",..: 14 14 13 13 11 11 5 4 3
5 ...

 $ treat   : Factor w/ 2 levels "1","2": 1 1 1 1 1 1 1 1 1 1 ...

 $ bed     : Factor w/ 8 levels "3","4","5","6",..: 8 8 8 8 8 8 8 8 8 8 ...

 $ range   : Factor w/ 2 levels "1","2": 2 2 1 1 2 2 1 2 2 1 ...

 $ week    : Factor w/ 4 levels "6","12","18",..: 1 1 1 1 1 1 1 1 1 1 ...

 $ w       : num  5.5 6.5 2.3 2.5 5.5 6 1.5 1.7 1 3 ...



My LMER model is

>lm<-lmer(w~treat*range+(1+treat|bed)+(1+treat|pop)+(1+treat*range|genus)+(1+treat*week|id),
data=wid)

 and my summary(lm) output is standard.

            However I want to investigate the interactions for all the
different levels in some of my factors (e.g. How does treatment affect
Genus 1, compared to Genus 2?). I found this code, which seems to do the
job nicely:

>lm.results<-list(ran=ranef(lm), fix=fixef(lm, vc=VarCorr(lm)))

>lm.results



But the output of this does not (to my knowledge) provide Std. Deviation
for the various interactions (this is just the genus section of the output):



$genus

            (Intercept)            treat2            range2
treat2:range2

1            0.84                        -0.07            1.86
-0.38

2            0.78                        -0.06            -0.75
0.24

6            -0.91                        0.08            -0.58
0.07

7            -1.34                        0.11            -0.65
0.05

10            -0.02                        0.00            -0.12
0.03

11            -0.89                        0.07            -0.89
0.14

12            2.42                        -0.20            1.65
-0.20

13            -0.87                        0.07            -0.52
0.06



I feel as though there must be something obvious I am missing, but for all
the googling and R-help searching, I have been unable to locate it.


Any help would be greatly appreciated,
Em

	[[alternative HTML version deleted]]

