Testowanie różnicowej ekspresji dla 3 grup
Teraz, gdy masz już zdefiniowaną macierz projektu i macierz kontrastów, możesz przystąpić do testowania różnicowej ekspresji.
To ćwiczenie jest częścią kursu
Analiza różnicowej ekspresji genów z limma w R
Instrukcje do ćwiczenia
Obiekt ExpressionSet eset z danymi dotyczącymi hipoksji, macierz projektu (design) oraz macierz kontrastów (cm) zostały wczytane do twojego środowiska pracy.
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)