Olá Inajara,<br><br>cole as colunas, ano, mes e dia e depois converta para a classe date.<br>exemplo:<br><br>x <- expand.grid(ano=2000,mes=1, dia=1:31)<br>## gerando dados de chuva<br>x$Precip <- sample(c(0,1), 31, replace=T)<br>
<br>## adicionando coluna "date" ao dataframe<br>x$date <- as.Date(do.call("paste", c(x[,1:3], sep="-")))<br><br>require(seas)<br><br>CCS <- interarrival(x, var="Precip",inv=T)<br>
## inv = T, a data de saída corresponde a data do primeiro dia seco da sequencia de dias secos<br>## a variável de interesse no seu caso é a "dry" <br>CCS<br> date dry wet<br>1 2000-01-02 1 NA<br>2 2000-01-04 2 1<br>
3 2000-01-07 2 1<br>4 2000-01-13 1 4<br>5 2000-01-17 2 3<br>6 2000-01-21 1 2<br>7 2000-01-23 2 1<br>8 2000-01-26 1 1<br><br><br>[ ]'s<br><br><br><div class="gmail_quote">2012/1/5 Inajara Welter Macedo <span dir="ltr"><<a href="mailto:inajarawm@hotmail.com">inajarawm@hotmail.com</a>></span><br>
<blockquote class="gmail_quote" style="margin:0pt 0pt 0pt 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div><div dir="ltr">
<div><span style="font-size:10pt">Ola,</span><br><div dir="ltr"><div><div dir="ltr"><div><br></div><div>Estou trabalhando com dados de precipitacao, de 10 anos, e quero obter resultados de periodos de seca. </div><div>Preciso contar quantos dias consecutivos sem chuva houve, ao todo, e por ano e depois por mes.</div>
<div><br></div><div>Meus dados eram uma matriz de 4 colunas, <span style="font-size:10pt">chamado "</span><b style="font-size:10pt">mod1", </b><span style="font-size:10pt">cada coluna corresponde ao ano, mes, dia, e precipitacao,</span><b style="font-size:10pt"> </b><span style="font-size:10pt">mas transformei para data.frame, para poder utilizar o pacote "seas": </span><span style="font-size:10pt">as.data.frame(mod1)</span></div>
<div><span style="font-size:10pt"><br></span></div><div><span style="font-size:10pt">Tentei usar a funcao "</span><span style="font-size:10pt">interarrival": </span><span style="font-size:10pt"> </span><span style="font-size:10pt">CCS<-interarrival (mod1, var = "Precip", p.cut = 0, inv = TRUE) # CCS= casos consecutivos secos</span></div>
<div><span style="font-size:10pt"><br></span></div><div><span style="font-size:10pt">Mas, ocorre um erro: </span><b><span style="font-size:10pt">Erro em seas.df.check(x, orig, var) : </span><span style="font-size:10pt"> a ‘date’ column must exist in ‘mod1’</span></b></div>
<div><b><span style="font-size:10pt"><br></span></b></div><div>Nao sei como colocar essa coluna de data que pede, porque existem as colunas de ano , mes, e dia. Como consigo converter essas colunas em uma unica coluna de data??</div>
<div><br></div><div><br></div><div>Obrigada, </div><div>Inajara - UFPEL/RS</div><div><br></div><div><span style="font-size:10pt"><br></span></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div>
<div><div><br></div><div><br></div></div> </div></div> </div></div> </div></div>
<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><br clear="all"><br>-- <br>#------------------------------------------------------------------------------#<br>
# Jônatan Dupont Tatsch #<br># Climate and Biosphere Laboratory #<br># Department of Atmospheric Sciences #<br>
# Institute of Astronomy, Geophysics and Atmospheric Sciences #<br># University of São Paulo #<br># Rua do Matão, 1226 #<br>
# Cid. Universitária, São Paulo, SP, Brazil, CEP: 05508-090 #<br># Phone:+55 11 3091-4772, Fax:+55 11 3091-4714 #<br># <a href="http://jonatandupont.weebly.com/index.html" target="_blank">http://jonatandupont.weebly.com/index.html</a> #<br>
#------------------------------------------------------------------------------#<br>