https://github.com/TRON-Bioinformatics/seq2HLA
Raw File
Tip revision: 15dc4ef3962f55f3e7060336991f22dbd3d9c5c6 authored by Patrick Sorn on 21 August 2023, 14:26:55 UTC
Added Python 3 support; Updated dependencies
Tip revision: 15dc4ef
command.R
args <- commandArgs(trailingOnly = TRUE)

n <- as.numeric(args[1])/3
offset <- 0
for(i in 1:n)
{
	x <- unlist(strsplit(args[i*2+1+offset], split=","))
      	x <- as.numeric(x)
      	print(args[i*2 + 2 + offset])
     	paril <- 1 - pnorm(as.numeric(args[i*2 + offset]), mean(x), sd(x))
	poutlier <- pbinom(0, length(x), paril)
     	print(1 - poutlier)
      	offset <- offset + 1
}
back to top