Error in `is_ref2` in `hetero` case

[dovivi]

Hi,

while running your package, I got an error Error: object 'is_ref2' not found. I then checked impute.R code and saw that there's a problem at line 429:

if(hetero){
        is_ref2
        is_ref <- is_ref[sort(rep(1:length(is_ref),2))]
      }

I guess it should be:

if(hetero){
        is_ref2 <- is_ref[sort(rep(1:length(is_ref),2))]
      }

Would be nice if you could fix it and update the package.

Thanks!

[tpook92]

Thanks for the report. You are 100% correct on this and the new version will fix this issue.

I should note that the algorithm as uploaded here is not really designed for use in non-inbreed material (hetero = TRUE) and will typically result in datasets very similar to the use method for calculation of the auxiliary dataset from Beagle.

An extension on this is currently in early development/plan but there is still a long way to go.