search for: fstation

...n models of higher orders (ARMA)). If I rename weeks as week 1-week 30 to get unique identifiers, I loose the seasonal and year-specific effect in the final model. M2c.gamm <- gamm(het ~ s(LN.DIN, k=3) + s(LN.totn, k=3) + s(Ncell, k=3) + s(LN.biom) + s(temp, k=3) + s(week, k=3) + fstation + fyear, method = "ML", weights = varIdent(form=~1 | fstation), data = data1, correlation = corCompSymm(form = ~ week|year)) #seems not to work in a gamm() Thank you for your time! Anna Zakrisson Braeunlich PhD Student Department of Systems Ecology Stockholm Universit...

.... I plot "het" against "temp" conditional on "year" and "station". library(lattice) MyLines <- function(xi, yi, ...){ #Draws line in the panels while avoiding spaghetti-plots I <- order(xi) panel.lines(xi[I], yi[I], col = 1)} xyplot(het ~ temp | fstation, data = data1, groups = fyear, xlab = "temp", ylab = "heterocyst freq", panel = panel.superpose, #ensures that lines from all years (defined by "groups") from the same stations are superimposed in he same panel panel.groups = MyLines) #specifies which task should...

...pregunta hacerca de la funcion offset era si para este caso en el cual > la variable independiente es una fraccion es valido. En mixed > effects models and extension in R (pag. 240) hay un caso similar en el que > se utilizan el numero de ind por volumen ( > L0<-glm(ind~offset(LVol)+fStation,family=poisson,data=Lice). Pero no estoy > seguro de que sea lo mismo dado que lo mio es una freq? > > Saludos, > Matias > > > > ------------------------------ > From: jorgeivanvelez@gmail.com > Date: Tue, 11 Jun 2013 18:40:03 +1000 > Subject: Re: [R-es] calculo de...

...la funcion offset era si para este caso en el cual > > la variable independiente es una fraccion es valido. En mixed > > effects models and extension in R (pag. 240) hay un caso similar en el que > > se utilizan el numero de ind por volumen ( > > L0<-glm(ind~offset(LVol)+fStation,family=poisson,data=Lice). Pero no estoy > > seguro de que sea lo mismo dado que lo mio es una freq? > > > > Saludos, > > Matias > > > > > > > > ------------------------------ > > From: jorgeivanvelez en gmail.com > > Date: Tue, 11 Jun 201...

...pregunta hacerca de la funcion offset era si para este caso en el cual > la variable independiente es una fraccion es valido. En mixed > effects models and extension in R (pag. 240) hay un caso similar en el que > se utilizan el numero de ind por volumen ( > L0<-glm(ind~offset(LVol)+fStation,family=poisson,data=Lice). Pero no estoy > seguro de que sea lo mismo dado que lo mio es una freq? > > Saludos, > Matias > > > > ------------------------------ > From: jorgeivanvelez en gmail.com > Date: Tue, 11 Jun 2013 18:40:03 +1000 > Subject: Re: [R-es] calculo...

...la funcion offset era si para este caso en el cual > > la variable independiente es una fraccion es valido. En mixed > > effects models and extension in R (pag. 240) hay un caso similar en el que > > se utilizan el numero de ind por volumen ( > > L0<-glm(ind~offset(LVol)+fStation,family=poisson,data=Lice). Pero no estoy > > seguro de que sea lo mismo dado que lo mio es una freq? > > > > Saludos, > > Matias > > > > > > > > ------------------------------ > > From: jorgeivanvelez@gmail.com > > Date: Tue, 11 Jun 2013 1...

...mal y link= log. Mi pregunta hacerca de la funcion offset era si para este caso en el cual la variable independiente es una fraccion es valido. En mixed effects models and extension in R (pag. 240) hay un caso similar en el que se utilizan el numero de ind por volumen ( L0<-glm(ind~offset(LVol)+fStation,family=poisson,data=Lice). Pero no estoy seguro de que sea lo mismo dado que lo mio es una freq? Saludos, Matias ------------------------------ From: jorgeivanvelez@gmail.com Date: Tue, 11 Jun 2013 18:40:03 +1000 Subject: Re: [R-es] calculo de poder estadistico en glm To: matutetote@hotmail.com...

...la funcion offset era si para este caso en el cual > > la variable independiente es una fraccion es valido. En mixed > > effects models and extension in R (pag. 240) hay un caso similar en el que > > se utilizan el numero de ind por volumen ( > > L0<-glm(ind~offset(LVol)+fStation,family=poisson,data=Lice). Pero no estoy > > seguro de que sea lo mismo dado que lo mio es una freq? > > > > Saludos, > > Matias > > > > > > > > ------------------------------ > > From: jorgeivanvelez@gmail.com > > Date: Tue, 11 Jun 2013 1...