
OI EDER, SOU EU ANANDA. EU ACHO QUE VC TERIA QUE RODAR ANALISE CIRCULAR. TENHO ALGUMAS INFORMAÇÕES SOBRE ESSE ASSUNTO. ANANDA AGUIAR Em 4 de março de 2012 11:53, Eder David Borges da Silva <eder@leg.ufpr.br>escreveu:
Bom dia pessoal, Tenho um conjunto de dados, de fenologia, (fases da planta: florescimento, frutificação, etc), e periodicamente é dado um nota de 0 a 5 para percentual da fase, por exemplo entre 40 a 60 % de florescimento recebe nota 3, e assim por diante dependendo do percentual temos uma nota. Oque me interessa saber é ter primeiramente um descritivo, por exemplo de valor médio da fase, e IC para a fase, do tipo, com 95% de confiança o florescimento com nota 1 vai de 3 de janeiro a 10 de fevereiro. Ai surge minha duvidas como fazer estes cálculos com data no R, pois não é simplesmente pedir um summary por exemplo pois oque me interessa e o período. Um exemplo do conjunto de dados é simulados abaixo:
ny <- 5 #numero de anos n <- 365*ny #dias inicio <- as.Date(c('01/01/2000'),'%d/%m/%Y') periodo <- seq(inicio,inicio+(n-1),by=1) y <- as.vector(sapply(1:ny,function(x)sort(ceiling(runif(n/ny,0,4))))) y1 <- as.vector(sapply(1:ny,function(x)sort(ceiling(runif(n/ny,0,4))))) dados <- data.frame(rep=rep(1:2,each=n),periodo,resp=c(y,y1)) head(dados) lattice::xyplot(resp~periodo,groups=rep,dados)
Gostaria de um descritivo por nota, alguma sugestão a mais será muito bem vinda. Obrigado. _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
-- Att, Fernando Rodrigo Bortolozo