小编Ano*_*non的帖子

我正在使用 WSL2。\n在 R 中，我尝试安装软件包。例如：

> install.packages("tidyverse")\n

但我收到Failed to create bus connection错误。

追溯：

Failed to create bus connection: No such file or directory\nWarning in system("timedatectl", intern = TRUE) :\n  running command 'timedatectl' had status 1\n** help\n*** installing help indices\n*** copying figures\n** building package indices\n** installing vignettes\n** testing if installed package can be loaded from temporary location\nFailed to create bus connection: No such file or directory\nWarning in system("timedatectl", intern = TRUE) :\n  running command 'timedatectl' had …

在代码的第二行中，我尝试按*.out.tab列合并所有文件。代码的第三行提取第一列和每4 个后续列（第 4、8、12、16...），这意味着每个文件的每第四列。

如果没有 for 循环，就会像......

paste 1.out.tab 2.out.tab 3.out.tab 4.out.tab \
awk '{for(i=1;i<=NF;i+=4){printf "%s ",$i;} print ""}' | \
tail -n +5 > tmpfile
cat tmpfile | sed "s/^ENSG*//" >gene_count.txt

但是，现在我想使用 for 循环来合并所有文件。

for f in `./alignments/repaired_reads/*ReadsPerGene.out.tab | sed 's/.ReadsPerGene.out.tab//'`;
paste "$f"\.out.tab | \
awk '{for(i=1;i<=NF;i+=4){printf "%s ",$i;} print ""}' | \
tail -n +5 > tmpfile
cat tmpfile | sed "s/^ENSG*://" > gene_count.txt

输入示例：

head ./alignments/repaired_reads/SRR9200814ReadsPerGene.out.tab

N_unmapped  18517   18517   18517
N_multimapping  1620 …

我想使用hic2coolfile1.hic将文件转换为用 Python 编写的.cool格式。我使用命令行转换文件：

hic2cool convert file1.hic file1.cool -r 10000

追溯：

Traceback (most recent call last):
  File "/home/melchua/.local/bin/hic2cool", line 5, in <module>
    from hic2cool.__main__ import main
  File "/home/melchua/.local/lib/python3.9/site-packages/hic2cool/__init__.py", line 2, in <module>
    from .hic2cool_utils import (
  File "/home/melchua/.local/lib/python3.9/site-packages/hic2cool/hic2cool_utils.py", line 27, in <module>
    import cooler
  File "/home/melchua/.local/lib/python3.9/site-packages/cooler/__init__.py", line 14, in <module>
    from .api import Cooler, annotate
  File "/home/melchua/.local/lib/python3.9/site-packages/cooler/api.py", line 12, in <module>
    from .core import (get, region_to_offset, region_to_extent, RangeSelector1D,
  File "/home/melchua/.local/lib/python3.9/site-packages/cooler/core.py", line 3, in <module>
    from pandas.api.types …

如何使用 for 循环合并下划线之前具有匹配字符/数字的两个文件？我的目录中有很多文件

输入：

SRR9200887_1.fastq
SRR9200887_2.fastq
SRR9200888_1.fastq
SRR9200888_2.fastq
SRR9200889_1.fastq
SRR9200889_2.fastq

预期输出：

SRR9200887.fastq
SRR9200888.fastq
SRR9200889.fastq

我的尝试：

for l in $(ls *.fastq | cut -d_ -f1 | sort |uniq); do cat ${l}*.fastq