Test per l'espressione differenziale con medie di gruppo
Ora che hai definito la matrice di disegno e la matrice dei contrasti, puoi testare l'espressione differenziale.
Questo esercizio fa parte del corso
Analisi dell'espressione differenziale con limma in R
Istruzioni dell'esercizio
L'oggetto ExpressionSet eset con i dati sulla leucemia, la matrice di disegno (design) e la matrice dei contrasti (cm) sono già caricati nel tuo workspace.
Stima i coefficienti del modello con
lmFit.Applica i contrasti con
contrasts.fit.Calcola le statistiche t con
eBayes.Riassumi i risultati con
decideTests. Non è necessario sotto-selezionarefit2come hai fatto con la parametrizzazione a contrasti di trattamento, perché nel modello a medie di gruppo non c'è il termine d'intercetta.
Esercizio pratico interattivo
Prova a risolvere questo esercizio completando il codice di esempio.
# 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)