[R-br] Estimar assimetria e curtose
Ricardo Fernando Zampieri
rfzampieri em gmail.com
Domingo Outubro 12 19:58:04 BRT 2014
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:
#A. CURTOSE
#Analisa a proximidade ou o afastamento de um conjunto de dados
quantitativos em relação à distribuição normal (Gaussiana)
#Caracteriza o achatamento ou o alongamento da curva da função de
distribuição de probabilidade
#Fórmula convencional para se calcular a curtose: C= Q3 - Q1 / 2(P90 - P10)
#mede o grau de achatamento ou alongamento da curva
#quanto > valor da curtose > achatamento da curva
#Formas da curva da função de distribuição de probabilidade:
#mesocúrtica: mesma forma que a distribuição normal (C= 0,263)
#platicúrtica: maior achatamento que a distribuição normal (C> 0,263)
#leptocúrtica: mais alta e concentrada que a distribuição normal (C< 0,263)
#Para se calcular a curtose: pacotes moments e e1071
#Comando kurtosis para se calcular a curtose (moments) ou o excesso de
curtose (e1071)
#Library(moments): comando library para se carregar pacotes; moments
necessário para se calcular a curtose (conforme a fórmula deste pacote);
moments: http://cran.r-project.org/web/packages/moments/moments.pdf
#pode-se obter: C=3 (mesocúrtica) C>3 (platicúrtica) C<3 (leptocúrtica)
#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: http://cran.r-project.org/web/packages/e1071/e1071.pdf
#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)
#kurtosis(Renda,type=1) retorna o mesmo valor de moments-3
#no pacote e1071, a distribuição normal padrão apresenta referencialmente
o valor de curtose=3 (em vez de 0,263)
#calcula-se então o excesso de curtose por meio do comando kurtosis:
#ExCur=C-3 (o programa calcula o valor de C, conforme a fórmula deste
pacote, e o subtrai de 3)
#pode-se obter: ExCur=0 (mesocúrtica) ExCur>0 (platicúrtica) ExCur<0
(leptocúrtica)
#Em distribuições normais, independentemente do número de dados amostrais,
o excesso de curtose terá valor próximo de 0
median(Renda) #Comando median para se calcular a mediana
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)
k1<-(1955-950)/(2*(2155-920)) #Fórmula convencional para se calcular a
curtose
k1 #Para se visualizar curtose
library(moments) #Comando library para se carregar pacotes; moments
necessário para se calcular a curtose (conforme a fórmula deste pacote);
moments: http://cran.r-project.org/web/packages/moments/moments.pdf
kurtosis(Renda) #Comando kurtosis para se calcular a curtose
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: http://cran.r-project.org/web/packages/e1071/e1071.pdf
kurtosis(Renda,type=1) #Comando kurtosis para se calcular o excesso de
curtose
kurtosis(Renda,type=2) #Comando kurtosis para se calcular o excesso de
curtose
kurtosis(Renda,type=3) #Comando kurtosis para se calcular o excesso de
curtose
#B. ASSIMETRIA
#Analisa a proximidade ou o afastamento de um conjunto de dados
quantitativos em relação à distribuição normal (Gaussiana)
#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
#Pearson descreveu duas fórmulas para se calcular o coeficiente de
assimetria:
#primeiro coeficiente de assimetria de Pearson: AS= média - moda / desvio
padrão (conjunto de dados não pode ser plurimodal)
#segundo coeficiente de assimetria de Pearson: AS= Q1 + Q3 - 2Q2 / Q3 - Q1
#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
#simétrica: AS=0 (média, mediana e moda coincidem)
#assimétrica negativa: AS<0 (cauda da distribuição do lado esquerdo do
gráfico); média deslocada para o lado esquerdo
#assimétrica positiva: AS>0 (cauda da distribuição do lado direito do
gráfico); média deslocada para o lado direito
#Comando skewness para se calcular a assimetria
#Library(moments): comando library para se carregar pacotes; moments
necessário para se calcular a assimetria (conforme a fórmula deste pacote);
moments: http://cran.r-project.org/web/packages/moments/moments.pdf
#Library(e1071): comando library para se carregar pacotes; e1071 necessário
para se calcular a assimetria (conforme a fórmula deste pacote); e1071:
http://cran.r-project.org/web/packages/e1071/e1071.pdf
#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)
median(Renda) #Comando median para se calcular a mediana
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)
asrenda<-(950+1955-(2*1120))/(1955-950) #Fórmula para se calcular o segundo
coeficiente de assimetria de Pearson
asrenda #Para se visualizar o segundo coeficiente de assimetria de Pearson
library(moments) #Comando library para se carregar pacotes; moments
necessário para se calcular a assimetria (conforme a fórmula deste pacote);
moments: http://cran.r-project.org/web/packages/moments/moments.pdf
skewness(Renda) #Comando skewness para se calcular a assimetria
library(e1071) #Comando library para se carregar pacotes; e1071 necessário
para se calcular a assimetria (conforme a fórmula deste pacote); e1071:
http://cran.r-project.org/web/packages/e1071/e1071.pdf
skewness(Renda,type=1) #Comando skewness para se calcular a assimetria
skewness(Renda,type=2) #Comando skewness para se calcular a assimetria
skewness(Renda,type=3) #Comando skewness para se calcular a assimetria
Em 12 de outubro de 2014 16:15, Daniel Marcelino <dmsilva.br em gmail.com>
escreveu:
> André, eu acho que há muitas opções disponíveis em que permite você
> escolher o tipo de curtose você quer. Veja por exemplo no pacote:
> e1071
>
> x=c(2,3,NA,5,6,8,9)
>
> require(e1071)
>
> > kurtosis(x, na.rm = TRUE,type=3)
> [1] -1.878889
> > kurtosis(x, na.rm = TRUE,type=2)
> [1] -1.541333
> > kurtosis(x, na.rm = TRUE,type=1)
> [1] -1.3856
>
> require(SciencesPo)
> Tipo 2 como padrão.
> > detail(x)
> vars obs mean sd var se median mad trimmed winsor range min max
> 1 1 6 5.5 2.7 7.5 1.1 5.5 3.7 5.5 5.5 7 2 9
> skew kurt
> 1 0 -1.5
>
> 2014-10-12 14:14 GMT-03:00 Andre Oliveira <andreolsouza em yahoo.com.br>:
> > Pessoal alguém tem alguma sugestão melhor para estimar assimetria e
> curtose?
> >
> > Este pacote não bate a definição de Ferreira, D, F. 2005.
> >
> > curtose=m_4/s^4
> > assimetria= m_3/s^3
> >
> >
> --------------------------------------------------------------------------------------------------------------------------------------------------------
> > x=c(2,3,,5,6,8,9)
> > library(fBasics)
> > basicStats(x)
> >
> >
> >
> > 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
> >
> >
> > _______________________________________________
> > R-br mailing list
> > R-br em listas.c3sl.ufpr.br
> > 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.
> _______________________________________________
> R-br mailing list
> R-br em listas.c3sl.ufpr.br
> 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.
>
--
\\\ ///
( @ @ )
...o00o.(_).o00o...
Ricardo
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20141012/4c5afc71/attachment.html>
Mais detalhes sobre a lista de discussão R-br