找到约 15 条结果

问题 转录组上游分析之定量比对Hisat2出现内部异常错误

转录组上游分析(Linux)之定量比对Hisat2出现内部异常错误,不知道哪里错了?如图下: 下面为上图的红圈error信息: Extra parameter(s) specified: "?/home/vip6/database/genome/Ensembl/Homo_sapiens/GRCh38_release104/Homo_sapiens.GRCh38_release104.genome" No...

文章 安装R包出现链接库缺失的问题

 unable to load shared object '/home/wanghj/R/x86_64-pc-linux-gnu-library/4.2/00LOCK-igraph/00new/igraph/libs/igraph.so':  libglpk.so.36: cannot open shared object file: No such file or directory 如上面的报错,这个时候最简单的办法是通过locate查看是否存在链接库(linux里)...

问题 vcftools过滤除了等位基因选项外,其他参数通通没用

...eles 2 |gzip - > all.clean.vcf.gz 结果是After filtering, kept 0 out of a possible 118111 Sites(中间很多Warning),于是把参数全部拆开一次只过滤一个参数,发现除了--min-alleles 2  --max-alleles 2,可以过滤成功:kept 106816 out of a possible 118111 Sites之...

问题 R语言,KNN补值函数impute.knn函数

...         rng.seed = 362436069)  {   rng.state <- NULL   if (exists(".Random.seed")) {     rng.state <- .Random.seed   }   set.seed(rng.seed)   x <- data   p <- nrow(x)   col.nas <- drop(rep(1, p) %*% is.na(x))  #matrix multiplication   #count NA number in eac...

文章 rocks7 集群安装笔记

...html 这里也提供度盘资源: 链接:https://pan.baidu.com/s/1nLjxZ0VnmbR_QUEfFS0rAA 提取码:87ds 2、镜像部署 镜像悉数下载完成并检查文件完整性无误后,我们来对镜像进行部署;下载的镜像文件存放:D:/rocks/roll-server/rolls 目录下; ...

问题 转录组kegg富集,deg emsembl id不存在

筛选出的deg的基因emsembl id无法进行转成entrezid,说不存在这些id Error in .testForValidKeys(x, keys, keytype, fks) :    None of the keys entered are valid keys for 'ENSEMBL'. Please use the keys method to see a listing of valid arguments. Calls: bitr Execution halted

文章 典型相关性分析CCA与冗余分析RDA

...(detrended correspondence analysis) 分析,看分析结果中Lengths of gradient 的第一轴的大小。如果大于4.0,就应该选CCA;如果3.0-4.0 之间,选RDA 和CCA均可;如果小于3.0,RDA 的结果要好于CCA。目前CCA要比RDA用得更普遍,有两个原因:大部...

文章 network.r绘制网络图并输出最大子图节点及网络图

...: Rscript network.r -h usage: network.r [-h] -e filepath -v filepath -f from -t to -w weight -g group [--isolated] [-l style] [--label.node] [--sub_label.node] [-o path] [-n prefix] [-H number] [-W number] network plot:https://www.omicsclass.com/article/1531 opti...

文章 cactus下载

...组的时候需要用到Progressive Cactus 参考文章:Dense sampling of bird diversity increases power of comparative genomic 下载: 软件github:https://github.com/ComparativeGenomicsToolkit/cactus 选择的是Pre-compiled Binaries Linux Tarball这个版本: wget https://github.com...

文章 BWA比对及Samtools提取mapped reads

...是序列末尾的READ1 flag和READ2 flag,flag有3类: 1:Only READ1 is set. 2:Only READ2 is set. O:Either both READ1 and READ2 are set;or neither is set. 对于PE测序reads,同时指定-1 R1.fq-2 R2.fq-s singleton.fq 时,samtools会将flag1和flag2配对的序列分别输出...

问题 是没有gene-family了嘛

...s → docker scout quickview omicsclass/gene-family:v1.0.1 Error response from daemon: pull access denied for omicsclass/gene-family, repository does not exist or may require 'docker login': denied: requested access to the resource is denied PS C:\Users\Administrator> docker search omicsclass ...

问题 LDdecay

...放了所有样品的集合,但是跑出来是如下这样 the Number of subPop samples[found in VCF] is 0 sub Group Population szie is too small, SubGroup sample size: 0 ##begin pair-wise R^2 cal. after filter Remain SNP Number :     0 #% number bin is        1 1%......-->100%........

问题 tassel新报错,求老师们看一下[Thread-8] ERROR net.maizegenetics.plugindef.AbstractPlugin - Index 0 out of bounds for length 0

... generated? (BLUE = best linear unbiased estimate) (Default: false) -saveToFile <true | false> : Should the results be saved to a file rather than stored in memory? It true, the results will be written to a file as each SNP is analyzed in order to reduce memory requirementsand the results wil...

问题 GWAS分析绘制进化树遇到下面问题是什么意思,请问怎么解决?

...istered S3 methods overwritten by 'treeio':   method              from   MRCA.phylo          tidytree   MRCA.treedata       tidytree   Nnode.treedata      tidytree   Ntip.treedata       tidytree   ancestor.phylo      tidytree   ancestor.treedata   tidytree ...

问题 T2T基因组调研图中jellyfish -h和genomescope.R -m 都是设置的10000000,设置的依据是什么,不同作物有差异,怎么选择参数呢

...-o illu.histo illu# histo输出k-mer频率直方图数据,统计出现了x(x=1,2,3…)次的kmer的数目y# -v 输出详细信息# -h 最高计数值genomescope.R -i illu.histo -o genomescope_out -p 2 -k 21 -m 10000000# -p PLOIDY,倍性# -m 最大 Kmer 覆盖阈值