Tester l'expression différentielle pour les moyennes de groupe

Maintenant que vous avez défini la matrice de plan expérimental et la matrice de contrastes, vous pouvez tester l'expression différentielle.

Cette activité fait partie du cours

Analyse d'expression différentielle avec limma dans R

Voir le cours

Instructions de l’exercice

L'objet ExpressionSet eset avec les données de leucémie, la matrice de plan expérimental (design) et la matrice de contrastes (cm) ont été chargés dans votre espace de travail.

Ajustez les coefficients du modèle avec lmFit.
Ajustez les contrastes avec contrasts.fit.
Calculez les statistiques t avec eBayes.
Résumez les résultats avec decideTests. Vous n'avez pas besoin de sous-ensemble fit2 comme dans la paramétrisation par contrastes de traitement, car il n'y a pas de terme d'ordonnée à l'origine (interception) dans le modèle à moyennes de groupe.

Exercice interactif pratique

Essayez cet exercice en complétant ce code d’exemple.

# 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)

Modifier et exécuter le code