MAP*_*APK 3 linux bash r system
我有一个像这样运行的 bash 命令:
esearch -db protein -query "AVA17449.1" | elink -target nuccore | efetch -format ft
Run Code Online (Sandbox Code Playgroud)
但我想在 R 中这样做(这不起作用)
output <- system("esearch -db protein -query "AVA17449.1" | elink -target nuccore | efetch -format ft")
Run Code Online (Sandbox Code Playgroud)
在 R 中调用此命令的正确方法是什么?
PS esearch 可以使用以下命令安装
cd ~
/bin/bash
perl -MNet::FTP -e \
'$ftp = new Net::FTP("ftp.ncbi.nlm.nih.gov", Passive => 1);
$ftp->login; $ftp->binary;
$ftp->get("/entrez/entrezdirect/edirect.tar.gz");'
gunzip -c edirect.tar.gz | tar xf -
rm edirect.tar.gz
builtin exit
export PATH=$PATH:$HOME/edirect >& /dev/null || setenv PATH "${PATH}:$HOME/edirect"
./edirect/setup.sh
Run Code Online (Sandbox Code Playgroud)
要么使用单引号:
cmd <- paste(c('esearch -db protein -query "AVA17449.1"',
'elink -target nuccore',
'efetch -format ft'),
collapse=" | ")
Run Code Online (Sandbox Code Playgroud)
(这些paste(...,collapse=...)东西不是必需的,但可能会让你的代码更具可读性......)
或者用反斜杠保护双引号:... -query \"AVA17449.\" ...
单引号可能更具可读性,但反斜杠可以在更复杂的情况下工作,在这种情况下,您需要在字符串中使用两种引号......
| 归档时间: |
|
| 查看次数: |
1925 次 |
| 最近记录: |