How exactly do you want to define your weights? varIdent returns a varFunc object, not a vector with the same length as the data variables, as required by the gam function.
Cheers Joris On Fri, Jun 18, 2010 at 7:35 PM, niallllll <n.mcgin...@nuigalway.ie> wrote: > > Hello, > > As I am relatively new to the R environment this question may be either > a) Really simple to answer > b) Or I am overlooking something relatively simple. > > I am trying to add a VarIdent structure to my gam model which is fitting > smoothing functions to the time variables year and month for a particular > species. When I try to add the varIdent weights to variable Month I get this > error returned. > > Error in model.frame.default(formula = Chel ~ 1 + M + YM + M + YM, weights = > Vf1, : > variable lengths differ (found for '(weights)') > > Here is the simple code I am using. > library(mgcv); library(lattice) > library(nlme) > Vf1=varIdent(form=~M) > Chel2=gam(Chel~s(M, fx=F, bs="tp")+ s(YM, fx=F, bs="tp")+s(M,YM, fx=F, > bs="tp"), family=poisson(link="log"), weights=Vf1) # second GAM with > interaction > > M=Month > YM=Year > > Regards, > Niall > -- > View this message in context: > http://r.789695.n4.nabble.com/varIdent-error-using-gam-function-in-mgcv-tp2260556p2260556.html > Sent from the R help mailing list archive at Nabble.com. > > ______________________________________________ > 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. > -- Joris Meys Statistical consultant Ghent University Faculty of Bioscience Engineering Department of Applied mathematics, biometrics and process control tel : +32 9 264 59 87 joris.m...@ugent.be ------------------------------- Disclaimer : http://helpdesk.ugent.be/e-maildisclaimer.php ______________________________________________ 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.
