bab7_levene2.R

Simpan=read.csv("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     1
## 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             60      466.6667          420
## 2     2             40      466.6667          280

qf(0.95, df1=1, df2=12)

## [1] 4.747225

library(lawstat)
data(Simpan)

## Warning in data(Simpan): data set 'Simpan' not found

levene.test(Simpan[,"Nilai"], Simpan[,"Kelas"], location="mean", correction.method="zero.correction")

## 
##  Classical Levene's test based on the absolute deviations from the mean
##  ( zero.correction not applied because the location is not set to median
##  )
## 
## data:  Simpan[, "Nilai"]
## Test Statistic = 1.4336e-32, p-value = 1