Repo ini menyediakan skrip kami untuk mereproduksi corrplot
korelasi genetik (Gbr. 2) berdasarkan hasil regresi skor LD bivariat (Kanai, M. et al ., Nat. Genet. 2018).
install.packages("dplyr")
devtools::install_github("mkanai/corrplot")
Kami memodifikasi corrplot untuk memvisualisasikan korelasi genetik berpasangan yang diperkirakan melalui regresi skor LD bivariat. Kotak yang lebih besar berhubungan dengan FDR yang lebih signifikan ( corrplot(method = 'psquare', p.mat = p.mat, ...)
). Korelasi yang signifikan (FDR < 0,05) ditunjukkan dengan tanda bintang ( 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
)File ini menyediakan daftar semua korelasi genetik berpasangan yang diperkirakan melalui perangkat lunak ldsc. Skrip mengharapkan semua baris bersifat unik ( yaitu, satu baris untuk setiap pasangan sifat). Bidang yang wajib diisi adalah sebagai berikut:
p1_category
: Kategori sifat dari sifat 1p1
: Sifat 1p2_category
: Kategori sifat dari sifat 2p2
: Sifat 2rg
: Korelasi genetikp
: Nilai-Pq
: nilai q FDRtraitlist.txt
)File ini menyediakan daftar ciri-ciri dan kategorinya. Ini mendefinisikan warna setiap kategori dalam sebuah gambar. Bidang yang wajib diisi adalah sebagai berikut:
CATEGORY
: Kategori sifatTRAIT
: Nama sifatCOLOR
: Warna kategoriContoh keluaran ditunjukkan di bawah ini. Untuk mendapatkan gambar yang dipublikasikan, kami mengedit keluaran pdf menggunakan Adobe Illustrator.
Paket corrplot
asli:
Contoh data dan gambar yang dipublikasikan:
Masahiro Kanai ([email protected])
http://mkanai.github.io/