[R-sig-ME] plot glmmPQL

Mario V. Balzan m.balzan at sssup.it
Thu Apr 26 13:15:05 CEST 2012


Dear all,

I am investigating the influence of some experimental treatments
(A,B,C,...) on fecundity of a lepidoptera. These data have been
collected over time from different experimental plots (A1, A2, ... ).
I am using glmmPQL function with quasipoisson errors since these are
count data and with overdispersion.

the model which I came up with is as follows:
glmmPQL(Fecundity~Treatment,random=~Time|Plot, family="quasipoisson")

I would like to plot a curve of fecundity against time for each
treatment from this model. Is that possible? Would appreciate if
someone could help me out with this.

this is a sample of my data:
  Time. Fecundity Plot Treatment
      2        40   A1         A
      4        10   A1         A
      6         0   A1         A
      8         0   A1         A
      2        21   A2         A
      4         0   A2         A
      6         0   A2         A
      2         9   A3         A
      4         0   A3         A
    2        23   A4         A
     4         0   A4         A
     6         2   A4         A
     2         0   B1         B
    4         0   B1         B
    6        10   B1         B
     8         0   B1         B
   11         0   B1         B
    2         1   B2         B
    4         2   B2         B
    6         0   B2         B
    8        25   B2         B
   11         6   B2         B
    13         4   B2         B
    15        10   B2         B
    18         8   B2         B
     2         5   B3         B
     4         6   B3         B
     6        20   B3         B
     8         2   B3         B
    11         2   B3         B
     2        23   C1         C
     4         8   C1         C
     6         4   C1         C
     8         0   C1         C
    10         2   C1         C
   13         0   C1         C
     2        15   C2         C
     4        15   C2         C
     6         4   C2         C
     8         0   C2         C
    10         0   C2         C
    13         0   C2         C
     2        65   C3         C
     4         2   C3         C
     6         8   C3         C
     8        10   C3         C
    10         6   C3         C
   13         2   C3         C
...


Thanks in advance,
Mario



More information about the R-sig-mixed-models mailing list