Matriks kontras untuk rata-rata grup
Karena Anda menggunakan parameterisasi rata-rata grup, kini ada satu langkah tambahan dalam pipeline limma. Koefisien sekarang merepresentasikan tingkat ekspresi rata-rata untuk masing-masing dari dua grup sampel, sehingga Anda perlu menentukan kontras kustom untuk menguji perbedaan antara leukemia progresif dan stabil.
Latihan ini merupakan bagian dari kursus
Analisis Ekspresi Diferensial dengan limma di R
Instruksi latihan
Objek eset bertipe ExpressionSet dengan data leukemia dan matriks desain yang baru saja Anda buat (design) telah dimuat ke ruang kerja Anda.
Gunakan
makeContrastsuntuk mendefinisikan kontrasstatus, yaitu selisih tingkat ekspresi rata-rata antara kanker progresif dan stabil.Berikan matriks desain
designke argumenlevels. Ini memungkinkan nama kolomdesign(yang sesuai dengan koefisien) digunakan dalam definisi kontras.
Latihan interaktif langsung praktik
Cobalah latihan ini dengan melengkapi kode contoh ini.
# Load package
library(limma)
# Create a contrasts matrix
cm <- ___(status = ___ - ___,
levels = ___)
# View the contrasts matrix
cm