Skip to content

Commit

Permalink
change path to local files
Browse files Browse the repository at this point in the history
Signed-off-by: Awa Synthia <[email protected]>
  • Loading branch information
awasyn committed Oct 26, 2024
1 parent 57671ce commit 36d0e4f
Showing 1 changed file with 4 additions and 4 deletions.
8 changes: 4 additions & 4 deletions inst/report/scripts/run_molevolvr_pipeline.R
Original file line number Diff line number Diff line change
Expand Up @@ -702,7 +702,7 @@ cleanup_blast <- function(infile_blast, acc2info, prefix, wblast = F) {

# TaxID to lineage mapping
cleanedup_blast$TaxID <- as.integer(cleanedup_blast$TaxID)
lineage_map <- fread("~/awasyn/new_trial/lineage_lookup.txt",
lineage_map <- fread("~/data/lineage_lookup.txt",
header = TRUE, fill = TRUE,
colClasses = lineage_map_cols)

Expand Down Expand Up @@ -841,7 +841,7 @@ run_interproscan <- function(query_file, prefix, outdir) {
# Run InterProScan command
# Construct the command
command <- paste(
"/home/simple/iprdir/interproscan-5.70-102.0/interproscan.sh -i",
"~/interproscan-5.70-102.0/interproscan.sh -i",
shQuote(query_file),
"-b", shQuote(outfile),
"-f TSV --cpu", Sys.getenv("INTERPROSCAN_CPUS", "4"),
Expand Down Expand Up @@ -875,7 +875,7 @@ ipr2lin <- function(ipr, acc2info, prefix) {
ipr_tax <- left_join(ipr_in, acc2info_out, by = "AccNum")

# read in lineage map
lineage_map <- fread("~/awasyn/new_trial/lineage_lookup.txt",
lineage_map <- fread("~/data/lineage_lookup.txt",
header = T, fill = T)

# merge ipr+info w/ lineage
Expand All @@ -889,7 +889,7 @@ ipr2lin <- function(ipr, acc2info, prefix) {
select(-Species.x, -Species.y)

# add lookup table to iprscan file
lookup_tbl <- fread(input = "~/awasyn/new_trial/cln_lookup_tbl.tsv",
lookup_tbl <- fread(input = "~/data/cln_lookup_tbl.tsv",
sep = "\t", header = T, fill = T) %>%
distinct()
if ("AccNum.x" %in% names(ipr_lin)) {
Expand Down

0 comments on commit 36d0e4f

Please sign in to comment.