Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

genomic_converter() does not work for with plink ,bed files #193

Open
Npaffen opened this issue Oct 20, 2024 · 0 comments
Open

genomic_converter() does not work for with plink ,bed files #193

Npaffen opened this issue Oct 20, 2024 · 0 comments

Comments

@Npaffen
Copy link

Npaffen commented Oct 20, 2024

When I try to convert my bed-files to the genind-format the converter fails, stating some column is missing. This looks like to be a problem connected to the strata-format.

genomic_converter(data = "data/merged_array_1kg_hgdp_Rfiltered_maf_pruned_nokins_hwe.bed", output = "genind" )
################################################################################
######################### radiator::genomic_converter ##########################
################################################################################
Execution date@time: 20241020@1129
Folder created: -8_radiator_genomic_converter_20241020@1129
Function call and arguments stored in: [email protected]
Filters parameters file generated: [email protected]
Reading PLINK bed file...ed.file

Data summary:
Number of individuals: 8490
Number of markers: 1837423
Error in dplyr::select(): 29.62GB/s, eta: 0s
! Can't select columns that don't exist.
✖ Column INDIVIDUALS doesn't exist.
Run rlang::last_trace() to see where the error occurred.

Computation time, overall: 96 sec

Computation time, overall: 96 sec
######################### completed genomic_converter ##########################
✖ Importing data: plink.bed.file [1m 35.9s]

rlang::last_trace()
<error/vctrs_error_subscript_oob>
Error in dplyr::select():
! Can't select columns that don't exist.
✖ Column INDIVIDUALS doesn't exist.


Backtrace:

  1. ├─radiator::genomic_converter(...)
  2. │ └─radiator::tidy_genomic_data(...)
  3. │ └─radiator::tidy_plink(...)
  4. │ └─radiator::read_plink(...)
  5. │ └─... %$% strata
  6. ├─base::with(., strata)
  7. ├─radiator::read_strata(strata = .)
  8. │ └─strata %<>% dplyr::select(INDIVIDUALS, STRATA)
  9. ├─dplyr::select(., INDIVIDUALS, STRATA)
  10. └─dplyr:::select.data.frame(., INDIVIDUALS, STRATA)
    Run rlang::last_trace(drop = FALSE) to see 17 hidden frames.
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant