[R-sig-ME] How to test if two gamm-predictions are significantly different?
Kenneth Frost
kfrost at wisc.edu
Wed May 30 02:54:47 CEST 2012
Hi, Karel-
Could you provide a longer sketch of your situation? It is not clear to me what your control treatment is and how it relates to concentration or time.
Ken
On 05/29/12, Karel Viaene wrote:
> Dear R community,
>
> A quick sketch of my situation:
>
> I have two continuous explanatory variables ("concentration" and "time")
> and a continuous response variable, "biomass".
>
> I've fitted a gamm model to these data using the package mgcv and want
> to predict at what concentration the biomass is significantly different
> from the control treatment (i.e. a concentration of 0) for a given point
> in time.
>
> I've done this by predicting the biomass for a series of 1000
> concentrations at a given point in time (using "predict"), constructing
> 95% confidence intervals for these predictions by adding and subtracting
> 1.96*SE and then selecting the lowest concentration where the two CI
> show no overlap.
>
> However I've realized that this technique is not adequate because two
> points can also be significantly different at the 5% significance level
> when the 95% CI do overlap and I want to calculate the lowest possible
> concentration. I've read some literature about this and am considering
> the following method:
>
> * Calculate the difference between the control (C0) and a predicted
> point (e.g. C1), thus C0-C1.
> * Construct a 95% CI for this difference by adding and subtracting
> 1.96*sqrt(SE0^2 + SE1^2).
> * Do this for all predictions.
> * Select the lowest concentration where the 95% CI does not include 0.
>
> Could you give me some feedback about this as I'm unsure if this method
> can be used for gamms. Any comments or suggestions are much appreciated.
>
> Many thanks in advance & kind regards
>
> Karel
>
> --
> Karel Viaene
> Ghent University
> Laboratory of Environmental Toxicology and Aquatic Ecology
> Plateaustraat 22
> 9000 Ghent, Belgium
> tel: +32 (0) 9 264 3779
>
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-mixed-models at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-sig-mixed-models
More information about the R-sig-mixed-models
mailing list