找到约 15 条结果

文章 PICRUSt2扩增子功能预测与差异比较分析

...ata, metadata = metadata, "type") 功能丰度差异比较分析 # Perform pathway DAA using ALDEx2 method# Please change column_to_rownames() to the feature column if you are not using example dataset# Please change group to "your_group_column" if you are not using example datasetdaa_resu...

文章 NCBI、ENA高通量测序数据下载-fastq数据

....gov/ ;NCBI是由美国国家生物技术信息中心(National Center for Biotechnology Information)开发并负责维护,隶属于美国国立卫生研究院(National Institutes of Health, NIH)。 ENA:https://www.ebi.ac.uk/ena/browser/home;欧洲核苷酸序列数据库(European...

文章 施一公:如何提高英文的科研写作能力

...纸, “This is the draft. Please let me know what you think. We can aim for a Science report.”(这是初稿,你看看如何,我们可以试试《科学》)我仔细一看,天啊!一共7页,四个多小时Jeremy已经把文章的整体写完了,只是缺少Method和 references。...

问题 GWAS分析排序后出现sortPositions: false keepDepth: false的报错并且没有结果输出是什么原因

...tassel/latest/lib/fastutil-8.2.2.jar:/share/work/biosoft/tassel/latest/lib/forester-1.039.jar:/share/work/biosoft/tassel/latest/lib/gs-core-1.3.jar:/share/work/biosoft/tassel/latest/lib/gs-ui-1.3.jar:/share/work/biosoft/tassel/latest/lib/guava-22.0.jar:/share/work/biosoft/tassel/latest/lib/htsjdk-2....

文章 awk正则表达式使用外部变量

...hell变量,先在变量外加双引号,再加单引号即可,如: for i in read write print; do cat file.txt | awk '{ if ($7~/'"$i"'/) {print $0} }'| wc -l ; done 在写shell脚本时,经常会使用到awk程序。但是有些复杂的逻辑,可能需要在awk中使用在shell中...

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

...lization(cor.matrix,stresscent)centralization(cor.matrix,evcent)sna: Tools for Social Network AnalysisVersion 2.3-2 created on 2014-01-13.copyright (c) 2005, Carter T. Butts, University of California-IrvineFor citation information, type citation("sna").Type help(package="sna") to get started.Error i...

文章 Logistic回归模型介绍及使用(附脚本)

...iate ...], --variate variate [variate ...] variate for cox analysis[required] -l pvalue, --pvalue pvalue pvalue cutoff to choose sig gene[default 0.05] --log2 whether do log2 transfrom for expression data[opti...

问题 gVCF文件转换成VCF格式时,已经在后台运行了很久了,突然程序关闭了,报错信息显示:

...llegalStateException: No likelihood sum exceeded zero -- method was called for variant data with no variant information.         at org.broadinstitute.hellbender.tools.walkers.genotyper.AlleleSubsettingUtils.calculateMostLikelyAlleles(AlleleSubsettingUtils.java:31>         at org.broadi...

问题 重复序列注释问题

...组重复序列注释时,出现如下错误该如何解决?  W tensorflow/stream_executor/platform/default/dso_loader.cc:64] Could not load dynamic library 'libcuda.so.1'; dlerror: libcuda.so.1: cannot open shared object file: No such file or directory; LD_LIBRARY_PATH: /share/work/biosoft/mpi...

问题 docker运行qiime1时按照教程进行到章节2课时6时,即聚类生成OTU-denovo时老是报错,查看类似问题后,升级了运行内存为70g,同时将教程内的30多个样本减成2个样本,FASTQ格式还有qimme格式检查都正确,依旧报错,报错内容如下:

...n2.7/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.   warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.') Traceback (most recent call last):   File "/biosoft/minico...

文章 bioperl批量计算蛋白质分子量等电点pI

...蛋白序列就可以:蛋白分子量分析:ExPASy (http://web.expasy.org/protparam/)来计算: 但是这个网站,只能一条一条的提交序列然后计算,如果后很多条序列要提交,其不支持批量提交,所有功能有限,这里我写了个perl脚本,利用b...

文章 muscle软件介绍

...log: MUSCLE v3.8.31 by Robert C. Edgarhttp://www.drive5.com/muscleThis software is donated to the public domain.Please cite: Edgar, R.C. Nucleic Acids Res 32(5), 1792-97.1 2 seqs, max length 19, avg  length 1500:00:00    11 MB(-1%)  Iter   1  100.00%  K-mer dist pass 100:00:00    11...

文章 一行代码检查GWAS分析中表型数据的错误

... 12 11 7 运行该代码 awk 'NR == 1 {print "Line", NR, $0; next} {for (i = 2; i <= NF; i++) if ($i !~ /^[0-9.]+$/ && $i != "NA") {print "Line", NR, $0; break}}' trait.txt 得到结果如下,会展示出表头及除了数字和NA之外的所有字符以及所在的行数,...

文章 Evolview:设置标签的装饰

...05)100:0.1)90:0.43;”。 [2].   点击Annotation upload → upload data for leaf label decorations,在data文本框输入“A triangle,red”并提交。 [3].   结果展示如下 命令格式说明: 设置语句至少有2列:location shape,color[:stroke_color],由制表符分...

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

...singles> file cannot also be specified.  Please run HISAT2 separately for mates and singles. Error: Encountered internal HISAT2 exception (#1) Command: /home/vip6/miniconda3/envs/rna/bin/hisat2-align-s --wrapper basic-0 -p 3 -S /home/vip6/project/human-64-Asthma-Trans/Mapping/Hisat//SRR103950...