#SeqSQC error: Problem while converting vcf format to gds format
0
0
Entering edit mode
Bhargavi • 0
@dd5d7bd4
Last seen 3.2 years ago
India

Enter the body of text here

An error occurred when I tried to convert our sample vcf to GDS format. Even after providing the read, write, and edit permissions for all the input files it is showing the warning. Any thoughts as to where I am going wrong?

Code should be placed in three backticks as shown below


# include your problematic code here with any corresponding output 
# please also include the results of running the following in an R session 

sessionInfo( )

Set working directory

setwd("D:/R_workspace/WIPRO3/Wipro-3/extdata/") library(vcfR) library(BiocGenerics) library(gdsfmt) library(SNPRelate) library(SeqSQC)

Load the input files

infile <- system.file("extdata", "examples_sub.vcf", package="SeqSQC") sample.annot <- system.file("extdata", "UIMA-F-SE8-2nd-S21.txt", package="SeqSQC") cr <- system.file("extdata", "SureSelectXT_V8_Covered.bed", package="SeqSQC") outdir <- "D:/R_workspace/WIPRO3/Wipro-3" outfile <- file.path(outdir, "testWrapUp") outfile [1] "D:/R_workspace/WIPRO3/Wipro-3/testWrapUp" seqfile <- sampleQC(vfile = infile, output = outfile, capture.region = cr, sample.annot = sample.annot, format.data = "NGS", format.file = "vcf", QCreport = FALSE) Load vcf file ... Start file conversion from VCF to SNP GDS ... Method: exacting biallelic SNPs Error in VCF_SampID(vcf.fn[i]) : Error VCF format: invalid sample id! In addition: Warning message: In file(filename, "rt") : file("") only supports open = "w+" and open = "w+b": using the former

SNPRelate gdsfmt SeqSQC • 882 views
ADD COMMENT

Login before adding your answer.

Traffic: 477 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6