snpnet icon indicating copy to clipboard operation
snpnet copied to clipboard

Error: Problem with `mutate()` column `stats_means`.

Open julibeg opened this issue 3 years ago • 0 comments

When trying to run snpnet, computeStats fails with

Error: Problem with `mutate()` column `stats_means`.
ℹ `stats_means = (HAP_ALT_CTS + HET_REF_ALT_CTS + 2 * TWO_ALT_GENO_CTS)/OBS_CT`.
✖ non-numeric argument to binary operator

My guess is that this is because of the presence of joined multiallelic sites causing PLINK to write cells with comma-separated numbers to some of the columns of the output of --geno-counts. In order to avoid the issue, I split the merged variants. However, then snpnet fails with

Error: --read-freq variant ID 'rs112422003' appears multiple times in main

So it looks like the user needs to make sure that the variants are split and have unique IDs. Is that correct? If so, then it would be nice if you could add this to the documentation somewhere.

julibeg avatar Jul 29 '21 14:07 julibeg