error : verbose && p.value.sd < early_A * p_smallest

[wujianming604]

Hi, Rick!

I am using exSTRa for detecting tandem repeat expansion with WES data. I used 15 samples, 13 locus, when I run
tsum <- tsum_test(str_score_1138)
, the error is as follows:

> str_score_1138
exstra_score object with 39221 observations of type named ($data),
  for 15 samples. ($samples)
  Includes associated STR database of 13 loci. ($db)

> tsum <- tsum_test(str_score_1138)
Working on locus DM1
    Reduced replicates to 124.
Working on locus DRPLA
    Reduced replicates to 124.
Working on locus HD
Working on locus HDL2
    Reduced replicates to 249.
Working on locus SBMA
    Reduced replicates to 62.
Working on locus SCA1
    Reduced replicates to 249.
Working on locus SCA2
Error in if (verbose && p.value.sd < early_A * p_smallest) { : 
  missing value where TRUE/FALSE needed
In addition: There were 14 warnings (use warnings() to see them)

Thanks for answering.
Kind regards,
jimmy

[trickytank]

Hi,
Sorry for the late reply. I've only gotten back to this project recently.

This might be hard for me to replicate. If you have the problem again, also run warnings() after.

[trickytank]

You may skip the locus with tsum_test(str_score_1138[locus != "SCA2",]).

I may not be able to replicate the issue without the data.