diff --git a/DESCRIPTION b/DESCRIPTION index c2e3a8d..fd664db 100644 --- a/DESCRIPTION +++ b/DESCRIPTION @@ -1,6 +1,6 @@ Package: antigen.garnish Title: Tumor neoantigen prediction -Version: 2.3.0 +Version: 2.3.1 Authors@R: c(person(given = c("Lee", "P."), family = "Richman", diff --git a/NEWS.md b/NEWS.md index f8ebd86..4ec2010 100644 --- a/NEWS.md +++ b/NEWS.md @@ -1,3 +1,9 @@ +# antigen.garnish 2.3.1 + +Bugfix: + +- errors parsing rare MHCII alleles for netMHC tools + # antigen.garnish 2.3.0 Bugfix: diff --git a/README.md b/README.md index 082a43a..0550966 100644 --- a/README.md +++ b/README.md @@ -26,9 +26,9 @@ Two methods exist to run `antigen.garnish`: ### Docker ```sh -docker pull andrewrech/antigen.garnish:2.3.0 +docker pull andrewrech/antigen.garnish:2.3.1 -cID=$(docker run -it -d andrewrech/antigen.garnish:2.3.0 /bin/bash) +cID=$(docker run -it -d andrewrech/antigen.garnish:2.3.1 /bin/bash) ``` [Download](https://services.healthtech.dtu.dk/software.php) netMHC binaries (academic license): NetMHC 4.0, NetMHCpan 4.1b, NetMHCII 2.3, NetMHCIIpan 4.0. diff --git a/antigen.garnish.pdf b/antigen.garnish.pdf new file mode 100644 index 0000000..ea6afdf Binary files /dev/null and b/antigen.garnish.pdf differ diff --git a/inst/extdata/src/antigen.garnish.flowchart.R b/inst/extdata/src/antigen.garnish.flowchart.R deleted file mode 100755 index 0142016..0000000 --- a/inst/extdata/src/antigen.garnish.flowchart.R +++ /dev/null @@ -1,29 +0,0 @@ - antigen.garnish_flowchart <- glue::glue("digraph boxes_and_circles {{ - - graph [overlap = true, bgcolor = transparent, rankdir = LR] - - node [margin = 0.2, shape = box, fontsize = 18, fontname = '{Font}', fontcolor = '#000000', color = '#000000', fillcolor = '#82b1ff:#8c9eff', style = radial] - - input [label = 'Somatic variants\n(VCF file)\nor peptides'] - epi [label = 'Top predicted\nneoepitopes'] - - node [margin = 0.2, face = bold, shape = box, fontname = '{Font}', fontcolor = '#000000', fontsize = 22, color = '#000000', fillcolor = '#ffd180:#ff9100', style = radial] - - ag [label = 'antigen.garnish'] - - input -> ag -> epi - - edge [penwidth = 1.5, color = '#000000', arrowsize = 1] - - }}") - - antigen.garnish_flowchart %<>% DiagrammeR::grViz(.) - write(DiagrammeRsvg::export_svg(antigen.garnish_flowchart), - file = paste0( - DROPBOX, - "S3/antigen.garnish_flowchart.svg" - ), - ncolumns = 1, - append = FALSE, - sep = " " - ) diff --git a/inst/extdata/src/antigen.garnish_transcript_db.R b/inst/extdata/src/antigen.garnish_transcript_db.R index 0f20d31..97f32e4 100755 --- a/inst/extdata/src/antigen.garnish_transcript_db.R +++ b/inst/extdata/src/antigen.garnish_transcript_db.R @@ -114,7 +114,7 @@ make_refseq_db <- function(txt) { 1:length(l) %>% lapply(function(i) { print(i) ids <- l[[i]] %>% paste(collapse = ",") - glue::glue("http 'https://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=nuccore&id={ids}&retmode=xml' > rs{i}.xml") %>% system() + paste0("http 'https://eutils.ncbi.nlm.nih.gov/entrez/eutils/efetch.fcgi?db=nuccore&id=", ids, "&retmode=xml' > rs", i, ".xml") %>% system() }) # extract RefSeq IDm cDNA, and peptide sequence from .xml @@ -140,7 +140,7 @@ make_refseq_db <- function(txt) { sub <- x %>% # move to sequence record ancestor corresponding to transcript ID - xml2::xml_find_all(glue::glue("//GBSeq_accession-version[text()='{id}']/ancestor::GBSeq")) + xml2::xml_find_all(paste0("//GBSeq_accession-version[text()='", id, "']/ancestor::GBSeq")) peptide <- sub %>% # find translation entry