Domanda:
Come selezionare una potenza per una rete con topologia senza scala
llrs
2017-05-22 20:02:28 UTC
view on stackexchange narkive permalink

In un'analisi di rete di coespressione genica ponderata (utilizzando WGCNA), la potenza della soglia morbida è consigliata come filtro del rumore. Consiste nell'elevare la correlazione a un certo numero. Per decidere questa potenza si stima la topologia senza scala per alcune potenze. La funzione per stimarlo stampa:

  Power SFT.R.sq pendenza troncata.R.sq mean.k. median.k. max k.1 1 0.9300 3.110 0.996 2960.0 3060.0 39702 2 0.7510 1.010 0.964 1750.0 1780.0 29003 3 0.1730 0.258 0.806 1170.0 1150.0 22804 4 0.0942 -0.183 0.713 833.0 782.0 18705 5 0.3800 -0.463 0.777 623.0 559.0 15806.656 0.5350- 7 0.6270 -0.797 0.872 381.0 312.0 11908 8 0.6870 -0.910 0.900 307.0 241.0 10509 9 0.7270 -1.000 0.918 252.0 189.0 93610 10 0.7490 -1.080 0.928 210.0 150.0 84111 12 0.7850 -1.190 0.948 150.0 98.0 69312 14 0.8090 -1.280 0.958 -1.360 0.968 84.0 45.6 49714 18 0.8410 -1.410 0.973 65.2 32.1 42915 20 0.8490 -1.450 0.977 51.6 23.0 375  

La raccomandazione s delle FAQ indicano che un valore SFT.R.sq deve essere superiore

0,8 per potenze ragionevoli (inferiore a 15 per reti ibride non firmate o firmate e inferiore di 30 per le reti firmate)

e la connettività media al di sotto delle centinaia.

Altri hanno utilizzato la potenza proprio come filtro del rumore senza preoccuparsi molto della topologia senza scala in forma. Sceglierei la prima potenza anche se la connettività media è dell'ordine di migliaia perché l'adattamento della topologia senza scala è piuttosto alto, tuttavia la pendenza mi lascia perplesso.

Come dovrebbe essere la potenza della soglia morbida selezionato?

Basato su una domanda

di esempio
Una risposta:
gringer
2017-06-08 20:27:39 UTC
view on stackexchange narkive permalink

Sarei preoccupato per i dati di input che hanno prodotto questi numeri. Questo set di dati ha una forma davvero strana con un calo all'inizio:

  data.df <- read.table ("soft_threshold.txt", header = TRUE); png ("scale-free -194.png "); plot (data.df $ SFT.R.sq); dummy <- dev.off ();  

scale-free power with an odd inflection

Ecco come mi aspetto che sia questo (dalle informazioni supplementari in questo documento):

normal scale-free topology

È più facile selezionare la potenza / soglia quando la trama sembra la seconda situazione. Con i dati così come presentati, non sembra che nessun valore sarebbe adatto. ma se fossi costretto a scegliere, sceglierei 12 perché è il primo numero sopra la soglia magica (ignorando la stranezza all'inizio). Tuttavia, ammetto che non capisco i concetti alla base del potere e sto solo seguendo le istruzioni consigliate, quindi non posso spiegare perché è quello che è.

Sì, sono a conoscenza di questa particolarità, penso sia dovuto a un bug nel pacchetto WGCNA, quando pickSoftThreshold non utilizza i dati di espressione. Probabilmente avrei dovuto scegliere un dato più tipico. Ma il semplice taglio a una certa soglia lascia fuori la considerazione sulla parte media della connettività. (Inoltre sembra trattare lo stesso che la pendenza della potenza = 1 è positiva e le altre diventano negative)


Questa domanda e risposta è stata tradotta automaticamente dalla lingua inglese. Il contenuto originale è disponibile su stackexchange, che ringraziamo per la licenza cc by-sa 3.0 con cui è distribuito.
Loading...