Quale test statistico dovrebbe essere usato per testare l'arricchimento degli elenchi genici?


10

Ho eseguito un esperimento per testare la sensibilità cellulare a un determinato agente di danno del DNA. Abbiamo trovato 270 geni specificamente sensibili al farmaco e il numero totale di geni analizzati era 3668. 38 dei 270 geni sensibili sono classificati come "geni di riparazione del DNA". Se il numero di "geni di riparazione del DNA" contenuti nel genoma è 112 e il numero totale di geni nel genoma è 3668, l'arricchimento dei geni sensibili è nei geni di riparazione del DNA? Quale test statistico dovrebbe essere usato? Gradirei se potessi anche dirmi qualche strumento per calcolare il valore p online.

Risposte:


17

La pratica standard per verificare l'arricchimento degli elenchi genetici consiste nell'eseguire un test ipergeometrico o, equivalentemente, un test esatto di Fisher unilaterale . Hai la seguente tabella di contingenza :2×2

DNA RepairOtherSensitive38232270Not Sensitive74332433981123556

È possibile eseguire il test Rcome segue:

fisher.test(matrix(c(38,74,232,3324),nrow=2,ncol=2),alternative="greater")

Il che dà un risultato molto significativo:

Fisher's Exact Test for Count Data

data:  matrix(c(38, 74, 232, 3324), nrow = 2, ncol = 2) 
p-value < 2.2e-16
alternative hypothesis: true odds ratio is greater than 1 
95 percent confidence interval:
5.062107      Inf 
sample estimates:
odds ratio 
7.34918

Si noti che mentre stiamo testando la sovrarappresentazione (anziché la sottorappresentazione) il alternativeparametro è impostato su "greater".


1
Grazie mille per la tua risposta. Ho anche pensato che il test esatto di Fisher potesse essere un buon metodo per l'analisi. Non ho alcun software statistico per eseguire i risultati per altre classi funzionali che vorrei testare anch'io. Conosci qualche strumento "online" per ottenere i valori con tutti i decimali?
Laura,

1
Puoi scaricare R gratuitamente. Vedi r-project.org Quindi non avere software è solubile (e pensare che tu abbia bisogno di un modo di calcolare online non è corretto). Ma ti preghiamo di fare un po 'di ricerca per scoprire queste cose da solo. Consulta i consigli su stats.stackexchange.com/help/how-to-ask su come porre una buona domanda.
Nick Cox,

@Nick Il tuo consiglio è buono, ma per favore non considerarlo come una caratterizzazione del poster: tale frase è fin troppo facilmente fraintesa come un attacco, che dubito tu abbia inteso. Pertanto ho rimosso la frase preliminare nel tuo commento (che non ha aggiunto alcuna informazione).
whuber

Un ottimo strumento online per questo è: mathcelebrity.com/fishers_exact_test.php

Potresti spiegare ulteriormente, la sovrarappresentazione di cosa viene calcolato esattamente?
sdgaw erzswer,
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.