<html><body><div style="color:#000; background-color:#fff; font-family:HelveticaNeue-Light, Helvetica Neue Light, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif;font-size:10px"><div id="yui_3_16_0_1_1413252330318_3036">Obrigado pessoal.</div><div id="yui_3_16_0_1_1413252330318_3044"><br></div><div id="yui_3_16_0_1_1413252330318_2550" dir="ltr">Me parece que: <br></div><div id="yui_3_16_0_1_1413252330318_2549" dir="ltr">require(e1071)<br style="" class="">skewness(x, type=1)<br style="" class="">kurtosis(x,type=1)</div><div id="yui_3_16_0_1_1413252330318_2548" dir="ltr"><br></div><div id="yui_3_16_0_1_1413252330318_3040" dir="ltr">atende aos estimadores..</div><div id="yui_3_16_0_1_1413252330318_3041" dir="ltr"><br></div><div id="yui_3_16_0_1_1413252330318_2547" dir="ltr">curtose=m_4/s^4<br style="" class="" clear="none">
assimetria= m_3/s^3</div><div id="yui_3_16_0_1_1413252330318_2539"><span></span></div><div id="yui_3_16_0_1_1413252330318_2540"> </div><div id="yui_3_16_0_1_1413252330318_2542"><div id="yui_3_16_0_1_1413252330318_2541">
<br>
<br>André Oliveira Souza.
<br>Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES
<br>
<br></div></div> <div id="yui_3_16_0_1_1413252330318_3042" class="qtdSeparateBR"><br><br></div><div id="yui_3_16_0_1_1413252330318_2882" style="display: block;" class="yahoo_quoted"> <div id="yui_3_16_0_1_1413252330318_2881" style="font-family: HelveticaNeue-Light, Helvetica Neue Light, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 10px;"> <div id="yui_3_16_0_1_1413252330318_2880" style="font-family: HelveticaNeue, Helvetica Neue, Helvetica, Arial, Lucida Grande, sans-serif; font-size: 16px;"> <div id="yui_3_16_0_1_1413252330318_3043" dir="ltr"> <font id="yui_3_16_0_1_1413252330318_3045" face="Arial" size="2"> Em Domingo, 12 de Outubro de 2014 19:58, Ricardo Fernando Zampieri <rfzampieri@gmail.com> escreveu:<br> </font> </div> <br><br> <div id="yui_3_16_0_1_1413252330318_2879" class="y_msg_container"><div id="yiv8462273685"><div id="yui_3_16_0_1_1413252330318_2878"><div id="yui_3_16_0_1_1413252330318_2877" dir="ltr">Cada pacote calcula de modo diferente (e1071, moments). Sugiro dar uma olhada nos pdfs que explicam os comandos de cada pacote. Veja estas anotações:<br clear="none"><br clear="none">#A. CURTOSE<br clear="none"><br clear="none">#Analisa a proximidade ou o afastamento de um conjunto de dados quantitativos em relação à distribuição normal (Gaussiana)<br clear="none">#Caracteriza o achatamento ou o alongamento da curva da função de distribuição de probabilidade<br clear="none">#Fórmula convencional para se calcular a curtose: C= Q3 - Q1 / 2(P90 - P10)<br clear="none"> #mede o grau de achatamento ou alongamento da curva<br clear="none"> #quanto > valor da curtose > achatamento da curva<br clear="none">#Formas da curva da função de distribuição de probabilidade:<br clear="none"> #mesocúrtica: mesma forma que a distribuição normal (C= 0,263)<br clear="none"> #platicúrtica: maior achatamento que a distribuição normal (C> 0,263)<br clear="none"> #leptocúrtica: mais alta e concentrada que a distribuição normal (C< 0,263)<br clear="none">#Para se calcular a curtose: pacotes moments e e1071<br clear="none">#Comando kurtosis para se calcular a curtose (moments) ou o excesso de curtose (e1071)<br clear="none">#Library(moments): comando library para se carregar pacotes; moments necessário para se calcular a curtose (conforme a fórmula deste pacote); moments: <a rel="nofollow" shape="rect" target="_blank" href="http://cran.r-project.org/web/packages/moments/moments.pdf">http://cran.r-project.org/web/packages/moments/moments.pdf</a><br clear="none"> #pode-se obter: C=3 (mesocúrtica) C>3 (platicúrtica) C<3 (leptocúrtica)<br clear="none">#Library(e1071): comando library para se carregar pacotes; e1071 necessário para se calcular o excesso de curtose (conforme a fórmula deste pacote); e1071: <a rel="nofollow" shape="rect" target="_blank" href="http://cran.r-project.org/web/packages/e1071/e1071.pdf">http://cran.r-project.org/web/packages/e1071/e1071.pdf</a><br clear="none"> #há três tipos de cálculos estatísticos no e1071: kurtosis(Renda,type=1); kurtosis(Renda,type=2); kurtosis(Renda,type=3); kurtosis(Renda)=kurtosis(Renda,type=3)<br clear="none"> #kurtosis(Renda,type=1) retorna o mesmo valor de moments-3<br clear="none"> #no pacote e1071, a distribuição normal padrão apresenta referencialmente o valor de curtose=3 (em vez de 0,263)<br clear="none"> #calcula-se então o excesso de curtose por meio do comando kurtosis:<br clear="none"> #ExCur=C-3 (o programa calcula o valor de C, conforme a fórmula deste pacote, e o subtrai de 3)<br clear="none"> #pode-se obter: ExCur=0 (mesocúrtica) ExCur>0 (platicúrtica) ExCur<0 (leptocúrtica)<br clear="none"> #Em distribuições normais, independentemente do número de dados amostrais, o excesso de curtose terá valor próximo de 0<br clear="none"><br clear="none">median(Renda) #Comando median para se calcular a mediana<br clear="none">quantile(Renda,prob=c(0.75,0.25,0.9,0.1,0.5),type=2) #Argumento prob para se visualizar valores correspondentes às probabilidades especificadas (dados necessários para se calcular o coeficiente de curtose pela fórmula convencional)<br clear="none">k1<-(1955-950)/(2*(2155-920)) #Fórmula convencional para se calcular a curtose<br clear="none">k1 #Para se visualizar curtose<br clear="none"><br clear="none">library(moments) #Comando library para se carregar pacotes; moments necessário para se calcular a curtose (conforme a fórmula deste pacote); moments: <a rel="nofollow" shape="rect" target="_blank" href="http://cran.r-project.org/web/packages/moments/moments.pdf">http://cran.r-project.org/web/packages/moments/moments.pdf</a><br clear="none">kurtosis(Renda) #Comando kurtosis para se calcular a curtose<br clear="none"><br clear="none">library(e1071) #Comando library para se carregar pacotes; e1071 necessário para se calcular o excesso de curtose (conforme a fórmula deste pacote); e1071: <a rel="nofollow" shape="rect" target="_blank" href="http://cran.r-project.org/web/packages/e1071/e1071.pdf">http://cran.r-project.org/web/packages/e1071/e1071.pdf</a><br clear="none">kurtosis(Renda,type=1) #Comando kurtosis para se calcular o excesso de curtose<br clear="none">kurtosis(Renda,type=2) #Comando kurtosis para se calcular o excesso de curtose<br clear="none">kurtosis(Renda,type=3) #Comando kurtosis para se calcular o excesso de curtose<br clear="none"> <br clear="none">#B. ASSIMETRIA<br clear="none"><br clear="none">#Analisa a proximidade ou o afastamento de um conjunto de dados quantitativos em relação à distribuição normal (Gaussiana)<br clear="none">#Mede o grau de afastamento de uma distribuição em relação a um eixo central (geralmente a média), podendo-se caracterizar curvas assimétricas à direita ou à esquerda<br clear="none">#Pearson descreveu duas fórmulas para se calcular o coeficiente de assimetria:<br clear="none"> #primeiro coeficiente de assimetria de Pearson: AS= média - moda / desvio padrão (conjunto de dados não pode ser plurimodal)<br clear="none"> #segundo coeficiente de assimetria de Pearson: AS= Q1 + Q3 - 2Q2 / Q3 - Q1<br clear="none">#Classificação das distribuições de acordo com o coeficiente de assimetria de Pearson: p. 16-18 Apostila Estatística - Graduação UEL Parte 01<br clear="none"> #simétrica: AS=0 (média, mediana e moda coincidem)<br clear="none"> #assimétrica negativa: AS<0 (cauda da distribuição do lado esquerdo do gráfico); média deslocada para o lado esquerdo<br clear="none"> #assimétrica positiva: AS>0 (cauda da distribuição do lado direito do gráfico); média deslocada para o lado direito<br clear="none">#Comando skewness para se calcular a assimetria<br clear="none">#Library(moments): comando library para se carregar pacotes; moments necessário para se calcular a assimetria (conforme a fórmula deste pacote); moments: <a rel="nofollow" shape="rect" target="_blank" href="http://cran.r-project.org/web/packages/moments/moments.pdf">http://cran.r-project.org/web/packages/moments/moments.pdf</a><br clear="none">#Library(e1071): comando library para se carregar pacotes; e1071 necessário para se calcular a assimetria (conforme a fórmula deste pacote); e1071: <a rel="nofollow" shape="rect" target="_blank" href="http://cran.r-project.org/web/packages/e1071/e1071.pdf">http://cran.r-project.org/web/packages/e1071/e1071.pdf</a><br clear="none"> #há três tipos de cálculos estatísticos no e1071: skewness(Renda,type=1); skewness(Renda,type=2); skewness(Renda,type=3); skewness(Renda)=skewness(Renda,type=3)<br clear="none"><br clear="none">median(Renda) #Comando median para se calcular a mediana<br clear="none">quantile(Renda,prob=c(0.25,0.5,0.75),type=2) #Argumento prob para se visualizar valores correspondentes às probabilidades especificadas (dados necessários para se calcular o segundo coeficiente de assimetria de Pearson)<br clear="none">asrenda<-(950+1955-(2*1120))/(1955-950) #Fórmula para se calcular o segundo coeficiente de assimetria de Pearson<br clear="none">asrenda #Para se visualizar o segundo coeficiente de assimetria de Pearson<br clear="none"><br clear="none">library(moments) #Comando library para se carregar pacotes; moments necessário para se calcular a assimetria (conforme a fórmula deste pacote); moments: <a rel="nofollow" shape="rect" target="_blank" href="http://cran.r-project.org/web/packages/moments/moments.pdf">http://cran.r-project.org/web/packages/moments/moments.pdf</a><br clear="none">skewness(Renda) #Comando skewness para se calcular a assimetria<br clear="none"><br clear="none">library(e1071) #Comando library para se carregar pacotes; e1071 necessário para se calcular a assimetria (conforme a fórmula deste pacote); e1071: <a rel="nofollow" shape="rect" target="_blank" href="http://cran.r-project.org/web/packages/e1071/e1071.pdf">http://cran.r-project.org/web/packages/e1071/e1071.pdf</a><br clear="none">skewness(Renda,type=1) #Comando skewness para se calcular a assimetria<br clear="none">skewness(Renda,type=2) #Comando skewness para se calcular a assimetria<br clear="none">skewness(Renda,type=3) #Comando skewness para se calcular a assimetria<br clear="none"><br clear="none"><br clear="none"><br clear="none"></div><div id="yui_3_16_0_1_1413252330318_2888" class="yiv8462273685gmail_extra"><br clear="none"><div id="yui_3_16_0_1_1413252330318_2887" class="yiv8462273685gmail_quote">Em 12 de outubro de 2014 16:15, Daniel Marcelino <span dir="ltr"><<a rel="nofollow" shape="rect" ymailto="mailto:dmsilva.br@gmail.com" target="_blank" href="mailto:dmsilva.br@gmail.com">dmsilva.br@gmail.com</a>></span> escreveu:<br clear="none"><blockquote id="yui_3_16_0_1_1413252330318_2886" class="yiv8462273685gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">André, eu acho que há muitas opções disponíveis em que permite você<br clear="none">
escolher o tipo de curtose você quer. Veja por exemplo no pacote:<br clear="none">
e1071<br clear="none">
<br clear="none">
x=c(2,3,NA,5,6,8,9)<br clear="none">
<br clear="none">
require(e1071)<br clear="none">
<br clear="none">
> kurtosis(x, na.rm = TRUE,type=3)<br clear="none">
[1] -1.878889<br clear="none">
> kurtosis(x, na.rm = TRUE,type=2)<br clear="none">
[1] -1.541333<br clear="none">
> kurtosis(x, na.rm = TRUE,type=1)<br clear="none">
[1] -1.3856<br clear="none">
<br clear="none">
require(SciencesPo)<br clear="none">
Tipo 2 como padrão.<br clear="none">
> detail(x)<br clear="none">
vars obs mean sd var se median mad trimmed winsor range min max<br clear="none">
1 1 6 5.5 2.7 7.5 1.1 5.5 3.7 5.5 5.5 7 2 9<br clear="none">
skew kurt<br clear="none">
1 0 -1.5<br clear="none">
<br clear="none">
2014-10-12 14:14 GMT-03:00 Andre Oliveira <<a rel="nofollow" shape="rect" ymailto="mailto:andreolsouza@yahoo.com.br" target="_blank" href="mailto:andreolsouza@yahoo.com.br">andreolsouza@yahoo.com.br</a>>:<br clear="none">
> Pessoal alguém tem alguma sugestão melhor para estimar assimetria e curtose?<br clear="none">
><br clear="none">
> Este pacote não bate a definição de Ferreira, D, F. 2005.<br clear="none">
><br clear="none">
> curtose=m_4/s^4<br clear="none">
> assimetria= m_3/s^3<br clear="none">
><br clear="none">
> --------------------------------------------------------------------------------------------------------------------------------------------------------<br clear="none">
> x=c(2,3,,5,6,8,9)<br clear="none">
> library(fBasics)<br clear="none">
> basicStats(x)<br clear="none">
><br clear="none">
><br clear="none">
><br clear="none">
> André Oliveira Souza.<br clear="none">
> Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal<br clear="none">
> de Educação, Ciência e Tecnologia do Espirito Santo. IFES<br clear="none">
><br clear="none">
><br clear="none">
> _______________________________________________<br clear="none">
> R-br mailing list<br clear="none">
> <a rel="nofollow" shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none">
> <a rel="nofollow" shape="rect" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">
> Leia o guia de postagem (<a rel="nofollow" shape="rect" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código<br clear="none">
> mínimo reproduzível.<div class="yiv8462273685yqt8823400124" id="yiv8462273685yqtfd35638"><br clear="none">
_______________________________________________<br clear="none">
R-br mailing list<br clear="none">
<a rel="nofollow" shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" target="_blank" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none">
<a rel="nofollow" shape="rect" target="_blank" href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br clear="none">
Leia o guia de postagem (<a rel="nofollow" shape="rect" target="_blank" href="http://www.leg.ufpr.br/r-br-guia">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.</div><br clear="none">
</blockquote></div><br clear="none"><br clear="all"><br clear="none">-- <br clear="none"> \\\ /// <br clear="none"> ( @ @ ) <br clear="none"> ...o00o.(_).o00o...<br clear="none"> Ricardo
</div></div></div><br><div class="yqt8823400124" id="yqtfd52924">_______________________________________________<br clear="none">R-br mailing list<br clear="none"><a shape="rect" ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br clear="none"><a shape="rect" 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 clear="none">Leia o guia de postagem (<a shape="rect" 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.</div><br><br></div> </div> </div> </div> </div></body></html>