PADK Prak 6

Contoh Syntax

# Membuat tabel kontingensi hasil rapid test
rapid_test <- matrix(c(10, 45, 50, 95), nrow = 2, byrow = TRUE,
                     dimnames = list("Prediksi" = c("Positif", "Negatif"),
                                    "Aktual" = c("Positif", "Negatif")))

# Menampilkan tabel kontingensi
print(rapid_test)

##          Aktual
## Prediksi  Positif Negatif
##   Positif      10      45
##   Negatif      50      95

# Menghitung metrik-metrik performa
TP <- rapid_test["Positif", "Positif"] # True Positive
FP <- rapid_test["Positif", "Negatif"] # False Positive
FN <- rapid_test["Negatif", "Positif"] # False Negative
TN <- rapid_test["Negatif", "Negatif"] # True Negative

# Akurasi
akurasi <- (TP + TN) / sum(rapid_test)
print(paste("Akurasi:", akurasi))

## [1] "Akurasi: 0.525"

# Sensitivitas (Recall)
sensitivitas <- TP / (TP + FN)
print(paste("Sensitivitas:", sensitivitas))

## [1] "Sensitivitas: 0.166666666666667"

# Presisi (Precision)
presisi <- TP / (TP + FP)
print(paste("Presisi:", presisi))

## [1] "Presisi: 0.181818181818182"

Latihan Soal

##1

# Membuat tabel kontingensi hasil rapid test
rapid_test1 <- matrix(c(50, 30, 20, 200), nrow = 2, byrow = TRUE,
                     dimnames = list("Prediksi" = c("Positif", "Negatif"),
                                    "Aktual" = c("Positif", "Negatif")))

# Menampilkan tabel kontingensi
print(rapid_test1)

##          Aktual
## Prediksi  Positif Negatif
##   Positif      50      30
##   Negatif      20     200

# Menghitung metrik-metrik performa
TP1 <- rapid_test1["Positif", "Positif"] # True Positive
FP1 <- rapid_test1["Positif", "Negatif"] # False Positive
FN1 <- rapid_test1["Negatif", "Positif"] # False Negative
TN1 <- rapid_test1["Negatif", "Negatif"] # True Negative

# Akurasi
akurasi1 <- (TP1 + TN1) / sum(rapid_test1)
print(paste("Akurasi:", akurasi1))

## [1] "Akurasi: 0.833333333333333"

# Sensitivitas (Recall)
sensitivitas1 <- TP1 / (TP1 + FN1)
print(paste("Sensitivitas:", sensitivitas1))

## [1] "Sensitivitas: 0.714285714285714"

# Spesifisitas
spesifisitas1 <- TN1 / (TN1 + FP1)
print(paste("Spesifisitas:", spesifisitas1))

## [1] "Spesifisitas: 0.869565217391304"

# Presisi (Precision)
presisi1 <- TP1 / (TP1 + FP1)
print(paste("Presisi:", presisi1))

## [1] "Presisi: 0.625"

2

# Membuat tabel kontingensi hasil rapid test
rapid_test2 <- matrix(c(10, 190, 50, 150), nrow = 2, byrow = TRUE,
                     dimnames = list("Prediksi" = c("Positif", "Negatif"),
                                    "Aktual" = c("Positif", "Negatif")))

# Menampilkan tabel kontingensi
print(rapid_test2)

##          Aktual
## Prediksi  Positif Negatif
##   Positif      10     190
##   Negatif      50     150

# Menghitung metrik-metrik performa
TP2 <- rapid_test2["Positif", "Positif"] # True Positive
FP2 <- rapid_test2["Positif", "Negatif"] # False Positive
FN2 <- rapid_test2["Negatif", "Positif"] # False Negative
TN2 <- rapid_test2["Negatif", "Negatif"] # True Negative

# Akurasi
akurasi2 <- (TP2 + TN2) / sum(rapid_test2)
print(paste("Akurasi:", akurasi2))

## [1] "Akurasi: 0.4"

# Sensitivitas (Recall)
sensitivitas2 <- TP2 / (TP2 + FN2)
print(paste("Sensitivitas:", sensitivitas2))

## [1] "Sensitivitas: 0.166666666666667"

# Spesifisitas
spesifisitas2 <- TN2 / (TN2 + FP2)
print(paste("Spesifisitas:", spesifisitas2))

## [1] "Spesifisitas: 0.441176470588235"

# Presisi (Precision)
presisi2 <- TP2 / (TP2 + FP2)
print(paste("Presisi:", presisi2))

## [1] "Presisi: 0.05"

3

# Membuat tabel kontingensi hasil rapid test
rapid_test3 <- matrix(c(30, 15, 10, 195), nrow = 2, byrow = TRUE,
                     dimnames = list("Prediksi" = c("Positif", "Negatif"),
                                    "Aktual" = c("Positif", "Negatif")))

# Menampilkan tabel kontingensi
print(rapid_test3)

##          Aktual
## Prediksi  Positif Negatif
##   Positif      30      15
##   Negatif      10     195

# Menghitung metrik-metrik performa
TP3 <- rapid_test3["Positif", "Positif"] # True Positive
FP3 <- rapid_test3["Positif", "Negatif"] # False Positive
FN3 <- rapid_test3["Negatif", "Positif"] # False Negative
TN3 <- rapid_test3["Negatif", "Negatif"] # True Negative

# Akurasi
akurasi3 <- (TP3 + TN3) / sum(rapid_test3)
print(paste("Akurasi:", akurasi3))

## [1] "Akurasi: 0.9"

# Sensitivitas (Recall)
sensitivitas3 <- TP3 / (TP3 + FN3)
print(paste("Sensitivitas:", sensitivitas3))

## [1] "Sensitivitas: 0.75"

# Spesifisitas
spesifisitas3 <- TN3 / (TN3 + FP3)
print(paste("Spesifisitas:", spesifisitas3))

## [1] "Spesifisitas: 0.928571428571429"

# Presisi (Precision)
presisi3 <- TP3 / (TP3 + FP3)
print(paste("Presisi:", presisi3))

## [1] "Presisi: 0.666666666666667"

4

# Membuat tabel kontingensi hasil rapid test
rapid_test4 <- matrix(c(40, 60, 20, 380), nrow = 2, byrow = TRUE,
                     dimnames = list("Prediksi" = c("Positif", "Negatif"),
                                    "Aktual" = c("Positif", "Negatif")))

# Menampilkan tabel kontingensi
print(rapid_test4)

##          Aktual
## Prediksi  Positif Negatif
##   Positif      40      60
##   Negatif      20     380

# Menghitung metrik-metrik performa
TP4 <- rapid_test4["Positif", "Positif"] # True Positive
FP4 <- rapid_test4["Positif", "Negatif"] # False Positive
FN4 <- rapid_test4["Negatif", "Positif"] # False Negative
TN4 <- rapid_test4["Negatif", "Negatif"] # True Negative

# Akurasi
akurasi4 <- (TP4 + TN4) / sum(rapid_test4)
print(paste("Akurasi:", akurasi4))

## [1] "Akurasi: 0.84"

# Sensitivitas (Recall)
sensitivitas4 <- TP4 / (TP4 + FN4)
print(paste("Sensitivitas:", sensitivitas4))

## [1] "Sensitivitas: 0.666666666666667"

# Spesifisitas
spesifisitas4 <- TN4 / (TN4 + FP4)
print(paste("Spesifisitas:", spesifisitas4))

## [1] "Spesifisitas: 0.863636363636364"

# Presisi (Precision)
presisi4 <- TP4 / (TP4 + FP4)
print(paste("Presisi:", presisi4))

## [1] "Presisi: 0.4"

5

# Membuat tabel kontingensi hasil rapid test
rapid_test5 <- matrix(c(30, 270, 70, 230), nrow = 2, byrow = TRUE,
                     dimnames = list("Prediksi" = c("Positif", "Negatif"),
                                    "Aktual" = c("Positif", "Negatif")))

# Menampilkan tabel kontingensi
print(rapid_test5)

##          Aktual
## Prediksi  Positif Negatif
##   Positif      30     270
##   Negatif      70     230

# Menghitung metrik-metrik performa
TP5 <- rapid_test5["Positif", "Positif"] # True Positive
FP5 <- rapid_test5["Positif", "Negatif"] # False Positive
FN5 <- rapid_test5["Negatif", "Positif"] # False Negative
TN5 <- rapid_test5["Negatif", "Negatif"] # True Negative

# Akurasi
akurasi5 <- (TP5 + TN5) / sum(rapid_test5)
print(paste("Akurasi:", akurasi5))

## [1] "Akurasi: 0.433333333333333"

# Sensitivitas (Recall)
sensitivitas5 <- TP5 / (TP5 + FN5)
print(paste("Sensitivitas:", sensitivitas5))

## [1] "Sensitivitas: 0.3"

# Spesifisitas
spesifisitas5 <- TN5 / (TN5 + FP5)
print(paste("Spesifisitas:", spesifisitas5))

## [1] "Spesifisitas: 0.46"

# Presisi (Precision)
presisi5 <- TP5 / (TP5 + FP5)
print(paste("Presisi:", presisi5))

## [1] "Presisi: 0.1"