Looking the code for the group.b.wcnt test, the weights was weights <- 1/sqrt(NS[vids]*MAF[vids]*(1-MAF[vids])), does this actually should be weights <- 1/sqrt(MAF[vids]*(1-MAF[vids])). As NS[vids] is the number of called samples for each marker, weight by this way will give more weight to high missing rate variant, seems doesn't make sense.
Hi Developers,
Looking the code for the group.b.wcnt test, the weights was
weights <- 1/sqrt(NS[vids]*MAF[vids]*(1-MAF[vids]))
, does this actually should beweights <- 1/sqrt(MAF[vids]*(1-MAF[vids]))
. As NS[vids] is the number of called samples for each marker, weight by this way will give more weight to high missing rate variant, seems doesn't make sense.Thank you very much! Best regards Wallace