Statistika Deskriptif

Ahmad Husain A

21 November 2017

Statistika Deskriptif

R menyediakan fungsi untuk mendapatkan gambaran ringkas statistika. Salah satu perintah atau command untuk memperoleh ialah sapply() yang menyediakan summary statistik dari data.

dataku = read.csv("E://BLOG/MATERI/STATISTIKA DESKRIPTIF/data.csv", sep = ";", header = TRUE)
sapply(dataku, mean)

## Warning in mean.default(X[[i]], ...): argument is not numeric or logical:
## returning NA

##             Provinsi Penduduk.2013..Ribu.        PDRB..Miliar. 
##                   NA            7539.9424          229639.9659 
##                  NTP 
##             103.2164

Fungsi yang dapat digunakan pada sapply diantaranya: mean, sd, var, min, max, median, range, dan quantile.

Namun juga ada fungsi yang dalam R yang mampu menampilkan keseluruhan hasil statistika dalam satu kali eksekusi, contohnya:

summary(dataku) #mean, median, quantile, min dan max.

##           Provinsi  Penduduk.2013..Ribu. PDRB..Miliar.    
##  Aceh         : 1   Min.   :  828.3      Min.   :   7725  
##  Bali         : 1   1st Qu.: 2360.4      1st Qu.:  50909  
##  Banten       : 1   Median : 3870.8      Median :  85558  
##  Bengkulu     : 1   Mean   : 7539.9      Mean   : 229640  
##  DI Yogyakarta: 1   3rd Qu.: 7828.7      3rd Qu.: 231683  
##  DKI Jakarta  : 1   Max.   :45340.8      Max.   :1255926  
##  (Other)      :27                                         
##       NTP        
##  Min.   : 88.93  
##  1st Qu.: 99.62  
##  Median :102.90  
##  Mean   :103.22  
##  3rd Qu.:105.90  
##  Max.   :124.70  
## 

Menggunakan Package Hmisc

library(Hmisc)

## Loading required package: lattice

## Loading required package: survival

## Loading required package: Formula

## Loading required package: ggplot2

## 
## Attaching package: 'Hmisc'

## The following objects are masked from 'package:base':
## 
##     format.pval, round.POSIXt, trunc.POSIXt, units

describe(dataku)

## dataku 
## 
##  4  Variables      33  Observations
## ---------------------------------------------------------------------------
## Provinsi 
##        n  missing distinct 
##       33        0       33 
## 
## lowest : Aceh              Bali              Banten            Bengkulu          DI Yogyakarta    
## highest: Sulawesi Tenggara Sulawesi Utara    Sumatera Barat    Sumatera Selatan  Sumatera Utara   
## ---------------------------------------------------------------------------
## Penduduk.2013..Ribu. 
##        n  missing distinct     Info     Mean      Gmd      .05      .10 
##       33        0       33        1     7540     8807     1108     1250 
##      .25      .50      .75      .90      .95 
##     2360     3871     7829    13163    35304 
## 
## lowest :   828.3  1098.0  1114.9  1234.3  1315.1
## highest: 11452.5 13590.3 33264.3 38363.2 45340.8
## ---------------------------------------------------------------------------
## PDRB..Miliar. 
##        n  missing distinct     Info     Mean      Gmd      .05      .10 
##       33        0       33        1   229640   298670    12648    18425 
##      .25      .50      .75      .90      .95 
##    50909    85558   231683   603448  1096637 
## 
## lowest :    7725.42   11752.20   13245.35   16184.01   27388.25
## highest:  522241.43  623749.62 1070177.14 1136326.87 1255925.78
## ---------------------------------------------------------------------------
## NTP 
##        n  missing distinct     Info     Mean      Gmd      .05      .10 
##       33        0       33        1    103.2    7.022    94.73    97.19 
##      .25      .50      .75      .90      .95 
##    99.62   102.90   105.90   109.87   112.79 
## 
## lowest :  88.93  94.23  95.07  97.01  97.93, highest: 109.53 109.95 110.06 116.89 124.70
## ---------------------------------------------------------------------------

Menggunaka Package psych

library(psych)

## 
## Attaching package: 'psych'

## The following object is masked from 'package:Hmisc':
## 
##     describe

## The following objects are masked from 'package:ggplot2':
## 
##     %+%, alpha

describe(dataku)

##                      vars  n      mean        sd   median   trimmed
## Provinsi*               1 33     17.00      9.67    17.00     17.00
## Penduduk.2013..Ribu.    2 33   7539.94  10665.77  3870.80   4770.69
## PDRB..Miliar.           3 33 229639.97 333378.04 85558.31 151220.97
## NTP                     4 33    103.22      6.58   102.90    102.82
##                           mad     min       max      range skew kurtosis
## Provinsi*               11.86    1.00      33.0      32.00 0.00    -1.31
## Penduduk.2013..Ribu.  3048.82  828.30   45340.8   44512.50 2.44     4.96
## PDRB..Miliar.        69123.96 7725.42 1255925.8 1248200.36 1.96     2.74
## NTP                      4.83   88.93     124.7      35.77 0.91     2.02
##                            se
## Provinsi*                1.68
## Penduduk.2013..Ribu.  1856.67
## PDRB..Miliar.        58033.67
## NTP                      1.15

describeBy(dataku)

## Warning in describeBy(dataku): no grouping variable requested

##                      vars  n      mean        sd   median   trimmed
## Provinsi*               1 33     17.00      9.67    17.00     17.00
## Penduduk.2013..Ribu.    2 33   7539.94  10665.77  3870.80   4770.69
## PDRB..Miliar.           3 33 229639.97 333378.04 85558.31 151220.97
## NTP                     4 33    103.22      6.58   102.90    102.82
##                           mad     min       max      range skew kurtosis
## Provinsi*               11.86    1.00      33.0      32.00 0.00    -1.31
## Penduduk.2013..Ribu.  3048.82  828.30   45340.8   44512.50 2.44     4.96
## PDRB..Miliar.        69123.96 7725.42 1255925.8 1248200.36 1.96     2.74
## NTP                      4.83   88.93     124.7      35.77 0.91     2.02
##                            se
## Provinsi*                1.68
## Penduduk.2013..Ribu.  1856.67
## PDRB..Miliar.        58033.67
## NTP                      1.15