Are the following lines meant to tabulate the counts on both alleles for each locus?
compute_counts = function(data,locus){
return(table(data[,1+2_locus, 2+2_locus]))
}
The code currently only tabulates the counts on the first allele only.
I believe the code should actually be
compute_counts = function(data,locus){
return(table(data[,1+2_locus]) + table(data[,2+2_locus]))
}
Are the following lines meant to tabulate the counts on both alleles for each locus? compute_counts = function(data,locus){ return(table(data[,1+2_locus, 2+2_locus])) }
The code currently only tabulates the counts on the first allele only. I believe the code should actually be compute_counts = function(data,locus){ return(table(data[,1+2_locus]) + table(data[,2+2_locus])) }