<div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Cassiano, olá novamente!<br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Atualizei o código, pois percebi que havia um erro na atribuição dos labels (meses) que estavam "counter-clock".<br><br><img alt="Imagem inline 1" src="cid:ii_15535a6eb5645c56" height="388" width="569"><br></div><div class="gmail_default" style="font-family:verdana,sans-serif"><br><span style="font-family:monospace,monospace">### <code r><br>library(circular)<br>nascente <- read.table(text=<br>"AMBIENTE    mes    angle    DISTANCIA    ABUND    RIQUEZA<br>NASC        1    0    0    2    2<br>NASC        2    30    0    3    1<br>NASC        3    60    0    4    3<br>NASC        4    90    0    8    5<br>NASC        5    120    0    0    0<br>NASC        6    150    0    0    0<br>NASC        7    180    0    31    1<br>NASC        8    210    0    10    3<br>NASC        9    240    0    18    3<br>NASC        10    270    0    7    5<br>NASC        11    300    0    2    2<br>NASC        12    330    0    1    1", head=T)<br><br>nascente$meses <- c("Mai","Jun","Jul","Ago","Set","Out","Nov","Dez","Jan","Fev","Mar","Abr")<br>nascente<br><br>#convertendo os dados para circular, usando radianos<br>nasc.c0 <- circular(rep(rad(nascente$angle),nascente$ABUND), units=c("radians"), rota="clock")<br>nasc.c0<br>plot(nasc.c0, stack=T)<br>deg(nasc.c0)<br><br># Angulo médio + Desvio padrão angular<br>media <- mean(nasc.c0)<br>dp    <- sd(nasc.c0)<br>n     <- length(nasc.c0)<br>err   <- dp/sqrt(n)<br>len   <- rho.circular(nasc.c0)<br>cbind(media, dp, n, err, len)<br><br>sd(nasc.c0); sqrt(-2*log(rho.circular(nasc.c0)))<br><br># Gráfico<br># par(mfrow=c(1,3))<br># par(mar=c(0,0,2,0))<br>plot(nasc.c0,axes=F, main=NA, stack=F, col=3)<br>axis.circular(at=circular(rad(nascente$angle), rota="clock"), labels=nascente$meses, units=c("radians"), col=4)<br>rose.diag(nasc.c0, add=T, axes=F, col=8, bor=8, lwd=3)<br># arrows.circular(median(nasc.c0),y=rho.circular(nasc.c0),lwd=3,lty=1,angle=30,len=0.1,col="green") #aqui eu inseri a mediana<br>arrows.circular(mean(nasc.c0),y=len,   lwd=2, lty=1, angle=90, len=0.1, col=2) <br>arrows.circular(media+err,     y=1,    lwd=1, lty=3, angle=0, len=0.1, col=2) <br>arrows.circular(media-err,     y=1,    lwd=1, lty=3, angle=0, len=0.1, col=2)<br><br>nascente<br>### </code></span><br><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><span><br>================================================<br>Éder Comunello<br>Researcher at Brazilian Agricultural Research Corporation (Embrapa)<br>DSc in Agricultural Systems Engineering (USP/Esalq)<br>MSc in Environ. Sciences (UEM), <span>Agronomist (UEM)</span><br>---<br>Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>|<br>================================================<br>GEO, -22.2752, -54.8182, 408m<br>UTC-04:00 / DST: UTC-03:00</span><div><div><div><br></div><div><br></div></div><div style="font-size:small"><br></div></div></div></div></div></div></div></div></div></div></div></div></div>
<br><div class="gmail_quote">Em 9 de junho de 2016 08:41, Éder Comunello <span dir="ltr"><<a href="mailto:comunello.eder@gmail.com" target="_blank">comunello.eder@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><div class="gmail_default" style="font-family:verdana,sans-serif">Cassiano, bom dia!<br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif">Imagino que o que você esteja tentando adicionar ao gráfico é uma linha da média + erro tal qual essa que aparece no gráfico abaixo:<br><img src="http://www.kovcomp.co.uk/oriana/oricirc.gif" alt="" style="margin-top:43px;margin-right:0px" height="289" width="326"><br><br></div><div class="gmail_default" style="font-family:verdana,sans-serif;display:inline">Infelizmente, não encontrei como reproduzir esse efeito no R. O melhor que consegui foi marcar os limites com novas linhas.​</div><div class="gmail_extra"><img alt="Imagem inline 1" src="cid:ii_155352a72e73312d" height="388" width="569"><br><br><br><span style="font-family:monospace,monospace">### <code r><br>library(circular)<br>nascente <- read.table(text=<br>"AMBIENTE    mes    angle    DISTANCIA    ABUND    RIQUEZA<span class=""><br>NASC        1    0    0    2    2<br>NASC        2    30    0    3    1<br>NASC        3    60    0    4    3<br>NASC        4    90    0    8    5<br>NASC        5    120    0    0    0<br>NASC        6    150    0    0    0<br>NASC        7    180    0    31    1<br>NASC        8    210    0    10    3<br>NASC        9    240    0    18    3<br>NASC        10    270    0    7    5<br>NASC        11    300    0    2    2<br></span>NASC        12    330    0    1    1", head=T)<br><br>nascente<br>nascente0=nascente[which(nascente$DISTANCIA==0),]; nascente0<br>meses <- c("Mai","Jun","Jul","Ago","Set","Out","Nov","Dez","Jan","Fev","Mar","Abr")<span class=""><br><br>#convertendo os dados para circular, usando radianos<br></span>nasc.c0 <- circular(rep(rad(nascente0$angle),nascente0$ABUND), units=c("radians"), rota="clock")<br>nasc.c0<br>plot(nasc.c0, stack=T)<br><br># Angulo médio + Desvio padrão angular<br>media <- mean(nasc.c0)<br>dp    <- sd(nasc.c0)<br>n     <- length(nasc.c0)<br>err   <- dp/sqrt(n)<br>len   <- rho.circular(nasc.c0)<br>cbind(media, dp, n, err, len)<br><br>sd(nasc.c0); sqrt(-2*log(rho.circular(nasc.c0)))<br><br># Gráfico<br># par(mfrow=c(1,3))<br># par(mar=c(0,0,2,0))<br>plot(nasc.c0,axes=F,main=NA)<br>axis.circular(at=circular(rad(nascente0$angle)), labels=meses, units=c("radians"), rota="clock")<br>rose.diag(nasc.c0, add=T, axes=F, col=8, bor=8, lwd=3)<br># arrows.circular(median(nasc.c0),y=rho.circular(nasc.c0),lwd=3,lty=1,angle=30,len=0.1,col="green") #aqui eu inseri a mediana<br>arrows.circular(mean(nasc.c0),y=len,   lwd=2, lty=1, angle=90, len=0.1, col=2) <br>arrows.circular(media+err,     y=1,    lwd=1, lty=3, angle=0, len=0.1, col=2) <br>arrows.circular(media-err,     y=1,    lwd=1, lty=3, angle=0, len=0.1, col=2) <br>### </code></span><span class=""><br><br><br clear="all"><div><div data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><span><br>================================================<br>Éder Comunello<br>Researcher at Brazilian Agricultural Research Corporation (Embrapa)<br>DSc in Agricultural Systems Engineering (USP/Esalq)<br>MSc in Environ. Sciences (UEM), <span>Agronomist (UEM)</span><br>---<br>Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>|<br>================================================<br>GEO, -22.2752, -54.8182, 408m<br>UTC-04:00 / DST: UTC-03:00</span><div><div><div><br></div><div><br></div></div><br></div></div></div></div></div></div></div></div></div></div></div></div><br></span></div></div>
</blockquote></div><br></div>