Amigos de la erre.

He creado mi primer bucle con for para entrenar unos modelos con GAM. La
respuesta es quasipoisson porque estoy trabajando con densidades de peces.
Sin embargo, tengo un problema, no se muy bien como aёadir el offset a la
formula siguiente cuando creo el bucle.

GAM.A1 <-gam ((DYO)~s(DMA,k=4)+ s(WOD,k=4)+s(CIN,k=4)+s(DRA,k=4)+s(DBR,k=4)
+offset(log(ARE)),family=quasipoisson(link =
"log"),data=as.data.frame(cbind(DYO,Conjunto)))

Os escribo el codigo del bucle:

for(i in (1:nrow(VarCombS))){
  Formula <- as.formula(paste("DYO",sep=" ~
",paste(as.matrix(VarCombS)[i,1:NVar],sep="",collapse=" + "))) ### Cѓmo
introduzco el offset??
  GAM <- gam(Formula,family=quasipoisson(link =
"log"),data=as.data.frame(cbind(DYO,Conjunto)))
  Summ <- summary(GAM)
  Var.p.values <- Summ$s.table[,"p-value"]
  MaxSignif <-c(MaxSignif,max(Var.p.values))
  Significance <- all(Var.p.values<0.1)  ### Especificar el p values que
queremos dentro del modelo
  GlobalSignificance <- c(GlobalSignificance,Significance)
  RSquare <- c(RSquare,Summ$r.sq)
}

Si teneis alguna sugerencia os estarэa muy agradecido

-- 
*Juan Diego Alcaraz Hernсndez
*
**

*http://juandiegoalcaraz.wordpress.com/*
*
Institut d'Investigaciѓ per a la Gestiѓ Integrada de Zones Costaneres (IGIC)
Universitat Politшcnica de Valшncia
C/ Paranimf, 1
46730 Grau de Gandia (Valшncia)
Tlf: (+34) 963.877.007 (ext: 43040)
<http://juandiegoalcaraz.wordpress.com/>***
ЗДЏ`З.ИИ..><((((К>.ЗДЏ`З.ИИ.ЗДЏ`З.И><((((К>`З.ИИ.ЗДЏ`З.И
 <К))))><`З.ИИ.ЗДЏ`З.И.<К))))><.И. , . .ЗДЏ`З.. <К))))><И.З

	[[alternative HTML version deleted]]

