On Jul 1, 2014, at 11:33 PM, agostinodicia...@tiscali.it wrote: > Dear David, many thanks for your reply. > The aim was to evaluate the relationship between daily mortality (variable > "mortality") and mean daily values of PM10 levels (variable "PM10") in a > specific geographic area, considering the effect of meteorological > confounders (temperature, umidity). > I used the following model (mgcv package): > gam(mortality ~ (PM10) + (Tmax) + (umidity), data = data, family = > quasipoisson) > > the question is: how can I obtain a plot of log-relative risk of daily > mortality vs. PM10 levels? > regards
Still rather unclear what the problem is. I see no error message in your question. Are you asking how to assign a model to a name and then call the plot method? mdl <- gam(mortality ~ (PM10) + (Tmax) + (umidity), data = data, family = quasipoisson) plot(mdl) -- david. > > Agostino > Il 02.07.2014 01:50 David Winsemius ha scritto: > >> On Jul 1, 2014, at 12:02 AM, adc wrote: >>> I performed the following GAM by the MGCV package: >> I think it's actually spelled in all lower case. >>> gam(mortality ~ (PM10) + (Tmax) + (umidity), data = data, family = >>> quasipoisson) How can I obtain a plot of Log-relative risk of mortality vs. >>> PM10 ? thanks >> Shouldn't we need to know more details about the experimental setup to >> answer that question? And what sort of comparisons you are requesting? And >> about what parts of ?mgcv::plot.gam you need further explanations to answer >> the question? >>> agostino >> snipped >>> -- >> snipped >>> Sent from the R help mailing list archive at Nabble.com. [[alternative HTML >>> version deleted]] >> snipped >>> PLEASE do read the posting guide >>> http://www.R-project.org/posting-guide.html and provide commented, minimal, >>> self-contained, reproducible code. >> David Winsemius >> Alameda, CA, USA >> >> ______________________________________________ >> >> R-help@r-project.org >> mailing list >> >> https://stat.ethz.ch/mailman/listinfo/r-help >> >> PLEASE do read the posting guide >> http://www.R-project.org/posting-guide.html >> >> and provide commented, minimal, self-contained, reproducible code. >> > > > > Scopri istella, il nuovo motore per il web italiano. > Istella garantisce risultati di qualità e la possibilità di condividere, in > modo semplice e veloce, documenti, immagini, audio e video. > Usa istella, vai su http://www.istella.it?wtk=amc138614816829636 > David Winsemius Alameda, CA, USA ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
