Skip to content

Commit

Permalink
minor corrections
Browse files Browse the repository at this point in the history
  • Loading branch information
roddypr committed Dec 6, 2021
1 parent 7a66389 commit 5aa541e
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 8 deletions.
Original file line number Diff line number Diff line change
Expand Up @@ -148,6 +148,8 @@ wc -l results/snp_frequencies.csv
grep -c supergene results/snp_frequencies.csv
# 39210

rm -rf results/by_scaff

```


Expand Down
Original file line number Diff line number Diff line change
Expand Up @@ -191,14 +191,6 @@ private_alleles <- function(freq_matrix, pop) {
}
private_to_invicta_bb <- private_alleles(allele_frequencies_ma, pop="invicta_bb")
private_to_richteri_bb <- private_alleles(allele_frequencies_ma, pop="richteri_bb")
allele_frequencies_ma <- snp_freq[, c("daf_invicta_bb", "daf_invicta_lb", "daf_richteri_bb", "daf_richteri_lb")]
allele_frequencies_ma <- as.matrix(allele_frequencies_ma)
snp_freq$private_to_richteri <- snp_freq$daf_invicta_bb == 0 & snp_freq$daf_invicta_lb == 0 &
(snp_freq$daf_richteri_bb > 0 | snp_freq$daf_richteri_lb > 0 )
Expand Down

0 comments on commit 5aa541e

Please sign in to comment.