2409106083_MuhammadRayhanSetiawan_TugasBab6.R

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

summary(x)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##      40      60      70      70      80     100

library(nortest) #install package 'nortest' terlebih dahulu
lillie.test(x)

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

nilai <- c(40, 50, 50, 60, 60, 60, 70, 70,
           70, 70, 80, 80, 80, 90, 90, 100)

cat("STATISTIK DESKRIPTIF\\n")

## STATISTIK DESKRIPTIF\n

cat("Jumlah data :", length(nilai), "\\n")

## Jumlah data : 16 \n

cat("Rata-rata   :", mean(nilai), "\\n")

## Rata-rata   : 70 \n

cat("Median      :", median(nilai), "\\n")

## Median      : 70 \n

cat("Standar deviasi:", sd(nilai), "\\n\\n")

## Standar deviasi: 16.32993 \n\n

alpha <- 0.05

hasil_lilliefors <- lillie.test(nilai)

cat("HASIL UJI LILLIEFORS\\n")

## HASIL UJI LILLIEFORS\n

print(hasil_lilliefors)

## 
##  Lilliefors (Kolmogorov-Smirnov) normality test
## 
## data:  nilai
## D = 0.125, p-value = 0.7235

if (hasil_lilliefors$p.value > alpha) {
  cat("Keputusan: Gagal menolak H0. Tidak terdapat bukti yang cukup")
  cat(" bahwa data tidak berdistribusi normal.\\n\\n")
} else {
  cat("Keputusan: Tolak H0. Terdapat bukti bahwa data tidak")
  cat(" berdistribusi normal.\\n\\n")
}

## Keputusan: Gagal menolak H0. Tidak terdapat bukti yang cukup bahwa data tidak berdistribusi normal.\n\n

library(tseries) 

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

hasil_jarque_bera <- jarque.bera.test(nilai)
cat("HASIL UJI JARQUE-BERA\n")

## HASIL UJI JARQUE-BERA

print(hasil_jarque_bera)

## 
##  Jarque Bera Test
## 
## data:  nilai
## X-squared = 0.30827, df = 2, p-value = 0.8572

if (hasil_jarque_bera$p.value > alpha) {
  cat("Keputusan: Gagal menolak H0. Tidak terdapat bukti yang cukup")
  cat(" bahwa data tidak berdistribusi normal.\\n\\n")
} else {
  cat("Keputusan: Tolak H0. Terdapat bukti bahwa data tidak")
  cat(" berdistribusi normal.\\n\\n")
}

## Keputusan: Gagal menolak H0. Tidak terdapat bukti yang cukup bahwa data tidak berdistribusi normal.\n\n

qqnorm(nilai,
       main = "Q-Q Plot Nilai Ujian Matematika",
       xlab = "Kuantil Teoretis",
       ylab = "Kuantil Sampel")
qqline(nilai, col = "red", lwd = 2)

cat("KESIMPULAN\\n")

## KESIMPULAN\n

if (hasil_lilliefors$p.value > alpha &&
    hasil_jarque_bera$p.value > alpha) {
  cat("Berdasarkan uji Lilliefors dan Jarque-Bera pada taraf")
  cat(" signifikansi 5%, kedua p-value lebih besar dari 0.05.")
  cat(" Dengan demikian, H0 gagal ditolak dan tidak terdapat")
  cat(" bukti yang cukup untuk menyatakan bahwa data tidak")
  cat(" berdistribusi normal.\\n")
} else {
  cat("Setidaknya salah satu pengujian menghasilkan p-value")
  cat(" kurang dari atau sama dengan 0.05. Periksa hasil")
  cat(" pengujian dan Q-Q Plot untuk menilai normalitas data.\\n")
}

## Berdasarkan uji Lilliefors dan Jarque-Bera pada taraf signifikansi 5%, kedua p-value lebih besar dari 0.05. Dengan demikian, H0 gagal ditolak dan tidak terdapat bukti yang cukup untuk menyatakan bahwa data tidak berdistribusi normal.\n