diff --git a/figure/4_ggplotgene_gwasPS_new-1.pdf b/figure/4_ggplotgene_gwasPS_new-1.pdf index f6aa18e4f761cd53ea7a203d24cd77bac4697cf3..3975dd94d78a24948b508c44e29deb1ed005956b 100644 Binary files a/figure/4_ggplotgene_gwasPS_new-1.pdf and b/figure/4_ggplotgene_gwasPS_new-1.pdf differ diff --git a/rnw_scripts/covid_comp_gwas-copy.Rnw b/rnw_scripts/covid_comp_gwas-copy.Rnw index 15c66d9abe781c0472eecefcd9f18be407a0b349..9204a90692e01ebebe130ce0810b9cb8938b6e18 100755 --- a/rnw_scripts/covid_comp_gwas-copy.Rnw +++ b/rnw_scripts/covid_comp_gwas-copy.Rnw @@ -317,7 +317,7 @@ PS<-PS*3 Updated positions under PS <<>>= -PS<-c(448, 472, 553, 930) +PS<-c(447, 471, 552, 928) PS<-PS*3 @ diff --git a/rnw_scripts/covid_comp_gwas-copy.pdf b/rnw_scripts/covid_comp_gwas-copy.pdf index 92a7d68880db5c854ae40855380e84880594cfc8..19fdf4e41add903c02ee20a298c97db6c819051a 100644 Binary files a/rnw_scripts/covid_comp_gwas-copy.pdf and b/rnw_scripts/covid_comp_gwas-copy.pdf differ diff --git a/rnw_scripts/covid_comp_gwas-copy.tex b/rnw_scripts/covid_comp_gwas-copy.tex index 388a3cf925ecee92ea6a3deb25633cfbcfb51bff..8f9ac032bc9b3395afb4ead791d5cc61bcfad565 100644 --- a/rnw_scripts/covid_comp_gwas-copy.tex +++ b/rnw_scripts/covid_comp_gwas-copy.tex @@ -479,7 +479,7 @@ Updated positions under PS \begin{knitrout} \definecolor{shadecolor}{rgb}{0.969, 0.969, 0.969}\color{fgcolor}\begin{kframe} \begin{alltt} -\hlstd{PS}\hlkwb{<-}\hlkwd{c}\hlstd{(}\hlnum{448}\hlstd{,} \hlnum{472}\hlstd{,} \hlnum{553}\hlstd{,} \hlnum{930}\hlstd{)} +\hlstd{PS}\hlkwb{<-}\hlkwd{c}\hlstd{(}\hlnum{447}\hlstd{,} \hlnum{471}\hlstd{,} \hlnum{552}\hlstd{,} \hlnum{928}\hlstd{)} \hlstd{PS}\hlkwb{<-}\hlstd{PS}\hlopt{*}\hlnum{3} \end{alltt} @@ -692,12 +692,12 @@ Conversion des coordonnées dans example\_fyco1 ## posref gwas1 gwas2 gwas3 pos_aln PS ## 48433 45966331 FALSE FALSE TRUE 3009 FALSE ## 48435 45966333 FALSE FALSE TRUE 3007 FALSE -## 48652 45966550 FALSE FALSE FALSE 2790 TRUE +## 48658 45966556 FALSE FALSE FALSE 2784 TRUE ## 48697 45966595 FALSE FALSE TRUE 2745 FALSE -## 49780 45967678 FALSE FALSE FALSE 1659 TRUE -## 50023 45967921 FALSE FALSE FALSE 1416 TRUE -## 50095 45967993 FALSE FALSE FALSE 1344 TRUE +## 49783 45967681 FALSE FALSE FALSE 1656 TRUE +## 50026 45967924 FALSE FALSE FALSE 1413 TRUE ## 50097 45967995 FALSE FALSE TRUE 1342 FALSE +## 50098 45967996 FALSE FALSE FALSE 1341 TRUE ## 50145 45968043 FALSE TRUE FALSE 1294 FALSE ## 50617 45968515 FALSE FALSE TRUE 819 FALSE \end{verbatim} @@ -705,24 +705,30 @@ Conversion des coordonnées dans example\_fyco1 \hlstd{gwastab[gwastab}\hlopt{$}\hlstd{all_inv_var_meta_p}\hlopt{<}\hlstd{thres} \hlopt{&} \hlstd{gwastab}\hlopt{$}\hlstd{POS} \hlopt{%in%} \hlkwd{c}\hlstd{(}\hlnum{45966331}\hlstd{,} \hlnum{45966333}\hlstd{,} \hlnum{45966595}\hlstd{,} \hlnum{45967995}\hlstd{,} \hlnum{45968515}\hlstd{),]} \end{alltt} \begin{verbatim} -## CHR POS REF ALT SNP all_meta_N all_inv_var_meta_beta all_inv_var_meta_sebeta -## 1428785 3 45966331 G T 3:45966331:G:T 23 0.38699 0.027678 -## 1428786 3 45966333 T C 3:45966333:T:C 23 0.37298 0.027352 -## 1428788 3 45966595 G A 3:45966595:G:A 23 0.38420 0.027682 -## 1428791 3 45967995 G A 3:45967995:G:A 23 0.37943 0.027617 -## 1428794 3 45968515 T C 3:45968515:T:C 23 0.37288 0.027339 -## all_inv_var_meta_p all_inv_var_meta_cases\ttotal all_inv_var_meta_controls -## 1428785 2.0106e-44 8475 999425 -## 1428786 2.4293e-42 8475 999425 -## 1428788 8.4910e-44 8475 999425 -## 1428791 5.9274e-43 8475 999425 -## 1428794 2.3510e-42 8475 999425 -## all_inv_var_meta_effective all_inv_var_het_p all_meta_AF rsid -## 1428785 7110 0.00081500 0.1221 rs13079478 -## 1428786 7110 0.00107150 0.1230 rs13059238 -## 1428788 7110 0.00155590 0.1222 rs13079869 -## 1428791 7110 0.00039328 0.1224 rs33910087 -## 1428794 7110 0.00111540 0.1230 rs13071283 +## CHR POS REF ALT SNP all_meta_N all_inv_var_meta_beta +## 1428785 3 45966331 G T 3:45966331:G:T 23 0.38699 +## 1428786 3 45966333 T C 3:45966333:T:C 23 0.37298 +## 1428788 3 45966595 G A 3:45966595:G:A 23 0.38420 +## 1428791 3 45967995 G A 3:45967995:G:A 23 0.37943 +## 1428794 3 45968515 T C 3:45968515:T:C 23 0.37288 +## all_inv_var_meta_sebeta all_inv_var_meta_p all_inv_var_meta_cases\ttotal +## 1428785 0.027678 2.0106e-44 8475 +## 1428786 0.027352 2.4293e-42 8475 +## 1428788 0.027682 8.4910e-44 8475 +## 1428791 0.027617 5.9274e-43 8475 +## 1428794 0.027339 2.3510e-42 8475 +## all_inv_var_meta_controls all_inv_var_meta_effective all_inv_var_het_p all_meta_AF +## 1428785 999425 7110 0.00081500 0.1221 +## 1428786 999425 7110 0.00107150 0.1230 +## 1428788 999425 7110 0.00155590 0.1222 +## 1428791 999425 7110 0.00039328 0.1224 +## 1428794 999425 7110 0.00111540 0.1230 +## rsid +## 1428785 rs13079478 +## 1428786 rs13059238 +## 1428788 rs13079869 +## 1428791 rs33910087 +## 1428794 rs13071283 \end{verbatim} \end{kframe} \end{knitrout} @@ -776,8 +782,8 @@ dbSNP (rs13079478) \hlstd{fyco1tab[fyco1tab}\hlopt{$}\hlstd{posref}\hlopt{==}\hlnum{45966556}\hlstd{,]} \end{alltt} \begin{verbatim} -## posref gwas1 gwas2 gwas3 pos_aln PS -## 48658 45966556 FALSE FALSE FALSE 2784 FALSE +## posref gwas1 gwas2 gwas3 pos_aln PS +## 48658 45966556 FALSE FALSE FALSE 2784 TRUE \end{verbatim} \begin{alltt} \hlkwd{dim}\hlstd{(gwastab)} @@ -803,18 +809,20 @@ dbSNP (rs13079478) ## [7] all_inv_var_meta_beta all_inv_var_meta_sebeta all_inv_var_meta_p ## [10] all_inv_var_meta_cases\ttotal all_inv_var_meta_controls all_inv_var_meta_effective ## [13] all_inv_var_het_p all_meta_AF rsid -## <0 rows> (or 0-length row.names) +## <0 lignes> (ou 'row.names' de longueur nulle) \end{verbatim} \begin{alltt} \hlstd{gwastab[(gwastab}\hlopt{$}\hlstd{POS}\hlopt{>}\hlnum{45966500} \hlopt{&} \hlstd{gwastab}\hlopt{$}\hlstd{POS}\hlopt{<}\hlnum{45966600}\hlstd{),]} \end{alltt} \begin{verbatim} -## CHR POS REF ALT SNP all_meta_N all_inv_var_meta_beta all_inv_var_meta_sebeta -## 1428788 3 45966595 G A 3:45966595:G:A 23 0.3842 0.027682 -## all_inv_var_meta_p all_inv_var_meta_cases\ttotal all_inv_var_meta_controls -## 1428788 8.491e-44 8475 999425 -## all_inv_var_meta_effective all_inv_var_het_p all_meta_AF rsid -## 1428788 7110 0.0015559 0.1222 rs13079869 +## CHR POS REF ALT SNP all_meta_N all_inv_var_meta_beta +## 1428788 3 45966595 G A 3:45966595:G:A 23 0.3842 +## all_inv_var_meta_sebeta all_inv_var_meta_p all_inv_var_meta_cases\ttotal +## 1428788 0.027682 8.491e-44 8475 +## all_inv_var_meta_controls all_inv_var_meta_effective all_inv_var_het_p all_meta_AF +## 1428788 999425 7110 0.0015559 0.1222 +## rsid +## 1428788 rs13079869 \end{verbatim} \end{kframe} \end{knitrout}