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