mirror of
https://github.com/dchakro/shared_Rscripts.git
synced 2026-05-17 09:12:16 -07:00
52554283f1
commit fe050d47c2dc9e61a13d4a362abf56dabfbecbd7
Author: Deepankar Chakroborty <deepankar.chakroborty@utu.fi>
Date: Fri Jul 31 12:57:11 2020 +0300
Include usage instruction inside function body
Include usage instruction inside function body, so on typing just the function name to display the R code, the instructions will appear.
Also stated things clearly in the T&C
commit fedbca4a7aa8fb436702accdf1f9b5ec7f066d60
Author: Deepankar Chakroborty <deepankar.chakroborty@utu.fi>
Date: Thu Jul 30 12:50:20 2020 +0300
Updates
1. Improve readability
2. Manage differing lengths of breaks and labels.
commit fb32fb9173c66c8862adb0ad41d82db9d128777b
Author: Deepankar Chakroborty <deepankar.chakroborty@utu.fi>
Date: Thu Jul 30 12:43:20 2020 +0300
skip.steps works as expected
skip.steps = 1, now skips 1 observation.
commit 00a8bbcf04731c010e0a5989a805274609560188
Author: Deepankar Chakroborty <deepankar.chakroborty@utu.fi>
Date: Thu Jul 30 12:32:16 2020 +0300
add script to calculate breaks for axes in ggplot2
calculates breaks and labels for axes in ggplot2 with user defined gaps
127 lines
4.2 KiB
R
127 lines
4.2 KiB
R
# Installing dependencies
|
|
dependencies <- c("stringi", "doParallel")
|
|
missing_packages <- dependencies[!(dependencies %in% installed.packages()[, "Package"])]
|
|
if(length(missing_packages)) install.packages(missing_packages)
|
|
rm(missing_packages,dependencies)
|
|
|
|
IsolateCanonicalVariant <- function (AAchangeAnnotations){
|
|
# #<---------------------------->
|
|
# # You must include this section when:
|
|
# # Distributing, Using and/or Modifying this code.
|
|
# # Please read and abide by the terms of the included LICENSE.
|
|
# # Copyright 2020, Deepankar Chakroborty, All rights reserved.
|
|
# #
|
|
# # Author : Deepankar Chakroborty (https://gitlab.utu.fi/deecha)
|
|
# # Report issues: https://gitlab.utu.fi/deecha/shared_scripts/-/issues
|
|
# # License: https://gitlab.utu.fi/deecha/shared_scripts/-/blob/master/LICENSE
|
|
# #<---------------------------->
|
|
|
|
|
|
# # PURPOSE:
|
|
# # From a given vector of annotations for a particular DNA change
|
|
# # this function selects the canonical variant (if present)
|
|
# # by cross referencing the MANE Select and RefSeq Select sets.
|
|
|
|
# # LOGIC FLOW:
|
|
# # - If there is only one annotation; that is selected
|
|
# # - If canonical transcript is not found in MANE Select + RefSeq select
|
|
# # or a matching transcript ID is not found in the annotation then;
|
|
# # The mutation with to the the highest position (residue number) is selected.
|
|
# # - If a match for canonical isoform is found then;
|
|
# # that particular mutation is selected
|
|
|
|
# importing resources
|
|
library(doParallel)
|
|
refseq <- readRDS(url("https://gitlab.utu.fi/deecha/shared_scripts/-/raw/master/asset/RefSeqSelect_Gene_Transcript.RDS"),"rb")
|
|
source("https://gitlab.utu.fi/deecha/shared_scripts/-/raw/master/MutSiteFind.R")
|
|
|
|
# initializing cluster
|
|
myCluster <- makeCluster(parallel::detectCores(),
|
|
type = "FORK",
|
|
useXDR=F,
|
|
.combine=cbind)
|
|
registerDoParallel(myCluster)
|
|
print(myCluster)
|
|
|
|
file.create("log.txt")
|
|
message(paste0("Logging processed mutations at: ",getwd(),"/log.txt"))
|
|
|
|
# computation
|
|
results <- foreach(MutInfo = AAchangeAnnotations,.combine = c) %dopar% {
|
|
GENE <- can.isoform <- l <- l2 <- NA
|
|
l=unique(unlist(
|
|
stringi::stri_split_fixed(
|
|
str = MutInfo,
|
|
pattern = ","),
|
|
use.names = F,
|
|
recursive = F))
|
|
|
|
if(length(l)==1){
|
|
l2 <- unlist(stringi::stri_split_fixed(
|
|
str = l,
|
|
pattern = ":"),
|
|
use.names = F,
|
|
recursive = F)
|
|
MUTATION <- stringi::stri_replace_first_fixed(
|
|
str = l2[stringi::stri_detect_regex(
|
|
str = l2,
|
|
use.names = F,
|
|
pattern = "^p\\.")],
|
|
pattern = "p.",
|
|
replacement = "")
|
|
} else {
|
|
GENE <- stringi::stri_sub(
|
|
str = l[1],
|
|
from = 1,
|
|
to = stringi::stri_locate_first_fixed(
|
|
str = l[1],
|
|
pattern = ":")[,"end"]-1)
|
|
|
|
can.isoform <- refseq$transcript_accession[refseq$gene_id==GENE]
|
|
|
|
if(length(can.isoform)==0 | !any(stringi::stri_detect_fixed(str = l,pattern = can.isoform))){
|
|
# Canonical isoform not found in RefSeq, or Match for Canonical isoform not found in AAchangeAnnotations
|
|
l2 <- unlist(
|
|
stringi::stri_split_fixed(
|
|
str = l,
|
|
pattern = ":"),
|
|
use.names = F,
|
|
recursive = F)
|
|
|
|
l3 <- stringi::stri_replace_all_fixed(
|
|
str = l2[stringi::stri_startswith_fixed(
|
|
str = l2,
|
|
"p.")],
|
|
pattern = "p.",
|
|
replacement = "")
|
|
|
|
MUTATION <- l3[which.max(MutSiteFind(l3))]
|
|
} else {
|
|
# Canonical isoform found
|
|
l=l[grep(can.isoform,l)]
|
|
|
|
l2 <- unlist(
|
|
stringi::stri_split_fixed(
|
|
str = l,
|
|
pattern = ":"),
|
|
use.names = F,
|
|
recursive = F)
|
|
|
|
MUTATION <- stringi::stri_replace_first_fixed(
|
|
str = l2[stringi::stri_detect_regex(
|
|
str = l2,
|
|
use.names = F,
|
|
pattern = "^p\\.")],
|
|
pattern = "p.",
|
|
replacement = "")
|
|
}
|
|
}
|
|
if(length(MUTATION)==0){
|
|
MUTATION <- NA
|
|
}
|
|
cat(paste(MUTATION,"\n"),file="log.txt", append=T)
|
|
return(MUTATION)
|
|
}
|
|
stopCluster(myCluster)
|
|
return(results)
|
|
}
|