Design matrix
การทดลอง doxorubicin เป็น factorial design แบบ 2x2 จึงต้องสร้างตัวแปรรวมเพื่อใช้กับ group-means parametrization
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
การวิเคราะห์ Differential Expression ด้วย limma ใน R
คำแนะนำการฝึกหัด
ExpressionSet object eset ที่บรรจุข้อมูล doxorubicin ถูกโหลดไว้ใน workspace แล้ว
รวมตัวแปร
genotype(WT vs. Top2b null) และtreatment(PBS vs. Dox) ให้เป็น factor variable ตัวเดียวใช้
model.matrixเพื่อสร้าง design matrix แบบไม่มี intercept
แบบฝึกหัดเชิงโต้ตอบแบบลงมือทำ
ลองทำแบบฝึกหัดนี้โดยเติมโค้ดตัวอย่างนี้ให้สมบูรณ์
# Create single variable
group <- with(___(eset), paste(___, ___, sep = "."))
group <- factor(group)
# Create design matrix with no intercept
design <- model.matrix(~___ + ___)
colnames(design) <- levels(group)
# Count the number of samples modeled by each coefficient
colSums(design)