UJI NORMALITAS POPULASI

normalitas = read.csv("normalitas.csv")
x = normalitas$Data


library(nortest)
lillie.test(x)

## 
##  Lilliefors (Kolmogorov-Smirnov) normality test
## 
## data:  x
## D = 0.13341, p-value = 0.6708

library(e1071)
a = moment(x, order = 3, center = TRUE)
b = moment(x, order = 4, center = TRUE)
a

## [1] 421.9259

b

## [1] 143706.1

c = a^2
c

## [1] 178021.5

kurtosis = b / c
kurtosis

## [1] 0.8072401

library(tseries)

## Registered S3 method overwritten by 'quantmod':
##   method            from
##   as.zoo.data.frame zoo

jarque.bera.test(x)

## 
##  Jarque Bera Test
## 
## data:  x
## X-squared = 0.16965, df = 2, p-value = 0.9187

qqnorm(x)
qqline(x)