Cassiano, olá novamente!
Atualizei o código, pois percebi que havia um erro na atribuição dos labels
(meses) que estavam "counter-clock".
[image: Imagem inline 1]
### <code r>
library(circular)
nascente <- read.table(text=
"AMBIENTE mes angle DISTANCIA ABUND RIQUEZA
NASC 1 0 0 2 2
NASC 2 30 0 3 1
NASC 3 60 0 4 3
NASC 4 90 0 8 5
NASC 5 120 0 0 0
NASC 6 150 0 0 0
NASC 7 180 0 31 1
NASC 8 210 0 10 3
NASC 9 240 0 18 3
NASC 10 270 0 7 5
NASC 11 300 0 2 2
NASC 12 330 0 1 1", head=T)
nascente$meses <-
c("Mai","Jun","Jul","Ago","Set","Out","Nov","Dez","Jan","Fev","Mar","Abr")
nascente
#convertendo os dados para circular, usando radianos
nasc.c0 <- circular(rep(rad(nascente$angle),nascente$ABUND),
units=c("radians"), rota="clock")
nasc.c0
plot(nasc.c0, stack=T)
deg(nasc.c0)
# Angulo médio + Desvio padrão angular
media <- mean(nasc.c0)
dp <- sd(nasc.c0)
n <- length(nasc.c0)
err <- dp/sqrt(n)
len <- rho.circular(nasc.c0)
cbind(media, dp, n, err, len)
sd(nasc.c0); sqrt(-2*log(rho.circular(nasc.c0)))
# Gráfico
# par(mfrow=c(1,3))
# par(mar=c(0,0,2,0))
plot(nasc.c0,axes=F, main=NA, stack=F, col=3)
axis.circular(at=circular(rad(nascente$angle), rota="clock"),
labels=nascente$meses, units=c("radians"), col=4)
rose.diag(nasc.c0, add=T, axes=F, col=8, bor=8, lwd=3)
#
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
arrows.circular(mean(nasc.c0),y=len, lwd=2, lty=1, angle=90, len=0.1,
col=2)
arrows.circular(media+err, y=1, lwd=1, lty=3, angle=0, len=0.1,
col=2)
arrows.circular(media-err, y=1, lwd=1, lty=3, angle=0, len=0.1,
col=2)
nascente
### </code>
================================================
Éder Comunello
Researcher at Brazilian Agricultural Research Corporation (Embrapa)
DSc in Agricultural Systems Engineering (USP/Esalq)
MSc in Environ. Sciences (UEM), Agronomist (UEM)
---
Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>|
================================================
GEO, -22.2752, -54.8182, 408m
UTC-04:00 / DST: UTC-03:00
Em 9 de junho de 2016 08:41, Éder Comunello <[email protected]>
escreveu:
> Cassiano, bom dia!
>
> 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:
>
>
> Infelizmente, não encontrei como reproduzir esse efeito no R. O melhor que
> consegui foi marcar os limites com novas linhas.
> [image: Imagem inline 1]
>
>
> ### <code r>
> library(circular)
> nascente <- read.table(text=
> "AMBIENTE mes angle DISTANCIA ABUND RIQUEZA
> NASC 1 0 0 2 2
> NASC 2 30 0 3 1
> NASC 3 60 0 4 3
> NASC 4 90 0 8 5
> NASC 5 120 0 0 0
> NASC 6 150 0 0 0
> NASC 7 180 0 31 1
> NASC 8 210 0 10 3
> NASC 9 240 0 18 3
> NASC 10 270 0 7 5
> NASC 11 300 0 2 2
> NASC 12 330 0 1 1", head=T)
>
> nascente
> nascente0=nascente[which(nascente$DISTANCIA==0),]; nascente0
> meses <-
> c("Mai","Jun","Jul","Ago","Set","Out","Nov","Dez","Jan","Fev","Mar","Abr")
>
> #convertendo os dados para circular, usando radianos
> nasc.c0 <- circular(rep(rad(nascente0$angle),nascente0$ABUND),
> units=c("radians"), rota="clock")
> nasc.c0
> plot(nasc.c0, stack=T)
>
> # Angulo médio + Desvio padrão angular
> media <- mean(nasc.c0)
> dp <- sd(nasc.c0)
> n <- length(nasc.c0)
> err <- dp/sqrt(n)
> len <- rho.circular(nasc.c0)
> cbind(media, dp, n, err, len)
>
> sd(nasc.c0); sqrt(-2*log(rho.circular(nasc.c0)))
>
> # Gráfico
> # par(mfrow=c(1,3))
> # par(mar=c(0,0,2,0))
> plot(nasc.c0,axes=F,main=NA)
> axis.circular(at=circular(rad(nascente0$angle)), labels=meses,
> units=c("radians"), rota="clock")
> rose.diag(nasc.c0, add=T, axes=F, col=8, bor=8, lwd=3)
> #
> 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
> arrows.circular(mean(nasc.c0),y=len, lwd=2, lty=1, angle=90, len=0.1,
> col=2)
> arrows.circular(media+err, y=1, lwd=1, lty=3, angle=0, len=0.1,
> col=2)
> arrows.circular(media-err, y=1, lwd=1, lty=3, angle=0, len=0.1,
> col=2)
> ### </code>
>
>
>
> ================================================
> Éder Comunello
> Researcher at Brazilian Agricultural Research Corporation (Embrapa)
> DSc in Agricultural Systems Engineering (USP/Esalq)
> MSc in Environ. Sciences (UEM), Agronomist (UEM)
> ---
> Embrapa Agropecuária Oeste, Dourados, MS, Brazil |<O>|
> ================================================
> GEO, -22.2752, -54.8182, 408m
> UTC-04:00 / DST: UTC-03:00
>
>
>
>
>
_______________________________________________
R-br mailing list
[email protected]
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.
