<div 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><br>#A. CURTOSE<br><br>#Analisa a proximidade ou o afastamento de um conjunto de dados quantitativos em relação à distribuição normal (Gaussiana)<br>#Caracteriza o achatamento ou o alongamento da curva da função de distribuição de probabilidade<br>#Fórmula convencional para se calcular a curtose: C= Q3 - Q1 / 2(P90 - P10)<br> #mede o grau de achatamento ou alongamento da curva<br> #quanto > valor da curtose > achatamento da curva<br>#Formas da curva da função de distribuição de probabilidade:<br> #mesocúrtica: mesma forma que a distribuição normal (C= 0,263)<br> #platicúrtica: maior achatamento que a distribuição normal (C> 0,263)<br> #leptocúrtica: mais alta e concentrada que a distribuição normal (C< 0,263)<br>#Para se calcular a curtose: pacotes moments e e1071<br>#Comando kurtosis para se calcular a curtose (moments) ou o excesso de curtose (e1071)<br>#Library(moments): comando library para se carregar pacotes; moments necessário para se calcular a curtose (conforme a fórmula deste pacote); moments: <a href="http://cran.r-project.org/web/packages/moments/moments.pdf">http://cran.r-project.org/web/packages/moments/moments.pdf</a><br> #pode-se obter: C=3 (mesocúrtica)  C>3 (platicúrtica)  C<3 (leptocúrtica)<br>#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 href="http://cran.r-project.org/web/packages/e1071/e1071.pdf">http://cran.r-project.org/web/packages/e1071/e1071.pdf</a><br> #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> #kurtosis(Renda,type=1) retorna o mesmo valor de moments-3<br> #no pacote e1071, a distribuição normal padrão apresenta referencialmente o valor de curtose=3 (em vez de 0,263)<br> #calcula-se então o excesso de curtose por meio do comando kurtosis:<br>  #ExCur=C-3 (o programa calcula o valor de C, conforme a fórmula deste pacote, e o subtrai de 3)<br>  #pode-se obter: ExCur=0 (mesocúrtica)  ExCur>0 (platicúrtica)  ExCur<0 (leptocúrtica)<br> #Em distribuições normais, independentemente do número de dados amostrais, o excesso de curtose terá valor próximo de 0<br><br>median(Renda) #Comando median para se calcular a mediana<br>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>k1<-(1955-950)/(2*(2155-920)) #Fórmula convencional para se calcular a curtose<br>k1 #Para se visualizar curtose<br><br>library(moments) #Comando library para se carregar pacotes; moments necessário para se calcular a curtose (conforme a fórmula deste pacote); moments: <a href="http://cran.r-project.org/web/packages/moments/moments.pdf">http://cran.r-project.org/web/packages/moments/moments.pdf</a><br>kurtosis(Renda) #Comando kurtosis para se calcular a curtose<br><br>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 href="http://cran.r-project.org/web/packages/e1071/e1071.pdf">http://cran.r-project.org/web/packages/e1071/e1071.pdf</a><br>kurtosis(Renda,type=1) #Comando kurtosis para se calcular o excesso de curtose<br>kurtosis(Renda,type=2) #Comando kurtosis para se calcular o excesso de curtose<br>kurtosis(Renda,type=3) #Comando kurtosis para se calcular o excesso de curtose<br>        <br>#B. ASSIMETRIA<br><br>#Analisa a proximidade ou o afastamento de um conjunto de dados quantitativos em relação à distribuição normal (Gaussiana)<br>#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>#Pearson descreveu duas fórmulas para se calcular o coeficiente de assimetria:<br> #primeiro coeficiente de assimetria de Pearson: AS= média - moda / desvio padrão (conjunto de dados não pode ser plurimodal)<br> #segundo coeficiente de assimetria de Pearson:  AS= Q1 + Q3 - 2Q2 / Q3 - Q1<br>#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> #simétrica: AS=0 (média, mediana e moda coincidem)<br> #assimétrica negativa: AS<0 (cauda da distribuição do lado esquerdo do gráfico); média deslocada para o lado esquerdo<br> #assimétrica positiva: AS>0 (cauda da distribuição do lado direito do gráfico); média deslocada para o lado direito<br>#Comando skewness para se calcular a assimetria<br>#Library(moments): comando library para se carregar pacotes; moments necessário para se calcular a assimetria (conforme a fórmula deste pacote); moments: <a href="http://cran.r-project.org/web/packages/moments/moments.pdf">http://cran.r-project.org/web/packages/moments/moments.pdf</a><br>#Library(e1071): comando library para se carregar pacotes; e1071 necessário para se calcular a assimetria (conforme a fórmula deste pacote); e1071: <a href="http://cran.r-project.org/web/packages/e1071/e1071.pdf">http://cran.r-project.org/web/packages/e1071/e1071.pdf</a><br> #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><br>median(Renda) #Comando median para se calcular a mediana<br>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>asrenda<-(950+1955-(2*1120))/(1955-950) #Fórmula para se calcular o segundo coeficiente de assimetria de Pearson<br>asrenda #Para se visualizar o segundo coeficiente de assimetria de Pearson<br><br>library(moments) #Comando library para se carregar pacotes; moments necessário para se calcular a assimetria (conforme a fórmula deste pacote); moments: <a href="http://cran.r-project.org/web/packages/moments/moments.pdf">http://cran.r-project.org/web/packages/moments/moments.pdf</a><br>skewness(Renda) #Comando skewness para se calcular a assimetria<br><br>library(e1071) #Comando library para se carregar pacotes; e1071 necessário para se calcular a assimetria (conforme a fórmula deste pacote); e1071: <a href="http://cran.r-project.org/web/packages/e1071/e1071.pdf">http://cran.r-project.org/web/packages/e1071/e1071.pdf</a><br>skewness(Renda,type=1) #Comando skewness para se calcular a assimetria<br>skewness(Renda,type=2) #Comando skewness para se calcular a assimetria<br>skewness(Renda,type=3) #Comando skewness para se calcular a assimetria<br><br><br><br></div><div class="gmail_extra"><br><div class="gmail_quote">Em 12 de outubro de 2014 16:15, Daniel Marcelino <span dir="ltr"><<a href="mailto:dmsilva.br@gmail.com" target="_blank">dmsilva.br@gmail.com</a>></span> escreveu:<br><blockquote class="gmail_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>
escolher o tipo de curtose você quer. Veja por exemplo no pacote:<br>
e1071<br>
<br>
x=c(2,3,NA,5,6,8,9)<br>
<br>
require(e1071)<br>
<br>
> kurtosis(x, na.rm = TRUE,type=3)<br>
[1] -1.878889<br>
> kurtosis(x, na.rm = TRUE,type=2)<br>
[1] -1.541333<br>
> kurtosis(x, na.rm = TRUE,type=1)<br>
[1] -1.3856<br>
<br>
require(SciencesPo)<br>
Tipo 2 como padrão.<br>
> detail(x)<br>
  vars obs mean  sd var  se median mad trimmed winsor range min max<br>
1    1   6  5.5 2.7 7.5 1.1    5.5 3.7     5.5    5.5     7   2   9<br>
  skew kurt<br>
1    0 -1.5<br>
<br>
2014-10-12 14:14 GMT-03:00 Andre Oliveira <<a href="mailto:andreolsouza@yahoo.com.br">andreolsouza@yahoo.com.br</a>>:<br>
> Pessoal alguém tem alguma sugestão melhor para estimar assimetria e curtose?<br>
><br>
> Este pacote não bate a definição de  Ferreira, D, F. 2005.<br>
><br>
> curtose=m_4/s^4<br>
> assimetria= m_3/s^3<br>
><br>
> --------------------------------------------------------------------------------------------------------------------------------------------------------<br>
> x=c(2,3,,5,6,8,9)<br>
> library(fBasics)<br>
> basicStats(x)<br>
><br>
><br>
><br>
> André Oliveira Souza.<br>
> Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal<br>
> de Educação, Ciência e Tecnologia do Espirito Santo.  IFES<br>
><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" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código<br>
> mínimo reproduzível.<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" 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" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br>
</blockquote></div><br><br clear="all"><br>-- <br>             \\\ /// <br>           ( @ @ )  <br>    ...o00o.(_).o00o...<br>           Ricardo
</div>