statistik-bab-8.R

data=read.csv("data_berpasangan.csv")
data

##    X   Y
## 1 78 100
## 2 75  95
## 3 67  90
## 4 77  90
## 5 70  90
## 6 72  90
## 7 78  89
## 8 74  90
## 9 77 100

t.test(data$Y, data$X, paired=TRUE)

## 
##  Paired t-test
## 
## data:  data$Y and data$X
## t = 12.769, df = 8, p-value = 1.333e-06
## alternative hypothesis: true mean difference is not equal to 0
## 95 percent confidence interval:
##  15.11355 21.77534
## sample estimates:
## mean difference 
##        18.44444

data=read.csv("data_berpasangan.csv")
data

##    X   Y
## 1 78 100
## 2 75  95
## 3 67  90
## 4 77  90
## 5 70  90
## 6 72  90
## 7 78  89
## 8 74  90
## 9 77 100

selisih = data$Y - data$X
selisih

## [1] 22 20 23 13 20 18 11 16 23

qqnorm(selisih)
qqline(selisih)

library(nortest)
lillie.test(selisih)

## 
##  Lilliefors (Kolmogorov-Smirnov) normality test
## 
## data:  selisih
## D = 0.19575, p-value = 0.41

library(tseries)

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

jarque.bera.test(selisih)

## 
##  Jarque Bera Test
## 
## data:  selisih
## X-squared = 0.8536, df = 2, p-value = 0.6526