Code Monkey home page Code Monkey logo

Comments (4)

Imogen-D avatar Imogen-D commented on June 20, 2024

Update:

Running on a smaller subset of data in an interactive session gives a clearer error that appears to be with the stacks input. Help would be much appreciated.

`################################################################################
######################### radiator::genomic_converter ##########################
################################################################################
Execution date@time: 20220623@1130
Folder created: 12_radiator_genomic_converter_20220623@1130
Function call and arguments stored in: [email protected]
Filters parameters file generated: [email protected]

Importing data: vcf.file

Reading VCF...

Data summary:
    number of samples: 17
    number of markers: 4

Filter monomorphic markers
Number of individuals / strata / chrom / locus / SNP:
    Blacklisted: 0 / 0 / 0 / 1 / 1

Strata with low sample size detected: fig <- FALSE


Filter common markers:
Number of individuals / strata / chrom / locus / SNP:
    Blacklisted: 0 / 0 / 0 / 2 / 2

Generating individual stats...
Generating markers stats...
[==================================================] 100%, completed, 0s
[==================================================] 100%, completed, 0s

Number of chromosome/contig/scaffold: 1
Number of locus: 1
Number of markers: 1
Number of strata: 3
Number of individuals: 17

Number of ind/strata:
Y = 15
S = 1
D = 1

Number of duplicate id: 0
radiator Genomic Data Structure (GDS) file: [email protected]

Genotypes formats generated with 1 SNPs:
    GT_BIN (the dosage of ALT allele: 0, 1, 2 NA): TRUE
    GT_VCF (the genotype coding VCFs: 0/0, 0/1, 1/1, ./.): TRUE
    GT_VCF_NUC (the genotype coding in VCFs, but with nucleotides: A/C, ./.): TRUE
    GT (the genotype coding 'a la genepop': 001002, 001001, 000000): TRUE
LOCUS field empty... adding unique id instead
HQ values are all missing: removing column

Tidy data file written: [email protected]

Computation time, overall: 6 sec

Writing tidy data set:
practising.rad

Preparing data for output

Data is bi-allelic
Generating rubias output
Error in data.table::melt.data.table(data = ., id.vars = cols, measure.vars = measure_vars,  :
  One or more values in 'id.vars' is invalid.
Calls: genomic_converter ... %<>% -> <Anonymous> -> %>% -> <Anonymous> -> <Anonymous>
In addition: Warning messages:
1: Unknown or uninitialised column: `LOCUS`.
2: Unknown or uninitialised column: `VARIANT_ID`.

Computation time, overall: 8 sec
######################### completed genomic_converter ##########################
Execution halted`

from radiator.

thierrygosselin avatar thierrygosselin commented on June 20, 2024

Sorry for the long delay
Do you mind sharing the vsf so that I can fix the bug ?
I'll reopen the issue when I have something to work with.

from radiator.

quinn-ca avatar quinn-ca commented on June 20, 2024

I'm getting the same error trying to convert to rubias; my vcf (generated in Stacks) and strata files are attached!
vcf_autosomal_genstr_filtered_maxMeanDP49_rmRelated_snps.vcf.zip
sampleIDs_strata_rmRelated_maxMeanDP49.txt

I've tried the following recommendations in the question: #177 (substituting rubias format). The errors I get depending on the conversion code are:

Error in data.table::melt.data.table(data = ., id.vars = cols, measure.vars = measure_vars, : One or more values in 'id.vars' is invalid.

Error in if (length(markers) != ncol(data)) { : argument is of length zero

from radiator.

quinn-ca avatar quinn-ca commented on June 20, 2024

I realized my strata file wasn't formatted correctly for rubias, so the revised file is attached. This doesn't change the error and I tried the conversion without a strata file (to no avail), but thought the correct file would be helpful if the conversion gets up and running. Thanks!
rubias_strata_file.txt

from radiator.

Related Issues (20)

Recommend Projects

  • React photo React

    A declarative, efficient, and flexible JavaScript library for building user interfaces.

  • Vue.js photo Vue.js

    🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.

  • Typescript photo Typescript

    TypeScript is a superset of JavaScript that compiles to clean JavaScript output.

  • TensorFlow photo TensorFlow

    An Open Source Machine Learning Framework for Everyone

  • Django photo Django

    The Web framework for perfectionists with deadlines.

  • D3 photo D3

    Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

  • javascript

    JavaScript (JS) is a lightweight interpreted programming language with first-class functions.

  • web

    Some thing interesting about web. New door for the world.

  • server

    A server is a program made to process requests and deliver data to clients.

  • Machine learning

    Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.

  • Game

    Some thing interesting about game, make everyone happy.

Recommend Org

  • Facebook photo Facebook

    We are working to build community through open source technology. NB: members must have two-factor auth.

  • Microsoft photo Microsoft

    Open source projects and samples from Microsoft.

  • Google photo Google

    Google ❤️ Open Source for everyone.

  • D3 photo D3

    Data-Driven Documents codes.