找到约 15 条结果

问题 RNAseq有参转录组数据自主分析课程 ,运行建立基因组索引命令时,报错,无法获取gene_length.txt文件,请问是什么原因?如何处理呢?

...genome.fasta Toona_sinensis.gene.gff 报错信息如下:get gene length and gene.bed from gtf: RUN CMD: python /work/TS/scripts/get_gene_length_from_gtf.py -g Toona_sinensis.gene.gtf -p gene_length Traceback (most recent call last):   File "/work/TS/scripts/get_gene_length_from_gtf.py", line ...

文章 lefse 运行报错python3中

...s no attribute 'rownames' 解决办法,降低版本: rpy2 (3.5.10) and numpy (1.24.3)

问题 请教老师,Gemma分析时的报错问题?

[root@iZ2e]# ./gemma-0.98.5-linux-static-AMD64 -bfile c -k output/result1.sXX.txt -lmm 1 -p p.txt -c c.txt GEMMA 0.98.5 (2021-08-25) by Xiang Zhou, Pjotr Prins and team (C) 2012-2021 ERROR: Enforce failed for file_kin output/result1.sXX.txt: open file in src/param.cpp at line 934 in CheckParam

文章 重测序群体遗传进化分析之—进化树构建

...多不同的模型计算初始简约树的逻辑概率,并产生Akaike information criterion (AIC),* corrected Akaike information criterion* (AICc), and* the Bayesian information criterion* (BIC)三个结果标准值,通常ModelFinder选择BIC分数最低的模型(当然也可以指定AIC和AI...

问题 WGCNA-MAD

...at is in the data set dim(datExpr) ##################check missing value and filter ################### datExpr = t(datExpr[order(apply(datExpr,1,mad), decreasing = T)[1:5000],]) dim(datExpr) 根据上面的代码,我想把1.8w个基因缩减下,结果在进行MAD筛选前5k个基因的...

文章 4个样品重测序发PlosOne

...高深度重测序,四个品种中三个高抗白粉病( Edisto, MR-1, and PMR5: resistant),一个易感白粉病(SCNU1154: susceptible;),得到测序reads与参考基因组比对,查找SNP,Indel,SV等变异信息,并对变异信息进行质控过滤,变异信息在染色体...

问题 我的Conda 根本安装不成功软件,比如UMI tools ,怎么办

...es.DistributionNotFound: The 'umi-tools==0.2.3' distribution was not found and is required by the application

问题 ANNOVAR变异注释报错:Erro: invalid record found in exonic_variant_function file (exonic format error)

...use --nopolish to change this behavior) NOTICE: Running with system command <convert2annovar.pl  -includeinfo -allsample -withfreq -format vcf4 /work/bdorsalis/4.snp_indel/var_qc/all.clean.snp.vcf.gz > /work/bdorsalis/5.var_ann/test/snp.avinput> NOTICE: Finished reading 840 lines from ...

问题 请问一下做WGCNA时,出现这种问题如何解决呢?用的是课程的数据

traitColors = numbers2colors(datTraits, signed = FALSE) Warning message: In numbers2colors(datTraits, signed = FALSE) : (some columns in) 'x' are constant. Their color will be the color of NA. > # Plot the sample dendrogram and the colors underneath.

文章 VG 泛基因组call 变异

#this script report the commands used to call variants from the 16 Illumina WGS barn swallow individuals.#for steps 1-3, this github solution was followed: https://github.com/vgteam/vg/issues/3411#for steps 4-9, the section "Calling the graph by first splitting into components" of this tutorial was ...

文章 R 安装github包报错:

Using github PAT from envvar GITHUB_PAT. Use `gitcreds::gitcreds_set()` and unset GITHUB_PAT in .Renviron (or elsewhere) if you want to use the more secure git credential store instead. Error: Failed to install 'unknown package' from GitHub:   HTTP error 401.   Bad credentials   Rate limit ...

文章 LENS:人的基因相互作用分析在线工具

... interactome analysis of human genes. LENS takes one or two lists of genes and determines the network of protein-protein interactions that connects them. LENS computes several statistics including associations of the genes in the interactome to diseases, drugs, pathways, and GWASs, as well as, stati...

问题 我用conda 根本安装不了软件,非常慢,比如UMI-tools

...es.DistributionNotFound: The 'umi-tools==0.2.3' distribution was not found and is required by the application

文章 基因组版本 Primary 和 Alternate的区别

...mpanies a primary assembly. It is not useful by itself as it is fragmented and incomplete.https://lh3.github.io/2021/04/17/concepts-in-phased-assemblies

文章 blast序列比对中E值是什么意思?

...匹配了 E-value - this tells you the probability of getting this score and homology by chance (from random sequence), that is, the lower (the higher the number behind e-), the more reliable the result is. As I have mentioned above, it's not only about homology, but also about the length, because ...