Re: [R] gam confidence interval (package mgcv)

2011-06-28 Thread Simon Wood
not sure if I'm missing something here, but since you are using a log link, isn't the ratio you are looking for given by the `treatmentB' parameter in the summary (independent of X) > summary(gfit) [snip] Parametric coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 1.

Re: [R] gam confidence interval (package mgcv)

2011-06-27 Thread Remko Duursma
But that just gives me the prediction of Y for treatment A or B, not the ratio. As I stated: # I am interested in the relationship: # Y(treatment =="B") / Y(treatment=="A") as a function of X, with a confidence interval! I can get the SE for either of them using predict.gam without a problem, b

Re: [R] gam confidence interval (package mgcv)

2011-06-27 Thread David Winsemius
On Jun 27, 2011, at 10:45 PM, Remko Duursma wrote: Dear R-helpers, I am trying to construct a confidence interval on a prediction of a gam fit. I have the Wood (2006) book, and section 5.2.7 seems relevant but I am not able to apply that to this, different, problem. Any help is appreciated!

[R] gam confidence interval (package mgcv)

2011-06-27 Thread Remko Duursma
Dear R-helpers, I am trying to construct a confidence interval on a prediction of a gam fit. I have the Wood (2006) book, and section 5.2.7 seems relevant but I am not able to apply that to this, different, problem. Any help is appreciated! Basically I have a function Y = f(X) for two different