bodkan/admixr

Another duplicated populations error

Closed this issue · 0 comments

Reported over email:

I just wanted to comment on an error that I've had twice already when using qpAdm with Admixr: As far as I can tell, it has to do with using populations with the same characters as source and outgroup.

For example, it happened when using the curated dataset from the Reichlab with these settings:

  1. target: Germany_Corded_Ware ; sources: Germany_Corded_Ware.SG

  2. target: Israel_MLBA; sources: Israel_MLBA_published.

It was solved by changing the name of one population.

This is the error message:

Backtrace:
     █
  1. └─admixr::qpAdm(...)
  2.   └─base::lapply(...)
  3.     └─admixr:::FUN(X[[i]], ...)
  4.       └─admixr::read_output(files[["log_file"]])
  5.         └─parsers[[cmd]](log_lines, ...)
  6.           └─`%>%`(...)
  7.             ├─base::withVisible(eval(quote(`_fseq`(`_lhs`)), env, env))
  8.             └─base::eval(quote(`_fseq`(`_lhs`)), env, env)
  9.               └─base::eval(quote(`_fseq`(`_lhs`)), env, env)
 10.                 └─admixr:::`_fseq`(`_lhs`)
 11.                   └─magrittr::freduce(value, `_function_list`)
 12.                     ├─base::withVisible(function_list[[k]](value))
 13.                     └─function_list[[k]](value)
 14.                       └─dplyr::mutate_at(., dplyr::var
In addition: Warning message:
In function_list[[k]](value) : NAs introduced by coercion
Execution halted