[R-sig-ME] How to test if two gamm-predictions are significantly different?

Chris Howden chris at trickysolutions.com.au
Wed May 30 00:59:45 CEST 2012


Rather than using 95% CI'S for the 2 curves U can try using a 90% CI
for them. If 2 90% CI'S don't overlap it's closer to a 5% t-test than
if 2 95% CI's don't.

But I think your other method for a 95% CI for the difference may be better.

Chris Howden
Founding Partner
Tricky Solutions
Tricky Solutions 4 Tricky Problems
Evidence Based Strategic Development, IP Commercialisation and
Innovation, Data Analysis, Modelling and Training

(mobile) 0410 689 945
(fax / office)
chris at trickysolutions.com.au

Disclaimer: The information in this email and any attachments to it are
confidential and may contain legally privileged information. If you are not
the named or intended recipient, please delete this communication and
contact us immediately. Please note you are not authorised to copy,
use or disclose this communication or any attachments without our
consent. Although this email has been checked by anti-virus software,
there is a risk that email messages may be corrupted or infected by
viruses or other
interferences. No responsibility is accepted for such interference. Unless
expressly stated, the views of the writer are not those of the
company. Tricky Solutions always does our best to provide accurate
forecasts and analyses based on the data supplied, however it is
possible that some important predictors were not included in the data
sent to us. Information provided by us should not be solely relied
upon when making decisions and clients should use their own judgement.

On 30/05/2012, at 0:09, Karel Viaene <karel.viaene at ugent.be> 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