[R] Mixed model for count data?

Chun (Jimmie) Ye yechun at broadinstitute.org
Tue Jun 21 01:23:54 CEST 2011


Hi all,

I have a rather peculiar dataset that I'm not sure how to model
properly. This is data from an instrument that measures the size of
particles but instead of giving a continuous value, it generates a
"histogram" of the counts for a particular bin size. So the data looks
like this:

                              Condition 1
           Condition 2
Dimension  Rep1.A Rep1.B Rep2.A Rep2.B   Rep1.A Rep1.B Rep2.A Rep2.B
2                 5            6           7            8
40         35        33         31
2.1              7            8           4            5            30
        30        31         29
2.2              10          11         10          12          20
    18        21         20
2.3              50          45         44          39          5
      8         7           7
2.4              80          90         75          77          8
      10       3           5
2.5              30          22         31          35          10
      5        7           9
...
50               0            0            0            0           0
          1          0          0

There are two biological replicates and two technical replicates for
each of the conditions of interest. If these were continuous, I would
just fit a nested mixed model to the data. But instead, I have counts.
I was thinking that this might be suitable for modeling using a
multinomial response. I could also just compare the distribution using
a KS test but I'm not sure how to incorporate the replicate
information (except for averaging). Finally, I thought perhaps I could
estimate directly from the histogram the variance components from the
data. Any thoughts is greatly appreciated!

~Jimmie



More information about the R-help mailing list