2 grup arasında fark ifade testi
Artık bir ExpressionSet nesnesine ve bir tasarım matrisine sahip olduğuna göre, modelini limma ile kurup test edebilirsin.
Bu egzersiz
R ile limma kullanarak Diferansiyel Ekspresyon Analizi
kursunun bir parçasıdırEgzersiz talimatları
Lösemi verilerini içeren eset adlı ExpressionSet nesnesi ve az önce oluşturduğun tasarım matrisi (design) çalışma alanına yüklendi.
Model katsayılarını
lmFitile uydur.t-istatistiklerini
eBayesile hesapla.Sonuçları
decideTestsile özetle.Sabit terimi (intercept) test etmek mantıklı olmadığından,
decideTestsçalıştırırken katsayı olarak"Diseasestable"belirt.
Uygulamalı interaktif egzersiz
Bu örnek kodu tamamlayarak bu egzersizi bitirin.
# Load package
library(limma)
# Fit the model
fit <- ___(eset, ___)
# Calculate the t-statistics
fit <- ___(fit)
# Summarize results
results <- ___(fit[, "___"])
summary(results)