ทำความเข้าใจผลลัพธ์จาก logistic regression

Generalized Linear Models ใน R

Richard Erickson

Instructor

การสื่อสารผลลัพธ์จาก logistic regression?

  • Linear regression เข้าใจง่าย:

    • บวก intercept
    • คูณ slope
  • Poisson regression:

    • ต้องแปลงด้วย exponential
    • หลังแปลงแล้วคล้าย linear regression
  • Logistic regression???
Generalized Linear Models ใน R

Odds-ratios

  • ไม่ตรงไปตรงมาเหมือน Poisson exponential
  • ใช้เปรียบเทียบโอกาสสัมพัทธ์ของสองเหตุการณ์
Generalized Linear Models ใน R

ตัวอย่าง odds-ratios

  • เหรียญที่ไม่ยุติธรรม:
    • เปรียบเทียบหัวกับก้อย
    • หัวออก 3 ครั้งต่อก้อย 1 ครั้ง
    • อัตราต่อรอง 3 ต่อ 1
    • Odds-ratio คือ 3.0
  • มักใช้ในกีฬา/การพนัน
  • และการศึกษาทางการแพทย์
Generalized Linear Models ใน R

การหา odds-ratio จาก logistic regression

Log-odds ("logit"):

$\phi(x) = \text{ln}(\frac{p(x)}{1 - p(x)}) = \beta_0 + \beta_1 x$

Odds นำ exponential ($e^{x}$):

$\frac{p(x)}{1 - p(x)} = e^{\beta_0 + \beta_1x}$

Generalized Linear Models ใน R

Odds-ratio สำหรับตัวแปรต่อเนื่อง

Odds-ratio (OR) สำหรับตัวแปรต่อเนื่อง:

$\text{OR} = \frac{e^{\beta_0 + \beta_1(x +1)}}{e^{\beta_0 + \beta_1x}} = e^{\beta_1}$

Generalized Linear Models ใน R

การตีความ

ค่า OR:

  • OR = 1: coefficient ไม่มีผล
  • OR < 1: coefficient ลดโอกาส
  • OR > 1: coefficient เพิ่มโอกาส
Generalized Linear Models ใน R

ตัวอย่างมะเร็ง

ชายไม่สูบบุหรี่เทียบกับสูบบุหรี่ (Pesch et al. 2012)

  • OR: 103.5 (95% CI 74.8-143.2)
  • $>$ โอกาสเป็นมะเร็งของชายที่สูบบุหรี่สูงกว่า 100 ต่อ 1!
  • วรรณกรรมทางการแพทย์มักรายงาน 95% confidence interval แทน p-value
  • แนวโน้มที่กว้างขึ้นในการหลีกเลี่ยง p-value
Generalized Linear Models ใน R

ดึงค่าจาก GLM

glm_out <- glm(y ~ x, family = 'binomial')

coef(glm_out)

exp(coef(glm_out))

confint(glm_out)

exp(confint(glm_out))
Generalized Linear Models ใน R

Tidyverse

library(broom)

glm_out <- glm(y ~ x, family = 'binomial')

tidy(glm_out, exponentiate = TRUE, conf.int= TRUE)
Generalized Linear Models ใน R

มาฝึกกันเถอะ!

Generalized Linear Models ใน R

Preparing Video For Download...