On Mon, 13 Oct 2014 08:54:51 PM Kate Ignatius wrote:
> Just an update to this:
>
> gtal <- function(d) {
> alleles <- sapply(d, function(.) strsplit(as.character(.), "/"))
> gt <- unlist(lapply(alleles, function(x)
>ifelse(identical(x[[1]], vcf[,3]) & identical(x[[2]], vcf[,3])
Just an update to this:
gtal <- function(d) {
alleles <- sapply(d, function(.) strsplit(as.character(.), "/"))
gt <- unlist(lapply(alleles, function(x)
ifelse(identical(x[[1]], vcf[,3]) & identical(x[[2]], vcf[,3]), 'RR',
ifelse(identical(x[[1]], vcf[,4]) & identical(
2 matches
Mail list logo