# BAB 7 STATISTIK #
# NAMA: ZASKIA LAURA ANANDA
# NIM: 2502056008
# PROGRAM STUDI: ILMU KOMUNIKASI 2025 (A)
###############################################################################
# UJI KESAMAAN VARIANS POPULASI
# uji kesamaan varians populasi dengan uji levene
qf(.95, df1=2, df2=15)
## [1] 3.68232
# penyelesaian dalam R untuk uji kesamaan varians populasi dengan uji levene
varians=read.csv(choose.files("varians.csv"))
varians
## Nilai Kelas
## 1 70 1
## 2 80 1
## 3 87 1
## 4 77 1
## 5 80 1
## 6 80 2
## 7 85 2
## 8 70 2
## 9 77 2
## 10 85 2
## 11 60 2
## 12 80 2
## 13 70 3
## 14 87 3
## 15 90 3
## 16 77 3
## 17 76 3
## 18 87 3
library(car)
## Loading required package: carData
leveneTest(varians$Nilai, varians$Kelas)
## Warning in leveneTest.default(varians$Nilai, varians$Kelas): varians$Kelas
## coerced to factor.
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 2 0.4267 0.6604
## 15
library(Rcmdr)
## Loading required package: splines
## Loading required package: RcmdrMisc
## Loading required package: sandwich
## Loading required package: effects
## Registered S3 method overwritten by 'lme4':
## method from
## na.action.merMod car
## lattice theme set by effectsTheme()
## See ?effectsTheme for details.
## The Commander GUI is launched only in interactive sessions
##
## Attaching package: 'Rcmdr'
## The following object is masked from 'package:base':
##
## errorCondition
leveneTest(varians$Nilai, varians$Kelas)
## Warning in leveneTest.default(varians$Nilai, varians$Kelas): varians$Kelas
## coerced to factor.
## Levene's Test for Homogeneity of Variance (center = median)
## Df F value Pr(>F)
## group 2 0.4267 0.6604
## 15
library(lawstat)
##
## Attaching package: 'lawstat'
## The following object is masked from 'package:car':
##
## levene.test
data(varians)
## Warning in data(varians): data set 'varians' not found
leveneTest(varians[,"Nilai"], varians[,"Kelas"], location="median",correction.method="zero.correction") #Sesuai dengan Minitab
## Warning in leveneTest.default(varians[, "Nilai"], varians[, "Kelas"], location
## = "median", : varians[, "Kelas"] coerced to factor.
## Levene's Test for Homogeneity of Variance (center = median: "median")
## Df F value Pr(>F)
## group 2 0.4267 0.6604
## 15
leveneTest(varians[,"Nilai"], varians[,"Kelas"], location="mean",correction.method="zero.correction") #Sesuai dengan SPSS
## Warning in leveneTest.default(varians[, "Nilai"], varians[, "Kelas"], location
## = "mean", : varians[, "Kelas"] coerced to factor.
## Levene's Test for Homogeneity of Variance (center = median: "mean")
## Df F value Pr(>F)
## group 2 0.4267 0.6604
## 15
# contoh kasus 2, uji kesamaan varians populasi dengan uji levene (contoh perhitungan dan penyelesaian dengan R)
simpan=read.csv(choose.files("varians2.csv"),header=TRUE, sep=",") # membaca data
simpan
## Nilai Kelas
## 1 30 1
## 2 40 1
## 3 50 1
## 4 60 1
## 5 70 1
## 6 80 1
## 7 90 2
## 8 10 2
## 9 20 2
## 10 30 2
## 11 40 2
## 12 50 2
## 13 60 2
## 14 70 2
library(doBy)
summaryBy(Nilai ~ Kelas, data = simpan, FUN = function(x){ c(ratarata = mean(x), varians = var(x), jumlah=sum(x) )})
## Kelas Nilai.ratarata Nilai.varians Nilai.jumlah
## 1 1 55.00 350.0 330
## 2 2 46.25 712.5 370
qf(0.95, df1=1, df2=12)
## [1] 4.747225
library(lawstat)
data(simpan)
## Warning in data(simpan): data set 'simpan' not found
leveneTest(simpan[,"Nilai"], simpan[,"Kelas"], location="mean", correction.method="zero.correction")
## Warning in leveneTest.default(simpan[, "Nilai"], simpan[, "Kelas"], location =
## "mean", : simpan[, "Kelas"] coerced to factor.
## Levene's Test for Homogeneity of Variance (center = median: "mean")
## Df F value Pr(>F)
## group 1 0.8991 0.3617
## 12