Ce dépôt fournit notre script pour reproduire corrplot
de corrélation génétique (Fig. 2) sur la base des résultats de régression bivariée du score LD (Kanai, M. et al ., Nat. Genet. 2018).
install.packages("dplyr")
devtools::install_github("mkanai/corrplot")
Nous avons modifié corrplot pour visualiser les corrélations génétiques par paires estimées via une régression bivariée du score LD. Les carrés plus grands correspondent à des FDR plus significatifs ( corrplot(method = 'psquare', p.mat = p.mat, ...)
). Les corrélations significatives (FDR < 0,05) sont indiquées par des astérisques ( 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
)Ce fichier fournit une liste de toutes les corrélations génétiques par paires estimées via le logiciel ldsc. Le script s'attend à ce que toutes les lignes soient uniques ( c'est-à-dire une ligne pour chaque paire de traits). Les champs obligatoires sont les suivants :
p1_category
: Catégorie de trait du trait 1p1
: Caractère 1p2_category
: Catégorie de trait du trait 2p2
: Caractère 2rg
: Corrélation génétiquep
: valeur Pq
: valeur q du FDRtraitlist.txt
)Ce fichier fournit une liste de traits et leurs catégories. Il définit une couleur de chaque catégorie dans une figure. Les champs obligatoires sont les suivants :
CATEGORY
: Catégorie de traitTRAIT
: Nom du caractèreCOLOR
: Couleur de la catégorieUn exemple de sortie est présenté ci-dessous. Pour obtenir la figure publiée, nous avons édité une sortie PDF à l'aide d'Adobe Illustrator.
Le package corrplot
d'origine :
Les exemples de données et le chiffre publié :
Masahiro Kanai ([email protected])
http://mkanai.github.io/