...are()`:
! Must group by variables found in `.data`.
• Column `cluster` is not found.
Backtrace:
▆
1. ├─obj.markers %>% group_by(cluster) %>% ...
2. ├─dplyr::top_n(., n = opt$top_n, wt = avg_log2FC)
3. │ └─dplyr::filter(...)
4. ├─dplyr::group_by(., cl...
...,在进行ka/ks计算的时候,发现输入指令/biosoft/KaKs_Calculator2.0/src/AXTConvertor dup_gene_paired1.phy dup_gene_paired1.axt,总是出现错误,zsh: no such file or directory: /biosoft/KaKs_Calculator2.0/src/AXTConvertor,应该怎么处理?我用的是咱们网站上下载的...
...ocker: Cannot connect to the Docker daemon at unix:///var/run/docker.sock. Is the docker daemon running?.
See 'docker run --help'.
(base) [libowen@manage gwas]$
...GenomicsDBImport - GenomicsDBImport cannot use multiple VCF reader threads for initialization when the number of intervals is greater than 1. Falling back to serial VCF reader initialization.
15:54:41.529 WARN IntelInflater - Zero Bytes Written : 0
15:54:41.529 WARN IntelInflater - Zero Bytes ...
for species in Anser_brachyrhynchus Equus_caballus Homo_sapiens Mus_musculus Oryctolagus_cuniculus Ovis_aries_rambouillet Sus_scrofa.Sscrofa;do
genome=$datadir/$species.genome.fa.gz
gff=$datadir/$species.gff.gz
#保留蛋白编码基因
agat_sp_filter_feature_by_attribute_value...
...MC10110842.
数据:
单细胞:Single-cell RNA-seq (scRNA-seq) data for 5 patients with EAC (GSE173682)
空间转录组:spatial transcriptomic data for one patient with EC (GSM6177623 in GSE203612) were downloaded from Gene Expression Omnibus (GEO).
TCGA数据:TCGA-UCEC project from TC...
...器,没有用docker镜像。出现
Error in library("getopt") : there is no package called ‘getopt’
Execution halted
检查R里已经安装了getopt,调用这个R包也是正常的,但是perl $scriptdir/vcf_stat.pl -vcf snp.unknown_multianno.vcf -type snp就会报错。