Skip to content

Commit

Permalink
Wrong delimiter specified for featureCounts CORE files
Browse files Browse the repository at this point in the history
  • Loading branch information
isaacvock committed Aug 26, 2024
1 parent 69ce3c8 commit 59c779e
Showing 1 changed file with 4 additions and 35 deletions.
39 changes: 4 additions & 35 deletions workflow/scripts/merge_features_and_muts.R
Original file line number Diff line number Diff line change
Expand Up @@ -82,16 +82,8 @@ if(opt$genes){
message("genes_file is:")
print(genes_file)

if(use_readr){
genes <- fread(genes_file)

genes <- read_csv(genes_file)


}else{

genes <- fread(genes_file)

}

colnames(genes) <- c("qname", "status", "nhits", "GF")

Expand All @@ -117,16 +109,8 @@ if(opt$exons){
message("exons_file is:")
print(exons_file)

if(use_readr){

exons <- read_csv(exons_file)


}else{

exons <- fread(exons_file)
exons <- fread(exons_file)

}

colnames(exons) <- c("qname", "status", "nhits", "XF")

Expand All @@ -148,15 +132,8 @@ if(opt$exonbins){
exonbins_file <- list.files("./results/featurecounts_exonbins/",
pattern = sample, full.names = TRUE)[1]

if(use_readr){
exonbins <- fread(exonbins_file)

exonbins <- read_csv(exonbins_file)

}else{

exonbins <- fread(exonbins_file)

}

colnames(exonbins) <- c("qname", "status", "nhits", "exon_bin")

Expand All @@ -181,16 +158,8 @@ if(opt$transcripts){
transcripts_file <- list.files("./results/featurecounts_transcripts/",
pattern = sample, full.names = TRUE)[1]

if(use_readr){

transcripts <- read_csv(transcripts_file)


}else{

transcripts <- fread(transcripts_file)
transcripts <- fread(transcripts_file)

}

colnames(transcripts) <- c("qname", "status", "nhits", "transcripts")

Expand Down

0 comments on commit 59c779e

Please sign in to comment.