R -- Rscript传递参数

文章目录

getopt包

getopt(spec = NULL, opt = NULL, command = get_Rscript_filename(),

usage = FALSE, debug = FALSE)

  • spec:
    内容由一个4-5列的matrix构成,matrix可以直接申明也可以如下申明 matrix(x,ncol=4,byrow=TRUE)
    每一列的构成要素如下:
    Column 1: the long flag name. A multi-character string.
    Column 2: short flag alias of Column 1. A single-character string.
    Column 3: Argument mask of the flag. An integer. Possible values: 0=no argument, 1=required argument, 2=optional argument.
    Column 4: Data type to which the flag's argument shall be cast using storage.mode.
    Column 5 (optional): A brief description of the purpose of the option.
R 复制代码
#!/path/to/Rscript
library('getopt')

# get options, using the spec as defined by the enclosed list.
# we read the options from the default: commandArgs(TRUE).
spec = matrix(c(
	'verbose', 'v', 2, "integer",
	'help' , 'h', 0, "logical",
	'count' , 'c', 1, "integer",
	'mean' , 'm', 1, "double",
	'sd' , 's', 1, "double"
	), byrow=TRUE, ncol=4)

# 脚本运行到这里时,外部传递的参数会赋值给opt
opt = getopt(spec)

# 这里设置帮助信息,并且如果接收到了帮助信息指令则打印帮助信息并且终止后续命令的执行
# if help was asked for print a friendly message
# and exit with a non-zero error code
if ( !is.null(opt$help) ) {
	cat(getopt(spec, usage=TRUE))
	q(status=1)
}

# 默认参数设置部分
# set some reasonable defaults for the options that are needed,
# but were not specified.
if ( is.null(opt$mean ) ) { opt$mean = 0 }
if ( is.null(opt$sd ) ) { opt$sd = 1 }
if ( is.null(opt$count ) ) { opt$count = 10 }
if ( is.null(opt$verbose ) ) { opt$verbose = FALSE }

# 设置脚本运行情况的信息反馈部分
# print some progress messages to stderr, if requested.
if ( opt$verbose ) { write("writing...",stderr()) }

# do some operation based on user input.
cat(paste(rnorm(opt$count,mean=opt$mean,sd=opt$sd),collapse="\n"))
cat("\n")
# signal success and exit.
# q(status=0)
R 复制代码
# Start ========================================================================
#!/path/to/Rscript


# 脚本参数设置和接收部分
library(getopt)

spec = matrix(c(
    'input',  'i', 1, "character",
    'output', 'o', 1, "character",
    'help',   'h', 0, "logical"), byrow=TRUE, ncol=4)
    
args <- getopt(spec)

if(!is.null(args$help) || is.null(args$input) || is.null(args$output)){
    cat(paste(getopt(spec, usage = T), collapse = "\n"))
    cat("   Description: This script is aimed to calculate the ratio in malignant cell! Example:\n
            Rscript cal_ratio.R -i BRCA.rds -o BRCA_output\n
    Written by XqChai, 2023-07-20. \n")
    q(status=1)
}

if (is.null(args$output)) {args$species = "out"}

# 操作系统或者数据的部分
library(tidyverse)
library(reshape2)
library(stringi)

#读取rds文件
input <- args$input 
output <- args$output
sce <- readRDS(input)

optparse包

make_option(
  opt_str, # 参数名,格式:--parameter_name
  action = NULL, # 描述optparse在遇到选项"store"、"store-true"、"store-false"或"callback"时应采取的操作的字符串。
 	  # 如果在命令字符串上找到选项,则"store"操作表示optparse应存储指定的以下值。
	  # 如果找到选项,"store_true"将存储true,如果找到选项,"store_false"将存储false。
	  # "callback"存储回调参数中指定的函数生成的返回值。如果callback不为空,则默认值为"callback"else"store"。
  type = NULL, # 参数数据类型
  dest = NULL,
  default = NULL, #设置参数的默认值
  help = "", # 该参数的描述信息
  metavar = NULL,
  callback = NULL,
  callback_args = NULL
)

option_list = list(make_option(), make_option(), make_option()) 
R 复制代码
#!/path/to/Rscript

# 参数接受和设置部分
library(optparse)

option_list = list(
    make_option(c("-f", "--file"), type="character", default=NULL, 
                help="rpkm/fpkm matrix file", metavar="input"),
    make_option(c("-o", "--out"), type="character", default="out.tpm.xls", 
                help="output file name [default= %default]", metavar="output")
); 

opt_parser = OptionParser(option_list = option_list, epilogue = "Written by Chai X. Q. 2019-04-22.")
opt = parse_args(opt_parser);

if (is.null(opt$file)){
    print_help(opt_parser)
    stop("At least one argument must be supplied (input file).", call.=FALSE)
}

# 操作系统或者数据部分
#program
input <- opt$file
output <- opt$out
exp <- covert(input)

argparse包

该部分内容完全来自这位作者以及他的文字:

https://zhuanlan.zhihu.com/p/212997106

然后本人通常使用第一种(getopt)方法,第二种方法也使用过,但是这一部分内容仅仅是记录和比较。

argparse包分为三个步骤:

第一步是使用ArgumentParser()创建一个参数解析对象,

第二步是使用add_argument函数为前面的对象增加参数,

第三步是通过parse_args()把解析对象赋值并在后续程序中调用。

R 复制代码
#! /usr/bin/env Rscript
library(argparse)

# 创建参数解析对象
parser <- ArgumentParser()

# 设置参数
# 设置第一个参数verbose,缩写为v,其作用是告诉脚本是否打印完整的计算过程,其缺省值为TRUE
parser$add_argument("-v", "--verbose", action="store_true", default=TRUE,
        help="Print extra output [default]")

# 设置第二个参数quietly,缩写为q,其作用是修改verbose参数,当调用改参数时,verbose被修改为FALSE,从而导致不再打印计算过程
parser$add_argument("-q", "--quietly", action="store_false", 
        dest="verbose", help="Print little output")

# 设置第三个参数count,缩写为c,这是一个整数参数,缺省值5,在后续的代码中被用作确定输出随机数的个数
parser$add_argument("-c", "--count", type="integer", default=5, 
        help="Number of random normals to generate [default %(default)s]",
        metavar="number")

# 设置第四个参数generator,无缩写,用于确定调用何种随机分布,缺省值rnorm对应于正态分布
parser$add_argument("--generator", default="rnorm", 
        help = "Function to generate random deviates [default \"%(default)s\"]")

# 设置第五个参数mean,无缩写,浮点数,用于确定正态分布的均值,缺省值为0
parser$add_argument("--mean", default=0, type="double",
        help="Mean if generator == \"rnorm\" [default %(default)s]")

# 设置第五个参数sd,无缩写,浮点数,用于确定正态分布的标准差,缺省值为1
parser$add_argument("--sd", default=1, type="double",
        metavar="standard deviation",
        help="Standard deviation if generator == \"rnorm\" [default %(default)s]")

# 调用解析器,此时args就被赋值为命令行参数输入的相应值
args <- parser$parse_args()

# 根据verbose确定是否打印计算过程
if ( args$verbose ) { 
        write("writing some verbose output to standard error...\n", stderr()) 
}

# 根据其他参数,确定输出随机数的类型与次数
if( args$generator == "rnorm") {
        cat(paste(rnorm(args$count, mean=args$mean, sd=args$sd), collapse="\n"))
} else {
        cat(paste(do.call(args$generator, list(args$count)), collapse="\n"))
}
cat("\n")
相关推荐
satan–013 小时前
R语言的下载、安装及环境配置(Rstudio&VSCode)
开发语言·windows·vscode·r语言
PhyliciaFelicia2 天前
基于R语言机器学习遥感数据处理与模型空间预测
开发语言·深度学习·随机森林·机器学习·数据分析·r语言
长安不及十里3 天前
PaddleDetection 自定义训练目标检测
人工智能·目标检测·r语言
WangYan20224 天前
ChatGPT+R语言强强联合,数据分析不再难!回归与混合效应模型、多元统计分析、结构方程模型(SEM)(lavaan)、Meta分析、贝叶斯回归等应用
chatgpt·数据分析·r语言·结构方程模型·多元统计分析·回归与混合效应模型
hakesashou5 天前
python和r语言的区别是什么
开发语言·python·r语言
青春不流名5 天前
centos72009源码编译R语言
开发语言·r语言
SofterICer5 天前
pysim-4-1.1.17 eUICC ISD-R commands
开发语言·r语言
让学习成为一种生活方式5 天前
基因共表达分析-R-脚本04
java·windows·r语言
卡卡_R-Python5 天前
简单线性回归分析-基于R语言
算法·r语言·线性回归
纪伊路上盛名在6 天前
如何初步部署自己的服务器,达到生信分析的及格线2(待更新)
linux·运维·服务器·python·学习·r语言