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]]
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]]
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
>