<html><head><meta http-equiv="Content-Type" content="text/html charset=windows-1252"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Olá pessoal,<div class=""><br class=""></div><div class="">Estou com um problema para fazer a predição de um modelo GAM com uma variável categórica para rasters. </div><div class=""><br class=""></div><div class="">Quando faço a predição do modelo para rasters, configurando a função predict() considerando a variável categórica, ela não funciona.</div><div class=""><br class=""></div><div class="">Segue abaixo um código exemplo ilustrando o meu problema com e sem a variável categórica. </div><div class=""><br class=""></div><div class="">Espero que possam ajudar.</div><div class=""><br class=""></div><div class="">Muito obrigado!</div><div class=""><br class=""></div><div class="">Thiago</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><span style="text-align: -webkit-auto; orphans: 2; widows: 2;" class="">__________________________________</span></div><div class=""><div class=""><span class="Apple-style-span" style="border-collapse: separate; border-spacing: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><span class="Apple-style-span" style="border-collapse: separate; orphans: 2; text-align: -webkit-auto; widows: 2; border-spacing: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><span class="Apple-style-span" style="border-collapse: separate; text-align: -webkit-auto; border-spacing: 0px;"><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Thiago Cesar Lima Silveira<br class="">Biólogo, Me. Biologia Animal<br class="">Doutorando do PPG Zoologia - PUCRS<br class="">e-mail: <a href="mailto:thiagoclsilveira@yahoo.com.br" class="">thiagoclsilveira@yahoo.com.br</a></div><div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Skype: thiagocesarls<br class=""><br class=""><br class=""><br class="">######</div></span></div></span></div></span></div></div><div class=""><div class=""><br class=""></div><div class="">library(mgcv)</div><div class="">library(raster)</div><div class=""><br class=""></div><div class=""># exemplo de layer raster</div><div class=""><br class=""></div><div class="">v1rst<-raster()</div><div class="">values(v1rst) <- 1:ncell(v1rst)</div><div class="">names(v1rst)<-'v1'</div><div class=""><br class=""></div><div class=""># Exemplos de variáveis resposta e preditores</div><div class=""> </div><div class="">y<-c(1,33,500,700, 334,320, 703,303,3030,3002,200,0,100,100,169)</div><div class="">v1<-c(12,33,544,600, 34,30, 03,3390,3030,302,20,108,170,101,2009)</div><div class="">v2<-c(0,0,0,0,0,0,0,0,1,1,1,1,1,1,1)</div><div class="">v3<-c(12,33,544,600, 34,30, 03,3390,3030,302,20,108,170,101,2009)</div><div class="">df<-data.frame(y, v1, v2)</div><div class=""><br class=""></div><div class="">#GAM com variável categórica</div><div class=""><br class=""></div><div class="">gam1<-gam(y~s(v1)+factor(v2), data=df)</div><div class="">summary(gam1)</div><div class=""><br class=""></div><div class="">#GAM model sem factor</div><div class=""><br class=""></div><div class="">gam2<-gam(y~s(v1), data=df)</div><div class="">summary(gam2)</div><div class=""><br class=""></div><div class=""># data.frame com o valor constante</div><div class="">#(da classe ’factor’) para ser usado na função predict( ). </div><div class="">#Esse comando foi retirado do documento do pacote brt (<a href="http://cran.r-project.org/web/packages/dismo/vignettes/brt.pdf" class="">http://cran.r-project.org/web/packages/dismo/vignettes/brt.pdf</a>) página 20.</div><div class=""># Tentei algumas variações sem sucesso.</div><div class=""><br class=""></div><div class="">Method<-factor(var.factor<-factor( '0' ,levels=levels(df$v2))) </div><div class="">df$add<-data.frame(var.factor)</div><div class=""><br class=""></div><div class=""> #Prediction with factor</div><div class=""> p<-predict(v1rst,gam1, type='response', const='add' )</div><div class=""> </div><div class=""> #Mensagem de ERRO - 'Error in `[.data.frame`(blockvals, , f[j]) : undefined columns selected</div><div class=""> </div><div class=""> #Prediction sem o fator</div><div class=""> p<-predict(v1rst,gam2, type='response')</div><div class=""> plot(p) #ok!</div><div class=""><br class=""></div><div class=""><br class=""></div></div><br class=""><br class=""></body></html>