Bab_7_Varianss.R

qf(.95, df1=2, df2=15)

## [1] 3.68232

varians = read.csv("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    NA

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.2323 0.7957
##       14

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.2323 0.7957
##       14

data(varians)

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

leveneTest(as.numeric(varians$Nilai), as.factor(varians$Kelas))

## Levene's Test for Homogeneity of Variance (center = median)
##       Df F value Pr(>F)
## group  2  0.2323 0.7957
##       14

leveneTest(as.numeric(varians$Nilai), as.factor(varians$Kelas), location="mean")

## Levene's Test for Homogeneity of Variance (center = median: "mean")
##       Df F value Pr(>F)
## group  2  0.2323 0.7957
##       14

simpan=read.csv("Data_fiktif_bab7.csv",header = TRUE, sep = ",") #membaca data       
simpan

##    x  y
## 1 30 10
## 2 40 20
## 3 50 30
## 4 60 40
## 5 70 50
## 6 80 60
## 7 90 70
## 8 NA NA

library(doBy)
summaryBy(y ~ x, data = simpan, FUN = function(x)
{ c(ratarata = mean(x), varians = var(x), jumlah = sum(x) ) } ) 

##    x y.ratarata y.varians y.jumlah
## 1 30         10        NA       10
## 2 40         20        NA       20
## 3 50         30        NA       30
## 4 60         40        NA       40
## 5 70         50        NA       50
## 6 80         60        NA       60
## 7 90         70        NA       70
## 8 NA         NA        NA       NA