Bevezetés az R-be 2.0 - Feladatgyűjtemény

Levene-próba

1. feladat. Energiabevitel variabilitásának összehasonlítása Levene-próbával.
Egy 22 elemű mintában sovány (lean) és túlsúlyos (obese) nők napi energiabevitelét vizsgáljuk. Vizsgáljuk meg a sovány és túlsúlyos nők csoportjaiban a napi energiabevitel változó varianciáinak egyezését Levene-próbával! Az adatok az ISwR csomag energy adattáblája tartalmazza.

Az állítás elenőrzésére Levene-próbát hajtunk végre, a populációbeli varianciák egyezését vizsgáljuk:

• \(H_0:\sigma^2_1=\sigma^2_2\)
• \(H_1:\sigma^2_1 \neq \sigma^2_2\)

Adatok beolvasása

data(energy, package = "ISwR")
str(energy)

'data.frame':    22 obs. of  2 variables:
 $ expend : num  9.21 7.53 7.48 8.08 8.09 ...
 $ stature: Factor w/ 2 levels "lean","obese": 2 1 1 1 1 1 1 1 1 1 ...

energy

   expend stature
1    9.21   obese
2    7.53    lean
3    7.48    lean
4    8.08    lean
5    8.09    lean
6   10.15    lean
7    8.40    lean
8   10.88    lean
9    6.13    lean
10   7.90    lean
11  11.51   obese
12  12.79   obese
13   7.05    lean
14  11.85   obese
15   9.97   obese
16   7.48    lean
17   8.79   obese
18   9.69   obese
19   9.68   obese
20   7.58    lean
21   9.19   obese
22   8.11    lean

Levene-próba végrehajtása

library(car)
leveneTest(expend ~ stature, data = energy, center = mean)

Levene's Test for Homogeneity of Variance (center = mean)
      Df F value Pr(>F)
group  1   1.002 0.3288
      20

A próba nem szignifikáns, azaz nem találtunk elég bizonyítékot arra, hogy az energiabevitel varianciája a két csoportban eltér (\(F(1,20)=1,002; p=0,329\)).

2. feladat. Gyártott gépek átmérőjének pontossága.
Egy üzemben 10 szakaszban gyártanak gépeket, melyeknek átmérő adatait gyűjtöttük össze: http://users.atw.hu/rdata/ber2/gep.txt. Vizsgáljuk meg Levene-próbával, hogy az egyes szakaszokban gyártott gépek átmérőinek a pontossága eltér-e egymástól!

Az állítás elenőrzésére Levene-próbát hajtunk végre, a populációbeli varianciák egyezését vizsgáljuk:

• $$H_0:\sigma^2_1=\sigma^2_2= \ldots =\sigma^2_{10}$$
• $H_1:$ Legalább két populációbeli variancia eltér egymástól

Adatok beolvasása

d <- read.table("http://users.atw.hu/rdata/ber2/gep.txt", sep=" ", dec=".", header=T)
str(d)

'data.frame':    100 obs. of  2 variables:
 $ atmero : num  1.006 0.996 0.998 1 0.992 ...
 $ szakasz: num  1 1 1 1 1 1 1 1 1 1 ...

d$szakasz <- factor(d$szakasz)
str(d)

'data.frame':    100 obs. of  2 variables:
 $ atmero : num  1.006 0.996 0.998 1 0.992 ...
 $ szakasz: Factor w/ 10 levels "1","2","3","4",..: 1 1 1 1 1 1 1 1 1 1 ...

Levene-próba végrehajtása

library(car)
leveneTest(atmero ~ szakasz, data = d, center = median)

Levene's Test for Homogeneity of Variance (center = median)
      Df F value  Pr(>F)  
group  9  1.7059 0.09908 .
      90                  
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

A próba nem szignifikáns, azaz nincs elegendő bizonyíték arra, hogy az egyes szakaszokban gyártott gépek átmérőjének varianciája eltér egymástól (\(F(9,90)=1,706; p=0,099\)).