ทดสอบ differential expression สำหรับ 3 กลุ่ม

เมื่อกำหนด design matrix และ contrasts matrix เรียบร้อยแล้ว ก็สามารถทดสอบ differential expression ได้เลย

แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร

การวิเคราะห์ Differential Expression ด้วย limma ใน R

คำแนะนำการฝึกหัด

อ็อบเจกต์ ExpressionSet eset ที่มีข้อมูล hypoxia, design matrix (design) และ contrasts matrix (cm) ถูกโหลดไว้ใน workspace แล้ว

ฟิตค่าสัมประสิทธิ์ของโมเดลด้วย lmFit
ฟิต contrasts ด้วย contrasts.fit
คำนวณ t-statistics ด้วย eBayes
สรุปผลลัพธ์ด้วย decideTests

แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ

ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์

# Load package
library(limma)

# Fit the model
fit <- ___(eset, ___)

# Fit the contrasts
fit2 <- ___(fit, contrasts = ___)

# Calculate the t-statistics for the contrasts
fit2 <- ___(fit2)

# Summarize results
results <- ___(fit2)
summary(results)

แก้ไขและรันโค้ด