Este repositorio proporciona nuestro script para reproducir corrplot
de correlación genética (Fig. 2) basado en resultados de regresión de puntuación LD bivariada (Kanai, M. et al ., Nat. Genet. 2018).
install.packages("dplyr")
devtools::install_github("mkanai/corrplot")
Modificamos corrplot para visualizar correlaciones genéticas por pares estimadas mediante regresión de puntuación LD bivariada. Los cuadrados más grandes corresponden a FDR más significativos corrplot(method = 'psquare', p.mat = p.mat, ...)
). Las correlaciones significativas (FDR < 0,05) se indican mediante asteriscos ( corrplot(sig = 'pch', sig.level = 0.05, pch = '*')
).
Rscript plot_corrplot_rg.R input_example/input_rg.txt input_example/traitlist.txt
input_rg.txt
)Este archivo proporciona una lista de todas las correlaciones genéticas por pares estimadas mediante el software ldsc. El script espera que todas las filas sean únicas ( es decir, una fila por cada par de rasgos). Los campos obligatorios son los siguientes:
p1_category
: categoría de rasgo del rasgo 1p1
: rasgo 1p2_category
: categoría de rasgo del rasgo 2p2
: rasgo 2rg
: correlación genéticap
: valor pq
: valor q de FDRtraitlist.txt
)Este archivo proporciona una lista de rasgos y sus categorías. Define un color de cada categoría en una figura. Los campos obligatorios son los siguientes:
CATEGORY
: Categoría de rasgoTRAIT
: Nombre del rasgoCOLOR
: Color de categoríaA continuación se muestra un resultado de ejemplo. Para obtener la figura publicada, editamos una salida en PDF con Adobe Illustrator.
El paquete original corrplot
:
Los datos de ejemplo y la figura publicada:
Masahiro Kanai ([email protected])
http://mkanai.github.io/