Testowanie różnicowej ekspresji dla układu czynnikowego 2x2
Nawet jeśli masz więcej kontrastów niż w poprzednich przykładach, testowanie różnicowej ekspresji przy użyciu limma przebiega według tego samego schematu.
To ćwiczenie jest częścią kursu
Analiza różnicowej ekspresji genów z limma w R
Instrukcje do ćwiczenia
Obiekt ExpressionSet eset z danymi Populus, macierz modelu (design) oraz macierz kontrastów (cm) zostały wczytane do twojego środowiska roboczego.
Dopasuj współczynniki modelu za pomocą
lmFit.Dopasuj kontrasty za pomocą
contrasts.fit.Oblicz statystyki t za pomocą
eBayes.Podsumuj wyniki za pomocą
decideTests.
Interaktywne ćwiczenie praktyczne
Spróbuj tego ćwiczenia, uzupełniając ten przykładowy kod.
# 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)