Question: Dose-Response: Regressiion and ANOVA
0
14.8 years ago by
Arne.Muller@aventis.com620 wrote:
Dear All, I'm interested in your experience of analyzing dose-response experiments. I've an unbalanced design with 3 studies and within each study 4 to 5 doses of a drug (the doses do not exactmy the same in each study). Usually I run an anova per gene anova(lm(Intensity ~ Study + Dose + Study:Dose)) where Study and Dose are factors, and for the linear model I use a treatment contrast with the control dose (no treatment) as the base. Actually I'm realy interested in a dose-response, i.e. consistent down- or up-regulation of a gene by dose. I have turned the Dose factor into a continous variable and taken the log10 of the dose. The model formular stays the same as above. For the factorial desgin I get 112 genes a with significant dose effect, and for the regression model (dose is the regressor and study stays a factor) I get 280 genes at the same significance level. Can I conclude that the regression model is more sensitive? Is there a way to find out where the differences come from and why there are so many differences? There a strong difference using un-transformed or log transformed values for dose, and since the control dose is 0.0mM which value to add for the log transformation. Are any rules of thumb for the value to add? I'd be happy to read about your experience of analyzing dose-response experiments and whether you use a factorial or regression desgin. kind regards, Arne
regression dose • 1.2k views
modified 14.8 years ago by Liaw, Andy360 • written 14.8 years ago by Arne.Muller@aventis.com620
0
14.8 years ago by
Liaw, Andy360
Liaw, Andy360 wrote:
> From: Arne.Muller@aventis.com > > Dear All, > > I'm interested in your experience of analyzing dose-response > experiments. I've an unbalanced design with 3 studies and > within each study 4 to 5 doses of a drug (the doses do not > exactmy the same in each study). > > Usually I run an anova per gene > > anova(lm(Intensity ~ Study + Dose + Study:Dose)) > > where Study and Dose are factors, and for the linear model I > use a treatment contrast with the control dose (no treatment) > as the base. > > Actually I'm realy interested in a dose-response, i.e. > consistent down- or up-regulation of a gene by dose. I have > turned the Dose factor into a continous variable and taken > the log10 of the dose. The model formular stays the same as above. > > For the factorial desgin I get 112 genes a with significant > dose effect, and for the regression model (dose is the > regressor and study stays a factor) I get 280 genes at the > same significance level. > > Can I conclude that the regression model is more sensitive? > Is there a way to find out where the differences come from > and why there are so many differences? > > There a strong difference using un-transformed or log > transformed values for dose, and since the control dose is > 0.0mM which value to add for the log transformation. Are any > rules of thumb for the value to add? > > I'd be happy to read about your experience of analyzing > dose-response experiments and whether you use a factorial or > regression desgin. > > kind regards, > > Arne Arne, I think I've mentioned this before. For questions like these, you're likely to be much better served by consulting a statistician in your organization. There's only so much that can be done over a mailing list, and it might not be sufficient for your full understanding. There are a few possible reasons why the regression' models seem more sensitive than the factorial' models. In the factorial' models, you use up d-1 degrees of freedom for the main effect of dose (assuming there are d doses), plus more for the interaction. In the regression' model, the dose only took up 1 df, thus you have more residual df for estimating the MSE. If the dose effect is close enough to being linear within the range being considered, the `regression' model will be more sensitive. That doesn't mean that the dose effect is linear, but that there's a significant linear component in the dose effects. Regarding transformation of dose, that is also related to the linearity of dose effects. If the effect is linear in log(dose), using log(dose) will be more powerful (i.e., more sensitive), obviously. The treatment of zero dose has been dealt with quite a bit in the statistics literature, but I'm not familiar with that. HTH, Andy