Tugas 5.9

Tugas 5.9: Ukuran Statistik dalam Epidemiologi

# ============================
# Soal 1: Incidence Rate
# ============================
kasus1 <- 8
person_time <- 1050
IR <- kasus1 / person_time
IR_per1000 <- IR * 1000
cat("Soal 1: Incidence Rate\n")

## Soal 1: Incidence Rate

cat("IR =", round(IR,5), "per person-year\n")

## IR = 0.00762 per person-year

cat("IR =", round(IR_per1000,2), "per 1.000 person-years\n")

## IR = 7.62 per 1.000 person-years

cat("Kesimpulan: Ada sekitar", round(IR_per1000,2), 
    "kasus baru per 1.000 orang per tahun.\n\n")

## Kesimpulan: Ada sekitar 7.62 kasus baru per 1.000 orang per tahun.

# ============================
# Soal 2: Cumulative Incidence
# ============================
kasus2 <- 25
pop_awal <- 500
CI <- kasus2 / pop_awal
CI_percent <- CI * 100
cat("Soal 2: Cumulative Incidence\n")

## Soal 2: Cumulative Incidence

cat("CI =", round(CI,3), "(", round(CI_percent,1), "% )\n")

## CI = 0.05 ( 5 % )

cat("Kesimpulan: Risiko hipertensi dalam 2 tahun adalah",
    round(CI_percent,1), "% dari populasi.\n\n")

## Kesimpulan: Risiko hipertensi dalam 2 tahun adalah 5 % dari populasi.

# ============================
# Soal 3: Prevalensi
# ============================
pop3 <- 20000
kasus3 <- 400
Prev <- kasus3 / pop3
Prev_percent <- Prev * 100
cat("Soal 3: Prevalensi\n")

## Soal 3: Prevalensi

cat("Prevalensi =", round(Prev,3), "(", Prev_percent, "% )\n")

## Prevalensi = 0.02 ( 2 % )

cat("Kesimpulan: Pada saat survei,", Prev_percent, 
    "% penduduk menderita diabetes.\n\n")

## Kesimpulan: Pada saat survei, 2 % penduduk menderita diabetes.

# ============================
# Soal 4: RR & AR
# ============================
exp_total <- 200; exp_cases <- 40
unexp_total <- 300; unexp_cases <- 15
CI_exp <- exp_cases / exp_total
CI_unexp <- unexp_cases / unexp_total
RR <- CI_exp / CI_unexp
AR <- CI_exp - CI_unexp
cat("Soal 4: Relative Risk & Attributable Risk\n")

## Soal 4: Relative Risk & Attributable Risk

cat("CI perokok =", round(CI_exp,3), "(", round(CI_exp*100,1), "% )\n")

## CI perokok = 0.2 ( 20 % )

cat("CI bukan perokok =", round(CI_unexp,3), "(", round(CI_unexp*100,1), "% )\n")

## CI bukan perokok = 0.05 ( 5 % )

cat("RR =", round(RR,2), "\n")

## RR = 4

cat("AR =", round(AR,3), "(", round(AR*100,1), "% )\n")

## AR = 0.15 ( 15 % )

cat("Kesimpulan: Perokok memiliki risiko", round(RR,2), 
    "kali lebih tinggi. Tambahan risiko akibat merokok adalah",
    round(AR*100,1), "%.\n\n")

## Kesimpulan: Perokok memiliki risiko 4 kali lebih tinggi. Tambahan risiko akibat merokok adalah 15 %.

# ============================
# Soal 5: Odds Ratio
# ============================
a <- 45; b <- 30; c <- 20; d <- 55
OR <- (a*d) / (b*c)
cat("Soal 5: Odds Ratio\n")

## Soal 5: Odds Ratio

cat("OR =", round(OR,2), "\n")

## OR = 4.12

cat("Kesimpulan: Odds penyakit pada yang terpapar",
    "adalah", round(OR,2), "kali lebih tinggi.\n\n")

## Kesimpulan: Odds penyakit pada yang terpapar adalah 4.12 kali lebih tinggi.

# ============================
# Soal 6: Case Fatality Rate
# ============================
kasus6 <- 250
death <- 10
CFR <- (death / kasus6) * 100
cat("Soal 6: Case Fatality Rate\n")

## Soal 6: Case Fatality Rate

cat("CFR =", CFR, "%\n")

## CFR = 4 %

cat("Kesimpulan: Tingkat keparahan penyakit ditunjukkan dengan",
    CFR, "% kasus berakhir kematian.\n")

## Kesimpulan: Tingkat keparahan penyakit ditunjukkan dengan 4 % kasus berakhir kematian.