[R-es] Fwd: offset en bucle
Carlos J. Gil Bellosta
cgb en datanalytics.com
Jue Jun 20 13:13:40 CEST 2013
Hola, ¿qué tal?
Pues entiendo que tendrás que abundar en la expresión
as.formula(paste("DYO",sep=" ~
",paste(as.matrix(VarCombS)[i,1:NVar],sep="",collapse=" + ")))
para componer usando paste() una fórmula que al final incluya la
coletilla " + offset( log( ARE )) ".
Si no entiendo mal, es todo lo que necesitas.
Un saludo,
Carlos J. Gil Bellosta
http://www.datanalytics.com
El día 18 de junio de 2013 22:19, Juan Diego Alcaraz-Hernández
<jdalcaraz en gmail.com> escribió:
> 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/>***
> ·´¯`·.¸¸..><((((º>.·´¯`·.¸¸.·´¯`·.¸><((((º>`·.¸¸.·´¯`·.¸
> <º))))><`·.¸¸.·´¯`·.¸.<º))))><.¸. , . .·´¯`·.. <º))))><¸.·
>
>
>
>
> --
> *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]]
>
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
Más información sobre la lista de distribución R-help-es