Skip to content
Snippets Groups Projects
Commit 52cadfbe authored by mariefbourdon's avatar mariefbourdon
Browse files

220615

parent 93e5fbeb
No related branches found
No related tags found
No related merge requests found
No preview for this file type
...@@ -251,8 +251,9 @@ pheno_before_plot ...@@ -251,8 +251,9 @@ pheno_before_plot
pheno_before_zoom <- qtl_plot(pheno_before,lod=data.frame(group = c("alpha=0.05", "alpha=0.1","alpha=0.63"), pheno_before_zoom <- qtl_plot(pheno_before,lod=data.frame(group = c("alpha=0.05", "alpha=0.1","alpha=0.63"),
lod = threshold_before[1:3]), lod = threshold_before[1:3]),
ylab="LOD score", ylab="LOD score",
chrs = "7", chrs = "13",
size=0.6) + size=0.6,
rug = TRUE) +
theme(legend.position = "none", theme(legend.position = "none",
strip.background = element_blank(), strip.background = element_blank(),
strip.text.x = element_blank()) + strip.text.x = element_blank()) +
......
...@@ -163,7 +163,6 @@ pheno_after_plot ...@@ -163,7 +163,6 @@ pheno_after_plot
## Number of markers kept after each function ## Number of markers kept after each function
```{r barplot} ```{r barplot}
# TO DO WHEN ESTMAP OK
none <- tab2 %>% nrow() none <- tab2 %>% nrow()
match <- tab2 %>% filter(exclude_match==0) %>% nrow() match <- tab2 %>% filter(exclude_match==0) %>% nrow()
allele <- tab2 %>% filter(exclude_match==0&exclude_allele==0) %>% nrow() allele <- tab2 %>% filter(exclude_match==0&exclude_allele==0) %>% nrow()
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment