Preliminary release of moimix

Current development version of moimix - not yet submitted to bioconductor.
v0.0.2 will be released within the next week or so with updated datasets.

moimix: an R package for evaluating multiplicity of infection in malaria parasites

http://bahlolab.github.io/moimix/

Features

• Estimate multiplicity of infection from massively parallel sequencing data
• Estimate heterzygosity and within-isolate diversity directly from read-counts
• Call major alleles within isolates from B-allele frequencies
• Prepare SNP barcode data for use in COIL
• Simulate single nucleotide variant data with known multiplicity of infection

How do I install moimix?

There are plans to put moimix on Bioconductor
in the future, however it is currently only available to install as a development
version from Github:

# install using devtools packages
# first install bioc dependencies
source("https://bioconductor.org/biocLite.R")
biocLite(c("graph", "Rgraphviz", "SeqArray", "SeqVarTools"))
devtools::install_github("bahlolab/moimix")

What data input does moimix require?

moimix makes use of the Genomic Data Storage (GDS) format used
by the Bioconductor package SeqArray
to provide fast access to VCF files in R.

To convert a VCF file to the GDS:

library(SeqArray)
seqVCF2GDS("isolate_snps.vcf.gz", "isolate_snps.gds")

It is also possible to estimate MOI from a matrix of read counts where
the first column is the number of reads supporting the reference allele and the
second column is the number of reads supporting the alternate allele.

How do I use moimix?

See the introduction vignette for usage examples.