It would be nice to have a function that takes multiple SNPs as input. This would allow users with less programming experience to use flexdog() because they wouldn't have to write the for-loop themselves. Updog already Imports foreach and doParallel (because of mupdog()), so I could use those packages to provide parallelization.
data.frame Format
A possible format would be to require users to have SNPs in a data frame
The output should be a list-like format with two data.frames.
The first data.frame contains summary information for each SNP x Individual (posterior means, posterior modes, maximum posterior probabilities, posterior probabilities, id's, snp names, counts, size).
The second data.frame contains summary information for each SNP (bias, sequencing error rate, overdispersion parameter, log-likelihood, estimated proportion missing).
Description
It would be nice to have a function that takes multiple SNPs as input. This would allow users with less programming experience to use
flexdog()
because they wouldn't have to write the for-loop themselves. Updog already Imports foreach and doParallel (because ofmupdog()
), so I could use those packages to provide parallelization.data.frame
FormatA possible format would be to require users to have SNPs in a data frame
We could call this
flexdog_df()
.VariantAnnotation
FormatWe could also allow the user to insert a
VariantAnnotation
(https://doi.org/doi:10.18129/B9.bioc.VariantAnnotation) object. We could call thisflexdog_va()
.Output
The output should be a list-like format with two
data.frame
s.The first
data.frame
contains summary information for each SNP x Individual (posterior means, posterior modes, maximum posterior probabilities, posterior probabilities, id's, snp names, counts, size).The second
data.frame
contains summary information for each SNP (bias, sequencing error rate, overdispersion parameter, log-likelihood, estimated proportion missing).