datanormalitas.R

normalitas=read.csv("datanormalitas.csv")
X=normalitas$Data

library(e1071)

a=moment(X, order=2, center=TRUE)
b=moment(X, order=4, center=TRUE)
a

## [1] 266.6667

b

## [1] 154666.7

c=a^2
c

## [1] 71111.11

kurtosis=b/c
kurtosis

## [1] 2.175

tseries::jarque.bera.test(X)

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

## 
##  Jarque Bera Test
## 
## data:  X
## X-squared = 0.42539, df = 2, p-value = 0.8084