<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:12pt"><div><span>Caros Izidro Pi, Luis Iván e Jônatan,</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span><br></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span>Obrigado pelas suas respostas. Analisei cada uma delas e todas foram úteis, mas o Jônatan foi um passo adiante e antecipou a reposta para a próxima pergunta que eu ia fazer na lista - como combinar dois gráficos mantendo um mesmo eixo x. </span></div><div></div><div> </div><div>Saudações,<br>--<br>Thiago V.
dos Santos<br>PhD student<br>Land and Atmospheric Science<br>University of Minnesota<br>http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm<br>Phone: (612) 323 9898</div> <div class="qtdSeparateBR"><br><br></div><div class="yahoo_quoted" style="display: block;"> <div style="font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-size: 12pt;"> <div style="font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-size: 12pt;"> <div dir="ltr"> <font size="2" face="Arial"> On Thursday, September 11, 2014 8:13 AM, Jônatan <jdtatsch@gmail.com> wrote:<br> </font> </div> <br><br> <div class="y_msg_container"><div id="yiv6293530324"><div><div dir="ltr"><div>Segue algumas alternativas:</div><div><br clear="none"></div><div>library(repmis) # reads text data directly from dropbox - no need to download any file</div><div># read data</div><div>url
<- '<a rel="nofollow" shape="rect" target="_blank" href="https://dl.dropboxusercontent.com/u/27700634/precip.txt">https://dl.dropboxusercontent.com/u/27700634/precip.txt</a>'</div><div>tmp <- repmis::source_data(url, sep = '', header = TRUE)</div><div><br clear="none"></div><div># convert julian day to date</div><div>date <- as.Date(tmp$julian, origin='2011-12-31')</div><div>data <- cbind(date, tmp)</div><div>head(data)</div><div># define vector with depth of soil layers</div><div>depths <- c(0.05,0.10,0.20,0.30,</div><div> 0.40,0.60,0.80,1.00,</div><div> 1.50,2.00,2.50)</div><div><br clear="none"></div><div># teste com lattice</div><div>cols <- ncol(data):4</div><div>levelplot(as.matrix(data[, cols]), </div><div> aspect = "fill", xlab = "date", ylab = "soil
Moist.",</div><div> scales = list(x = list(at = which(date %in% pretty(date)), </div><div> labels = pretty(date)),</div><div> y = list(at = length(cols):1, labels = depths)),</div><div> contour = T)</div><div><br clear="none"></div><div># Plot soil moisture profile</div><div><br clear="none"></div><div># como são profundidades (abaixo da superfície) use a convenção de z negativo </div><div># acrescentei a chuva, precisa ajustar ylabels, mtext() resolve</div><div><br clear="none"></div><div>x11()</div><div>op <- par(mar = c(3,4,12,3))</div><div>plot <- filled.contour(date, sort(-depths), as.matrix(data[, cols]), ylab =
expression(theta(m^3~~m^-3)))</div><div>par(new=T,mar=c(23,4,1,8))</div><div> with(data,</div><div> plot(date, precip, type = "h", xaxt = "n", xpd = T, xaxs="i", ylab = "Precip (mm)", xlab = ""))</div><div>par(op)</div><div><br clear="none"></div><div><br clear="none"></div><div><br clear="none"></div></div><div class="yiv6293530324gmail_extra"><br clear="none"><div class="yiv6293530324gmail_quote">2014-09-10 16:15 GMT-03:00 Thiago V. dos Santos <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:thi_veloso@yahoo.com.br" target="_blank" href="mailto:thi_veloso@yahoo.com.br">thi_veloso@yahoo.com.br</a>></span>:<br clear="none"><blockquote class="yiv6293530324gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;"><div class="yiv6293530324yqt3440217636" id="yiv6293530324yqt20058"><div><div style="color: rgb(0, 0, 0); font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial,
'Lucida Grande', sans-serif; font-size: 12pt; background-color: rgb(255, 255, 255);"><div><span>Olá pessoal,</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span><br clear="none"></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span>Estou tentando fazer um gráfico de contorno no R mas ainda não encontrei a função certa para isso. </span><span style="background-color:transparent;">Basicamente, eu tenho uma simulação de umidade em 11 camadas do solo
para cada dia do ano, e gostaria de mostrar isso na forma de um perfil de solo.</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span style="background-color:transparent;"><br clear="none"></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span style="background-color:transparent;">Por favor leia os comentários ao reproduzir os meus dados com o código abaixo:</span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span style="background-color:transparent;"><br clear="none"></span></div><div
style="background-color:transparent;">library(repmis) # reads text data directly from dropbox - no need to download any file</div><div style="background-color:transparent;"><br clear="none"></div><div style="background-color:transparent;"># read data</div><div style="background-color:transparent;">url <- '<a rel="nofollow" shape="rect" target="_blank" href="https://dl.dropboxusercontent.com/u/27700634/precip.txt">https://dl.dropboxusercontent.com/u/27700634/precip.txt</a>'</div><div style="background-color:transparent;">tmp <- repmis::source_data(url, sep = '', header = TRUE)</div><div style="background-color:transparent;"><br clear="none"></div><div style="background-color:transparent;"># convert julian day to date</div><div style="background-color:transparent;">date <- as.Date(tmp$julian, origin='2011-12-31')</div><div style="background-color:transparent;">data <-
cbind(date, tmp)</div><div style="background-color:transparent;">head(data)</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><br clear="none"></div><div style="background-color:transparent;"># now, convert soil layers to matrix and transpose it</div><div style="background-color:transparent;">mat <- t(as.matrix(data[, 4:14]))</div><div style="background-color:transparent;"># essa é a matriz que eu quero plotar. Reparem que ela já está organizada como um "perfil",</div><div style="background-color:transparent;"><span style="background-color:transparent;"></span></div><div style="background-color:transparent;"># com as linhas representando as camadas de solo e as colunas representando os dias do ano</div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue',
Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span style="background-color:transparent;"><br clear="none"></span></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span style="background-color:transparent;">A minha primeira tentativa foi usar a função </span><span style="font-size:12pt;">filled.contour</span><span style="font-size:12pt;"> do pacote graphics. Eu defino um vetor com o nome das camadas do solo e então eu tento plotar, mas a função diz que a matriz está com as dimensoes invertidas:</span><br clear="none"></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span
style="font-size:12pt;"><br clear="none"></span></div><div style="background-color:transparent;"># define vector with depth of soil layers</div><div style="background-color:transparent;">depths <- c(0.05,0.10,0.20,0.30,</div><div style="background-color:transparent;"> 0.40,0.60,0.80,1.00,</div><div style="background-color:transparent;"> 1.50,2.00,2.50)</div><div style="background-color:transparent;"><br clear="none"></div><div style="background-color:transparent;"># Plot soil moisture profile</div><div style="background-color:transparent;"><span></span></div><div style="background-color:transparent;">plot <-
filled.contour(data$julian, depths, mat)</div><div><div><br clear="none"></div><div>#Error in .filled.contour(x, y, z, levels, col) : dimension mismatch</div><div><br clear="none"></div></div><div>Obviamente eu consigo transpor a matrix para forçar o gráfico, mas a figura final não é o que eu preciso - o perfil fica de "cabeça para baixo".</div><div><div style="margin-top:0.1em;margin-bottom:0.1em;background-color:transparent;">plot <- filled.contour(data$julian, depths, t(mat))</div><div style="margin-top:0.1em;margin-bottom:0.1em;background-color:transparent;"><br clear="none"></div><div style="margin-top: 0.1em; margin-bottom: 0.1em; color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;">Como essa função só aceita os eixos em ordem crescente, eu não consigo inverter o eixo y para mostrar a primeira camada do solo
no topo do gráfico.</div><div></div></div><div><br clear="none"></div><div style="color: rgb(0, 0, 0); font-size: 16px; font-family: HelveticaNeue, 'Helvetica Neue', Helvetica, Arial, 'Lucida Grande', sans-serif; font-style: normal; background-color: transparent;"><span>Poranto, eu aprecio qualquer sugestão de funções e/ou pacotes para construir o gráfico que eu preciso. </span></div><div></div><div> </div><div>Saudações,<br clear="none">--<br clear="none">Thiago V. dos Santos<br clear="none">PhD student<br clear="none">Land and Atmospheric Science<br clear="none">University of Minnesota<br clear="none"><a rel="nofollow" shape="rect" target="_blank" href="http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm">http://www.laas.umn.edu/CurrentStudents/MeettheStudents/ThiagodosSantos/index.htm</a><br clear="none">Phone: <a rel="nofollow" shape="rect" href="">(612) 323 9898</a></div></div></div></div><br
clear="none">_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none">
<a rel="nofollow" shape="rect" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">
Leia o guia de postagem (<a rel="nofollow" shape="rect" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br clear="none"></blockquote></div><br clear="none"><br clear="all"><div><br clear="none"></div>-- <br clear="none"><div dir="ltr"><div><font style="font-family: arial;">###############################################################</font><br clear="none"><span style="font-family: arial; font-size: small;">## </span><span style="font-size: 11px; text-align: center; color: rgb(150, 150, 150); font-family: Verdana, Arial, Helvetica, sans-serif; line-height: 13px; background-color: rgb(239, 239, 239);">Jônatan Dupont Tatsch</span></div><div><span style="font-family: arial; font-size: small;">## </span><span style="color: rgb(150, 150, 150); font-family: Verdana, Arial, Helvetica, sans-serif; font-size: 10.9090909957886px; line-height: 12.9971590042114px;
text-align: center; background-color: rgb(239, 239, 239);">Professor do Departamento de Física</span><font face="arial" style="font-size:small;"><br clear="none"></font></div><div><span style="font-size: small; font-family: arial;">## </span><span style="color: rgb(150, 150, 150); font-family: Verdana, Arial, Helvetica, sans-serif; font-size: 10.9090909957886px; line-height: 12.9971590042114px; text-align: center; background-color: rgb(239, 239, 239);">Centro de Ciências Exatas e Naturais (CCNE)</span><br clear="none"></div><div><span style="font-size: small; font-family: arial;">## </span><span style="color: rgb(150, 150, 150); font-family: Verdana, Arial, Helvetica, sans-serif; font-size: 10.9090909957886px; line-height: 12.9971590042114px; text-align: center; background-color: rgb(239, 239, 239);">Universidade Federal de Santa Maria</span></div><div><span style="font-family: arial; font-size: small;">## </span><span
style="font-size: 11px; text-align: center; color: rgb(150, 150, 150); font-family: Verdana, Arial, Helvetica, sans-serif; line-height: 13px; background-color: rgb(239, 239, 239);">Faixa de Camobi, Prédio 13 - Campus UFSM - </span><span style="font-size: 11px; text-align: center; color: rgb(150, 150, 150); font-family: Verdana, Arial, Helvetica, sans-serif; line-height: 13px; background-color: rgb(239, 239, 239);">Santa Maria, RS, Brasil - 97105-900</span><br clear="none"></div><div style="font-family: arial; font-size: small;">## <span style="font-size: 11px; text-align: center; color: rgb(150, 150, 150); font-family: Verdana, Arial, Helvetica, sans-serif; line-height: 13px; background-color: rgb(239, 239, 239);">Telefone: +55(55)33012083</span><br clear="none"></div><div style="font-family: arial; font-size: small;">## <span style="font-size: 11px; text-align: center; color: rgb(150, 150, 150); font-family: Verdana, Arial, Helvetica,
sans-serif; line-height: 13px; background-color: rgb(239, 239, 239);"><a rel="nofollow" shape="rect" target="_blank" href="http://www.ufsm.br/meteorologia">www.ufsm.br/meteorologia</a></span></div><div style="font-family: arial; font-size: small;">###############################################################</div></div>
</div></div></div><br><br></div> </div> </div> </div> </div></body></html>