<div dir="ltr">Segue:<div><br></div><div>Dados</div><div><div>AMBIENTE<span class="" style="white-space:pre"> </span>mess<span class="" style="white-space:pre"> </span>angle<span class="" style="white-space:pre"> </span>DISTANCIA<span class="" style="white-space:pre"> </span>ABUND<span class="" style="white-space:pre"> </span>RIQUEZA</div><div>NASC <span class="" style="white-space:pre"> </span>1<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>2<span class="" style="white-space:pre"> </span>2</div><div>NASC <span class="" style="white-space:pre"> </span>2<span class="" style="white-space:pre"> </span>30<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>3<span class="" style="white-space:pre"> </span>1</div><div>NASC <span class="" style="white-space:pre"> </span>3<span class="" style="white-space:pre"> </span>60<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>4<span class="" style="white-space:pre"> </span>3</div><div>NASC <span class="" style="white-space:pre"> </span>4<span class="" style="white-space:pre"> </span>90<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>8<span class="" style="white-space:pre"> </span>5</div><div>NASC <span class="" style="white-space:pre"> </span>5<span class="" style="white-space:pre"> </span>120<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>0</div><div>NASC <span class="" style="white-space:pre"> </span>6<span class="" style="white-space:pre"> </span>150<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>0</div><div>NASC <span class="" style="white-space:pre"> </span>7<span class="" style="white-space:pre"> </span>180<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>31<span class="" style="white-space:pre"> </span>1</div><div>NASC <span class="" style="white-space:pre"> </span>8<span class="" style="white-space:pre"> </span>210<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>10<span class="" style="white-space:pre"> </span>3</div><div>NASC <span class="" style="white-space:pre"> </span>9<span class="" style="white-space:pre"> </span>240<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>18<span class="" style="white-space:pre"> </span>3</div><div>NASC <span class="" style="white-space:pre"> </span>10<span class="" style="white-space:pre"> </span>270<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>7<span class="" style="white-space:pre"> </span>5</div><div>NASC <span class="" style="white-space:pre"> </span>11<span class="" style="white-space:pre"> </span>300<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>2<span class="" style="white-space:pre"> </span>2</div><div>NASC <span class="" style="white-space:pre"> </span>12<span class="" style="white-space:pre"> </span>330<span class="" style="white-space:pre"> </span>0<span class="" style="white-space:pre"> </span>1<span class="" style="white-space:pre"> </span>1</div><div><br></div><div>###############################################################</div><div>################### início da análise ##############################</div><div>###############################################################</div><div><br></div><div><div>library(circular)</div><div>nascente=read.table("nasc.csv",h=T)</div><div>nascente</div><div>nascente0=nascente[which(nascente$DISTANCIA==0),]</div><div>nascente0</div><div><br></div><div><br></div><div>meses<-c("Mai","Jun","Jul","Ago","Set","Out","Nov","Dez","Jan","Fev","Mar","Abr")</div><div><br></div><div><br></div><div><br></div><div><br></div><div>#convertendo os dados para circular, usando radianos</div><div>nasc.c0<-circular(rep(rad(nascente0$angle),nascente0$ABUND),units=c("radians"),rotation="clock")</div><div><br></div><div><br></div><div>#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% Angulo médio + Desvio padrão angular%%%%%%%%%%%%%%%%%%%%</div><div> # %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%</div><div> media=mean(nasc.c0)</div><div>sqrt(-2*log(rho.circular(nasc.c0)))</div><div>#%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%<br></div></div><div><br></div><div><br></div><div>#### Gráfico ############</div><div><div>par(mfrow=c(1,3))</div><div>par(mar=c(0,0,2,0))</div><div>plot(nasc.c0,axes=F,main="0")</div><div>axis.circular(at=circular(rad(nascente0$angle)),labels=meses,units=c("radians"),</div><div> rotation="clock")</div><div>rose.diag(nasc.c0,add=T,axes=F)</div><div><div>arrows.circular(median(nasc.c0),y=rho.circular(nasc.c0),lwd=4,lty=1,angle=90,len=0.0,col="green") #aqui eu inseri a mediana</div><div>arrows.circular(mean(nasc.c0),y=rho.circular(nasc.c0),lwd=4,lty=1,angle=90,len=0.0,col="red") #aqui eu inseri a média</div></div></div><div><br></div><div>#agora preciso inserir os desvios, ou erros padrão.</div><div><br></div><div>########################### FIM ######################################</div><div><br></div><div>Obrigado</div><div><br></div><div><br></div><div><br></div><div class="gmail_extra"><br><div class="gmail_quote">Em 8 de junho de 2016 12:00, <span dir="ltr"><<a href="mailto:r-br-request@listas.c3sl.ufpr.br" target="_blank">r-br-request@listas.c3sl.ufpr.br</a>></span> escreveu:<br><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex">Enviar submissões para a lista de discussão R-br para<br>
<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
<br>
Para se cadastrar ou descadastrar via WWW, visite o endereço<br>
<a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
ou, via email, envie uma mensagem com a palavra 'help' no assunto ou<br>
corpo da mensagem para<br>
<a href="mailto:r-br-request@listas.c3sl.ufpr.br">r-br-request@listas.c3sl.ufpr.br</a><br>
<br>
Você poderá entrar em contato com a pessoa que gerencia a lista pelo<br>
endereço<br>
<a href="mailto:r-br-owner@listas.c3sl.ufpr.br">r-br-owner@listas.c3sl.ufpr.br</a><br>
<br>
Quando responder, por favor edite sua linha Assunto assim ela será<br>
mais específica que "Re: Contents of R-br digest..."<br>
<br>
<br>
Tópicos de Hoje:<br>
<br>
1. Re: Estatística circular - gráfico (Éder Comunello)<br>
2. Erro modelo - glm (Andre Oliveira)<br>
3. Re: Erro modelo - glm (Marcus Nunes)<br>
4. Re: Erro modelo - glm (Leonardo Ferreira Fontenelle)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Tue, 7 Jun 2016 14:39:28 -0400<br>
From: Éder Comunello <<a href="mailto:comunello.eder@gmail.com">comunello.eder@gmail.com</a>><br>
To: "<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>" <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: Re: [R-br] Estatística circular - gráfico<br>
Message-ID:<br>
<CABmC8gn-f8sdoHs0-Hc=zAGTmd=iYjMHizFe=<a href="mailto:X_HK6JYa4Vobw@mail.gmail.com">X_HK6JYa4Vobw@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Cassiano, boa tarde!<br>
<br>
Talvez seja o caso de postar seu código com um dataset de exemplo.<br>
<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), Agronomist (UEM)<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<br>
<br>
<br>
<br>
<br>
Em 7 de junho de 2016 10:11, Cassiano <<a href="mailto:cassianosr@gmail.com">cassianosr@gmail.com</a>> escreveu:<br>
<br>
> Obrigado pela ajuda.<br>
> Realmente o ggplot2 produz gráficos mais interessantes. Mas ainda não<br>
> achei a opção para plotar a média e os desvios angulares.<br>
> Se puderem me dar mais instruções, ficarei agradecido.<br>
><br>
> Abraços<br>
> Cassiano<br>
><br>
><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160607/2bc39a6a/attachment-0001.html" rel="noreferrer" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160607/2bc39a6a/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 2<br>
Date: Tue, 7 Jun 2016 23:26:58 +0000 (UTC)<br>
From: Andre Oliveira <<a href="mailto:andreolsouza@yahoo.com.br">andreolsouza@yahoo.com.br</a>><br>
To: R-br Lista <<a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>><br>
Subject: [R-br] Erro modelo - glm<br>
Message-ID:<br>
<<a href="mailto:1694617424.92644.1465342018577.JavaMail.yahoo@mail.yahoo.com">1694617424.92644.1465342018577.JavaMail.yahoo@mail.yahoo.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
boa noite, alguém poderia ajudar descobrir onde está o erro?<br>
<br>
remove(list=ls())diurno=c(68,36,40,30)<br>
noturno=c(46,26,55,22)<br>
quadra= c("A", "B", "G","D")<br>
dados=data.frame(diurno,noturno,quadra)<br>
fit = glm(cbind(diurno, noturno) ~ factor(quadra), family = poisson, data=dados)<br>
Error in x[good, , drop = FALSE] :<br>
(subscript) subscrito lógico muito longo<br>
obrigado<br>
André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES<br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160607/7d2e5cef/attachment-0001.html" rel="noreferrer" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160607/7d2e5cef/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 3<br>
Date: Tue, 7 Jun 2016 22:37:17 -0300<br>
From: Marcus Nunes <<a href="mailto:marcus.nunes@gmail.com">marcus.nunes@gmail.com</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a>, Andre Oliveira<br>
<<a href="mailto:andreolsouza@yahoo.com.br">andreolsouza@yahoo.com.br</a>><br>
Subject: Re: [R-br] Erro modelo - glm<br>
Message-ID:<br>
<<a href="mailto:CA%2BQGQvvGr8Rg5w%2BmO4tNmRiG1du%2BnF0Z83XCZUrDfYyhgZdTww@mail.gmail.com">CA+QGQvvGr8Rg5w+mO4tNmRiG1du+nF0Z83XCZUrDfYyhgZdTww@mail.gmail.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Ajuste os modelos individualmente:<br>
<br>
diurno=c(68,36,40,30)<br>
noturno=c(46,26,55,22)<br>
quadra= c("A", "B", "G","D")<br>
dados=data.frame(diurno,noturno,quadra)<br>
fit1 = glm(diurno ~ factor(quadra), family = poisson, data=dados)<br>
fit2 = glm(noturno ~ factor(quadra), family = poisson, data=dados)<br>
<br>
Isto deve resolver o teu problema.<br>
<br>
<br>
<br>
2016-06-07 20:26 GMT-03:00 Andre Oliveira <<a href="mailto:andreolsouza@yahoo.com.br">andreolsouza@yahoo.com.br</a>>:<br>
<br>
> boa noite, alguém poderia ajudar descobrir onde está o erro?<br>
><br>
> remove(list=ls())<br>
> diurno=c(68,36,40,30)<br>
> noturno=c(46,26,55,22)<br>
> quadra= c("A", "B", "G","D")<br>
> dados=data.frame(diurno,noturno,quadra)<br>
> fit = glm(cbind(diurno, noturno) ~ factor(quadra), family = poisson,<br>
> data=dados)<br>
><br>
> *Error in x[good, , drop = FALSE] : *<br>
> * (subscript) subscrito lógico muito longo*<br>
><br>
> obrigado<br>
><br>
> André Oliveira Souza. Graduação em Matemática, mestrado em estatística<br>
> aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito<br>
> Santo. IFES<br>
><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" rel="noreferrer" 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" rel="noreferrer" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça<br>
> código mínimo reproduzível.<br>
><br>
<br>
<br>
<br>
--<br>
Marcus Nunes<br>
<a href="http://marcusnunes.me/" rel="noreferrer" target="_blank">http://marcusnunes.me/</a><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160607/a1cd75b1/attachment-0001.html" rel="noreferrer" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160607/a1cd75b1/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Message: 4<br>
Date: Wed, 08 Jun 2016 11:16:04 -0300<br>
From: Leonardo Ferreira Fontenelle <<a href="mailto:leonardof@leonardof.med.br">leonardof@leonardof.med.br</a>><br>
To: <a href="mailto:r-br@listas.c3sl.ufpr.br">r-br@listas.c3sl.ufpr.br</a><br>
Subject: Re: [R-br] Erro modelo - glm<br>
Message-ID:<br>
<<a href="mailto:1465395364.4116631.631632273.0B1F3BA2@webmail.messagingengine.com">1465395364.4116631.631632273.0B1F3BA2@webmail.messagingengine.com</a>><br>
Content-Type: text/plain; charset="utf-8"<br>
<br>
Em Ter 7 jun. 2016, às 20:26, Andre Oliveira escreveu:<br>
> boa noite, alguém poderia ajudar descobrir onde está o erro?<br>
><br>
> remove(list=ls())<br>
> diurno=c(68,36,40,30)<br>
> noturno=c(46,26,55,22)<br>
> quadra= c("A", "B", "G","D")<br>
> dados=data.frame(diurno,noturno,quadra)<br>
> fit = glm(cbind(diurno, noturno) ~ factor(quadra), family = poisson,<br>
> data=dados)<br>
> *Error in x[good, , drop = FALSE] :<br>
*<br>
> * (subscript) subscrito lógico muito longo*<br>
<br>
Andre, conforme está escrito na documentação do glm, a resposta pode ser<br>
fornecida como duas colunas (sucessos e fracassos) *apenas para<br>
regressões binomiais e quasi-binomiais*; não para Poisson.<br>
<br>
Supondo que seus dados sejam o resumo de 68 + 46 = 114 observações para<br>
a quadra A, 36 + 26 = 62 observações para a quadra B, e por aí em<br>
diante; e que você queira calcular uma razão de prevalência / risco /<br>
outra proporção entre as quadras, o jeito mais óbvio de fornecer seus<br>
dados para a função é em formato longo:<br>
<br>
dados <- data.frame(horariodiurno = c(1, 0, 1, 0, 1, 0, 1, 0), quadra =<br>
c("A", "A", "B", "B", "C", "C", "D", "D"))<br>
dados <- dados[rep(1:nrow(dados), c(68, 46, 36, 26, 40, 55, 30, 22)), ]<br>
fit = glm(horariodiurno ~ quadra, family = quasipoisson, data = dados)<br>
<br>
Repare também que, se você quer trabalhar com razão de proporções, a<br>
distribuição Poisson em princípio não se aplica, daí a necessidade de<br>
você especificar uma quasipoisson ou então um estimador de variância<br>
consistente com heterocedasticidade (pacote sandwich).<br>
<br>
Por favor, use o código com cuidado porque não estou com meu computador<br>
aqui para testá-lo.<br>
<br>
Atenciosamente,<br>
<br>
Leonardo Ferreira Fontenelle[1]<br>
<br>
Links:<br>
<br>
1. <a href="http://lattes.cnpq.br/9234772336296638" rel="noreferrer" target="_blank">http://lattes.cnpq.br/9234772336296638</a><br>
-------------- Próxima Parte ----------<br>
Um anexo em HTML foi limpo...<br>
URL: <<a href="http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160608/7af199f3/attachment-0001.html" rel="noreferrer" target="_blank">http://listas.inf.ufpr.br/pipermail/r-br/attachments/20160608/7af199f3/attachment-0001.html</a>><br>
<br>
------------------------------<br>
<br>
Subject: Legenda do Digest<br>
<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" rel="noreferrer" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>
<br>
<br>
------------------------------<br>
<br>
Fim da Digest R-br, volume 66, assunto 7<br>
****************************************<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div style="color:rgb(255,0,0)"><span style="color:rgb(0,0,0)">=======================================</span></div>Cassiano S. Rosa<div><div style="font-family:arial;font-size:small">CRBio: 104117/04-D</div><div style="font-family:arial;font-size:small"><br></div><div dir="ltr" style="font-family:arial;font-size:small"><div style="font-family:arial,sans-serif;font-size:12.8px;line-height:21.3px;color:rgb(68,68,68)"><span style="line-height:22.72px;text-align:center"><font style="line-height:normal"><div style="line-height:21.3px;display:inline!important"><span style="line-height:22.72px"><font face="trebuchet ms, sans-serif" style="line-height:normal"><div style="line-height:21.3px;display:inline!important">Universidade Federal do Triângulo Mineiro (UFTM), </div></font></span></div></font></span></div><div style="font-family:arial,sans-serif;font-size:12.8px;line-height:21.3px;color:rgb(68,68,68)"><span style="line-height:22.72px;text-align:center"><font style="line-height:normal"><div style="line-height:21.3px;display:inline!important"><span style="line-height:22.72px"><font face="trebuchet ms, sans-serif" style="line-height:normal"><div style="line-height:21.3px;display:inline!important"><i>campus</i> Iturama-MG</div></font></span></div></font></span></div><div style="font-family:arial,sans-serif;font-size:12.8px;line-height:21.3px;color:rgb(68,68,68);text-align:center"><font face="trebuchet ms, sans-serif"><span style="line-height:22.72px"><font style="line-height:normal"><div style="text-align:left;line-height:21.3px"><span style="line-height:21.3px;font-size:12.8px">Avenida Rio Paranaíba, 1241, Centro, 38280-000</span><br></div><div style="text-align:left;line-height:21.3px"><span style="line-height:21.3px;font-size:12.8px">Tel: (34)3415-2512</span></div></font></span></font></div></div><div style="text-align:justify;text-indent:30px"><font color="#000000" face="Verdana, Arial, Helvetica, sans-serif"><span style="font-size:11.8182px"><br></span></font></div><div><div><div><span style="font-family:Arial,Helvetica,sans-serif"><span style="text-align:left;color:rgb(102,102,102);vertical-align:top"><a href="http://lattes.cnpq.br/3030825329622014" target="_blank"><font size="2">Lattes</font></a></span></span></div><div><a href="http://scholar.google.com/citations?user=JmELzOMAAAAJ" target="_blank"><font size="2">Scholar Google</font></a></div></div><div><font size="2"><a href="https://uftm.academia.edu/CassianoRosa" target="_blank">Academia.edu</a></font><br></div><div><a href="https://www.researchgate.net/profile/Cassiano_Rosa" target="_blank">Researchgate</a><br></div><div><div><br>=======================================<br></div></div></div></div></div></div></div></div></div></div></div></div>
</div></div></div>