找到约 15 条结果

问题 mummer比对syri进行SV检测报错

...s do not match. Reading Coords - ERROR - Unequal number of chromosomes in the genomes. Exiting out.anno

问题 GSVA Error: useNames = NA is defunct. Instead, specify either useNames = TRUE or useNames = FALSE.

...rixStats 1.2 that released a couple of days ago. Can you downgrade to 1.1 and try?remotes::install_version("matrixStats", version="1.1.0") # restart your session and run previous scripts

文章 合并不同sample的vcf文件,通过bcftools

通过GATK calling出来的SNP如果使用UnifiedGenotype获得的SNP文件是分sample的,但是如果使用vcftools或者ANGSD则需要Vcf文件是multi-sample的,这里就需要我们将不同samples的文件进行合并,可以通过vcftools的perl模块进行,但是这种方式对perl的...

文章 跟着文献学做图 | 用蛋白质序列做主成分分析(PCA)

... for each locus   @ploidy: ploidy of each individual  (range: 1-1)   @type:  codom   @call: alignment2genind(x = mase.res) // Optional content   @other: a list containing: com 四、对genind对象执行主成分分析 genind对象中包含的表格可以提交给主成分分析(PCA) ...

问题 使用plink做pca分析时报错

...' on line 30616892 of .vcf file.(This is disallowed for humans.  Check if the problem is with your data, or ifyou forgot to define a different chromosome set with e.g. --chr-set.) 使用的命令如下: plink --vcf  $workdir/2022_clean.vcf.gz --pca 10 --out  plink_pca   \    --allow-extr...

问题 Perl GD模块官网例子报错

...e(255,0,0);       $blue = $im->colorAllocate(0,0,255);   # make the background transparent and interlaced $im->transparent($white); $im->interlaced('true');   # Put a black frame around the picture $im->rectangle(0,0,99,99,$black);   # Draw a blue oval $im->arc(50,50...

问题 老师我想问一下,MENA在analysis the network 报错的问题

... Carter T. Butts, University of California-IrvineFor citation information, type citation("sna").Type help(package="sna") to get started.Error in `rownames<-`(x, value) :length of 'dimnames' [1] not equal to array extentCalls: row.names<- -> row.names<-.default -> rownames<-Executio...

问题 顺势作用原件功能评估

...ss.com/article/1154  Fig4中描述的比如顺势原件 MYB binding site involved in flavonoid biosynthetic genes regulation; LTR, cis-acting element involved in low-temperature responsiveness。请问老师如何评估这些元件的功能,文中也没有引用参考。我觉得评估这些功...

文章 组装好的基因组通过agp文件得到contig序列

...ibao/jcvi/wiki/ALLMAPS%3A-How-to-lift-over-gene-annotations But, you want the reverse mapping here so we'd need to swap the AGP file. python -m jcvi.formats.agp swap assembly.agp This generates a new AGP assembly.swapped.agp. You'll see that it basically swaps the components/objects from the or...

问题 docker toolbox启动报错

...) Check network to re-create if needed... (default) Windows might ask for the permission to configure a dhcp server. Sometimes, such confirmation window is minimized in the taskbar. (default) Waiting for an IP... Machine "default" was started. Waiting for SSH to be available... Detecting the pr...

问题 利用KaKs_Calculator2.0计算两个基因的kaks时出现:Error.

The size of two sequences in 'EVM0055308&EVM0007023' is not equal. 出现问题后也按照课堂上和视频课程中的方法进行了解决,也检查了ATG,序列是否为3的倍数等问题,均符合要求,不知道是什么问题造成的,

文章 如何用R语言绘制进化树+热图

... required.", 'height' , 'H', 1, "integer","the height of graph ,unit is inches [optional, default: 8]", 'width' , 'W', 1, "integer","the width of graph ,unit is inches [optional, default: 6]", 'outdir', '...

问题 群体结构分析时,执行代码报错

...vironment value of R_HOME --> Q&A for bioinformatics, please visit the website: https://www.omicsclass.com/ --> R beginners ? I suggest your  learning  R language: https://study.omicsclass.com/index pophelper v2.3.1 ready. Error in plotQ(qlist, imgoutput = "join", returnplot = T, exp...

文章 MutMap+和mutmap分析差异(BSA)

...tructed by replacing nucleotides in Nipponbare with those of Hitomebore at the 124,968 SNP positions identified between the two cultivars by alignment of 12.25 Gb of Illumina Hitomebore short reads to the Nipponbare reference genome要使用mutmap的定位方法:1.一般为隐性突变,最好是...

问题 gatk call SNP样品多采取方法2,报错如下是怎么回事啊

A USER ERROR has occurred: Argument emit-ref-confidence has a bad value: Can only be used in single sample mode currently. Use the --sample-name argument to run on a single sample out of a multi-sample BAM file