[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