[R-sig-ME] How to handle measurement uncertainities with lmer
Frömmel, Ulrike
Ulrike.Froemmel at HS-Lausitz.de
Wed Aug 8 21:12:56 CEST 2012
Dear list members,
a result set obtained from an automated video microscope contains
adhesion density values of a particular bacterial isolate per cell
surface area (germs/mm^2). Because the device determines this density
several times, the result set also contains an associated measurement
uncertainity, given as coefficient of variation Cv = sd/µ (standard deviation / mean).
In addition, each isolate was determined repeatedly on the same plate,
and also on different plates.
Independently from these experiments, the presence or absence of around 40
genes had been determined for each of the 316 bacterial isolates.
To investigate how much a particular gene influences bacterial adhesion,
the following lmer model was used:
lmer(log(`Germs/mm^2`+1) ~ -1 + GeneA + GeneB + .... GeneZ + (1|Isolate) + (1|Plate), data=df)
> str(df):
'data.frame': 4101 obs. of 67 variables:
$ Adhesion : num 5.67 5.55 5.33 6.85 6.94 ...
$ Isolate : Factor w/ 316 levels "2164","2166",..: 183 183 183 12 12 107 107 107 246 76 ...
$ Plate : Factor w/ 59 levels "5637 BP01 m1",..: 1 1 1 1 1 1 1 1 1 1 ...
$ Germs/mm^2 : num 288 256 206 946 1031 ...
$ Germs_cv : num 48.49 5.51 4.56 9.07 6.78 ...
$ GeneA : Factor w/ 2 levels "-","+": 1 1 1 1 1 1 1 1 1 1 ...
$ GeneB : Factor w/ 2 levels "-","+": 1 1 1 1 1 1 1 1 1 1 ...
$ GeneC : Factor w/ 2 levels "-","+": 2 2 2 1 1 1 1 1 1 1 ...
[...] <- ~ 40 Genes plus some other data)
The log transformation was applied because the residual vs. fitted diagram
strongly indicated heteroscedasticity. With the log transformation applied,
the density values are approximately normal distributed.
My questions are:
1/ How to feed the uncertainity measures into lmer?
2/ What would an effect of zero actually mean, considering that the
data had to be log transformed in order to cope with
heteroscedascity?
3/ How can I do the inverse transform of the effects?
Best regards
More information about the R-sig-mixed-models
mailing list