Design matrix
The doxorubicin experiment is a 2x2 factorial design, so you will need to create a combined variable to use in the group-means parametrization.
แบบฝึกหัดนี้เป็นส่วนหนึ่งของหลักสูตร
Differential Expression Analysis with limma in R
คำแนะนำการฝึกหัด
The ExpressionSet object eset with the doxorubicin data has been loaded in your workspace.
Combine the variables
genotype(WT vs. Top2b null) andtreatment(PBS vs. Dox) into a single factor variable.Use
model.matrixto create a design matrix with no 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)