找到约 15 条结果

问题 全基因组关联分析GWAS教程中的beagle填充问题

...galArgumentException: Sample Line68 has an inconsistent number of alleles. The first genotype is diploid, but the genotype at position Chr2:35902177 is haploid

文章 shuf 随机抽取txt文件的行到新文件 linux命令

...数]... 或者:  shuf -i LO-HI [选项]...Write a random permutation of the input lines to standard output. 如果没有指定文件,或者文件为"-",则从标准输入读取。 必选参数对长短选项同时适用。  -e, --echo                treat each ARG as an input line...

问题 老师您好,我在做物种间基因组共线性分析时没有得到.anchors文件并且出现list index out of range报错,请问一下是什么原因?

所选的基因组是线粒体基因组测序的结果。

问题 在docker中用cleandata与宿主序列比对得到sam文件这一步出现错误,

...tie.map.log  # Filter and output as SAM(or BAM),-f 4 to Remove reads from the alignment,  # if you want to keep the -F on the alignment (12 for both ends), add -b before the output to output in BAM format samtools view --threads 6 \ -f 12 -h \ $workdir/1.host_del_and_qc/cleandata/${i}.bowtie.sam \...

文章 BUSCO评估基因组组装和完整性详解

...O评估结果文件。 格式如下: # BUSCO version is: 5.3.2# The lineage dataset is: eukaryota_odb10 (Creation date: 2020-09-10, number of genomes: 70, number of BUSCOs: 255)# Summarized benchmarking in BUSCO notation for file /share/nas5/wangq/project/genome_ass...

问题 染色体挂载,bwa比对问题

您好,根据基基因组装与注释课程,我在做自己物种染色体挂载时,使用bwa aln 和bwa sampe 参数,提示 Unmatched SAI magic. Please re-run `aln' with the same version of bwa. 重装软件也不管用,这种情况怎么解决?

问题 老师好!我使用JCVI按照官网流程处理两个物种,每次都提示ValueError: A total of 0 anchor was found. Aborted.麻烦您帮我看看,谢谢!

文章 再挖TCGA,发篇SCI

...金思路。接下来以一篇2018年7月30日发表在Cellular Physiology and Biochemistry(IF=5.5)上题为Genome-Wide Analysis of Prognostic lncRNAs, miRNAs, and mRNAs Forming a Competing Endogenous RNA Network in Hepatocellular Carcinoma 的文章为例,详细介绍下该类文章的写作...

问题 变异结果统计与绘图时,最后一步出现Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : line 17 did not have 34 elements的问题,查看文件发现并没有缺失的部分,怎么解决

...ellow,red" -t "SNP Density"这条代码后出现Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec,  :   line 17 did not have 34 elements的问题, 查看了第17行,34个元素都在,并没有缺失

文章 安装R包ggpubr 报错

...ect.hpp:186:5: warning: control reaches end of non-void function [-Wreturn-type]     }     ^make: *** [cpp11.o] Error 1ERROR: compilation failed for package ‘tidyr’ 这个是因为R更新4版本以后的报错,而你的R版本是3.* 所以报错,建议安装老板版本...

文章 8.28接收-毛竹bHLH基因家族分析文章解读及原文献下载

...其中16个亚家族中(except 1,7, 13, 14, and 20)。 图注: a: The 137 PebHLH genes are clustered into 16 subfamilies. b: Comparison of bHLH family members from Moso bamboo, rice and Arabidopsis. Different colors represent the different species, and the percentage of each subfamily in e...

问题 老师您好,我在做物种间共线性分析绘图时,显示list index out of range,我检查了准备的bed文件和cds文件没问题。请问这种情况怎么处理?

问题 在运行table_annovar.pl脚本的时候报错

...sample_INS_DEL.vcf.gz和unknow相关文件见附件 归档.zip NOTICE: the --polish argument is set ON automatically (use --nopolish to change this behavior) NOTICE: Running with system command <convert2annovar.pl  -includeinfo -allsample -withfreq -format vcf4 all_sample_INS_DEL.vcf.gz &...

问题 非模式物种建立注释库问题

...rieve the full column specification for this data. ℹ Specify the column types or set `show_col_types = FALSE` to quiet this message. Error in file(file, "rt") : cannot open the connection Calls: read.table -> file In addition: Warning message: In file(file, "rt") :   cannot open file 'NA...

文章 clusterProfiler非模式物种 KEGG与GO富集分析

...Glutathione metabolism ko00510 N-Glycan biosynthesis ko00513 Various types of N-glycan biosynthesis ko00512 Mucin type O-glycan biosynthesis 5.富集分析 利用clusterProfiler中的enricher这个通用函数进行富集分析: library(clusterProfiler)ko2name <- read.delim('ko2...