abaixo...<br><br><br clear="all">---<br>Fernando Mayer<br>Universidade Federal de Santa Catarina - UFSC<br>Departamento de Ecologia e Zoologia - ECZ/CCB<br>URL: <a href="http://sites.google.com/site/fernandomayer" target="_blank">http://sites.google.com/site/fernandomayer</a><br>
e-mail: fernandomayer [@] <a href="http://gmail.com" target="_blank">gmail.com</a><br>
<br><br><div class="gmail_quote">2012/3/16 Thiago Veloso <span dir="ltr"><<a href="mailto:thi_veloso@yahoo.com.br">thi_veloso@yahoo.com.br</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font:inherit" valign="top"><div><font face="arial"> Pessoal,</font></div><div><font face="arial"><br></font></div><div><font face="arial"> Possuo três mapas (em formato raster) de desmatamento na Amazônia em diferentes anos. Essas taxas de desmatamento variam de 0 a 1. A minha intenção é mostrar a variação temporal do desmatamento através de histogramas que mostram a contagem de células em diferentes classes (p. ex., de 0 a 0.3, de 0.3 a 0.6 e de 0.6 a 1) em cada ano. Aqui está uma ideia conceitual do que preciso fazer (apenas ilustrativa, criada com dados hipotéticos): <a href="http://img24.imageshack.us/img24/6416/examplezxt.png" target="_blank">http://img24.imageshack.us/img24/6416/examplezxt.png</a> e aqui disponibilizo um arquivo com amostra do arquivo de dados <a href="http://db.tt/XXDtRee1" target="_blank">http://db.tt/XXDtRee1</a>. Trata-se um netcdf com três níveis de tempo.</font></div>
<div><font face="arial"><br></font></div><div><font face="arial"> Na lista de SIG do R eu já recebi algumas dicas para fazer isso. Entretanto, estou recebendo um erro ao tentar criar fatores nos dados. Será que alguém aqui consegue me dar uma auxílio? Segue o código que estou utilizando e a mensagem de erro:</font></div>
<div><font face="arial"><br></font></div><div><font face="arial">--------------------------------------------</font></div><div><font face="arial">> #Load required packages</font></div><div><font face="arial">> library (raster)</font></div>
<div><font face="arial">Loading required package: sp</font></div><div><font face="arial">raster version 1.9-70 (27-February-2012)</font></div><div><font face="arial">> library (ggplot2)</font></div><div><font face="arial">> library (reshape2)</font></div>
<div><font face="arial">> </font></div><div><font face="arial">> #Load
data</font></div><div><font face="arial">> deforest <- brick("/home/thiago/Dropbox/<a href="http://defor-3-years.nc" target="_blank">defor-3-years.nc</a>")</font></div><div><font face="arial">Loading required package: ncdf</font></div>
<div><font face="arial">> </font></div><div><font face="arial">> #Probe layer names and plot raster data</font></div><div><font face="arial">> layerNames(deforest)</font></div><div><font face="arial">[1] "X0.1" "X0.2" "X0.3"</font></div>
<div><font face="arial">> plot(deforest)</font></div><div><font face="arial">> </font></div><div><font face="arial">> #Convert spatial object to data frame (omitting NAs)</font></div><div><font face="arial">> meltR <- melt(as.array(deforest),na.rm=TRUE)</font></div>
<div><font face="arial">> </font></div><div><font face="arial">> #Define classes used to count grids</font></div><div><font face="arial">> meltR$class <- cut(meltR$value, c(0,0.3,0.6,1))</font></div><div><font face="arial">> #Define factors to plot individual years</font></div>
<div><font face="arial">> meltR$X3 <- factor(meltR$X3, labels=c('2000','2001','2002'))</font></div><div><font face="arial">Error in factor(meltR$X3, labels = c("2000", "2001", "2002")) : </font></div>
<div><font face="arial"> invalid labels; length 3 should be 1 or 0</font></div></td></tr></tbody></table></blockquote><div><br><br>No objeto meltR não existe nenhum X3 como componente, por isso não há como criá-lo. Qual era o objetivo aqui?<br>
<br> </div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><table border="0" cellpadding="0" cellspacing="0"><tbody><tr><td style="font:inherit" valign="top">
<div><font face="arial">--------------------------------------------</font></div><div><font face="arial"><br></font></div><div><font face="arial"> O código continuaria assim:</font></div><div><font face="arial"><br></font></div>
<div><font face="arial">--------------------------------------------</font></div><div><font face="arial">#Plot data</font></div><div><font face="arial">p =
ggplot(meltR,aes(x=class))</font></div><div><font face="arial">p + geom_histogram(aes(fill=class))+facet_grid(. ~ X3)</font></div><div><font face="arial">--------------------------------------------</font></div><div><font face="arial"><br>
</font></div><div><font face="arial"> Caso seja útil, segue o sumário do sessionInfo():</font></div><div><font face="arial"><br></font></div><div><font face="arial">--------------------------------------------</font></div>
<div><font face="arial">> sessionInfo()</font></div><div><font face="arial">R version 2.14.1 (2011-12-22)</font></div><div><font face="arial">Platform: x86_64-pc-linux-gnu (64-bit)</font></div><div><font face="arial"><br>
</font></div><div><font face="arial">locale:</font></div><div><font face="arial"> [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
</font></div><div><font face="arial"> [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8 </font></div><div><font face="arial"> [5] LC_MONETARY=en_US.UTF-8 LC_MESSAGES=en_US.UTF-8 </font></div><div><font face="arial"> [7] LC_PAPER=C LC_NAME=C </font></div>
<div><font face="arial"> [9] LC_ADDRESS=C LC_TELEPHONE=C </font></div><div><font face="arial">[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C </font></div><div><font face="arial"><br></font></div>
<div><font face="arial">attached base packages:</font></div><div><font face="arial">[1] stats
graphics grDevices utils datasets methods base </font></div><div><font face="arial"><br></font></div><div><font face="arial">other attached packages:</font></div><div><font face="arial">[1] ncdf_1.6.6 reshape2_1.2.1 ggplot2_0.9.0 raster_1.9-70 sp_0.9-97 </font></div>
<div><font face="arial"><br></font></div><div><font face="arial">loaded via a namespace (and not attached):</font></div><div><font face="arial"> [1] colorspace_1.1-1 dichromat_1.2-4 digest_0.5.2 grid_2.14.1 </font></div>
<div><font face="arial"> [5] lattice_0.20-6 MASS_7.3-17 memoise_0.1 munsell_0.3 </font></div><div><font face="arial"> [9]
plyr_1.7.1 proto_0.3-9.2 RColorBrewer_1.0-5 scales_0.2.0 </font></div><div><font face="arial">[13] stringr_0.6 </font></div><div><font face="arial">--------------------------------------------</font></div>
<div><font face="arial"><br></font></div><div><font face="arial"> Um verdadeiro desafio - estou desde 4a feira quebrando a cabeça com isso. Agradeço desde já qualquer contribuição!</font></div><span class="HOEnZb"><font color="#888888"><div>
<font face="arial"> Thiago.</font></div></font></span></td></tr></tbody></table><br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br></blockquote></div><br>