找到约 15 条结果

问题 重测序

...################################## # 每个样品分开运行,生成gvcf for i in $(cat $workdir/data/data.txt); do   echo "RUN CMD: gatk --java-options '-Xmx100g' HaplotypeCaller -R $REF   \     -I $workdir/3.map/result/${i}.sorted.dedup.bam \     -O ${i}.g.vcf.gz --max-alternate-alle...

文章 meme安装报错

...libxml/tree.h:1124:3: note: expected ‘xmlBufferPtr’ but argument is of type ‘xmlBufPtr’   xmlBufferWriteQuotedString(xmlBufferPtr buf,   ^xsltutils.c: In function ‘xsltSaveResultToString’:xsltutils.c:1758:26: error: dereferencing pointer to incomplete type  *do...

文章 split_sample.r根据生存数据对样本随机分组

...potion] [-o outdir] [-n name] The complete data set is randomly divided into training set and test set optional arguments: -h, --help show this help message and exit -i input, --input input input matrix data [required] -e event, --even...

文章 Argument list too long报错

...文件,这可能效率较低,但可以确保不超过参数限制。 for file in /要链接的文件路径/*; do    ln -s "$file" /目标目录路径/done

文章 遗传图绘制:

...t;- seq(0, maxpos)## put labels on ruler at every 10 cMaxlab <- vector()for (lab in 0:maxpos) {  if (!lab %% 10) {    axlab <- c(axlab, lab)  }  else {    axlab <- c(axlab, NA)  }}outfile="genetic_map.pdf"lmv.linkage.plot(carrot,outfile,denmap=TRUE, cex.axis = 1, at.axis = at.axis, ...

文章 Evolview:进化树设置标签的颜色

...((mouse,rat),(chimp,human)));”。点击Annotation upload → upload data for coloring leaf labels,填写相关数据后提交。 命令格式说明 设置语句包括3列:location,color,commands,由制表符分隔。第三列可有可无。l  第一列:可以是一个叶节点...

文章 cutree对pheatmap返回结果实现聚类cluster划分

...merge,height, and labels, of appropriate content each.kan integer scalar or vector with the desired number of groupshnumeric scalar or vector with heights where the tree should be cut. 所以将热图结果中的list$tree_row作为tree输入即可: list=pheatmap(mat,scale = "row")row_cluste...

文章 SNP分析中由于基因组序列过长导致报错解决方法

...d "main" htsjdk.samtools.SAMException: Exception when processing alignment for BAM index A00253:355:H75GLDSX2:2:1334:20989:8625 1/2 150b aligned read.        at htsjdk.samtools.BAMFileWriter.writeAlignment(BAMFileWriter.java:124)        at htsjdk.samtools.SAMFileWriterImpl.addAlignment(SAMFi...

问题 合并提取后的domain序列之后,linux系统中的clustalw不能读出蛋白信息

...个,依次类推 perl script/domain_xulie.pl WRKY_hmm_out.txt Arabidopsis_thaliana.TAIR10.pep.all.fa WRKY_domain.fa 1.2e-28 把WRKY_domain.fa和WRKY_domain2.fa的序列copy到同一个fasta文件中(用Editplus软件,复制两个文件的序列粘贴到一个新建文件中),clustalw不...

文章 R语言中读入数据时指定每一列数据类型

...类型。如下所示: data <- read.table(filename, head=FALSE, as.is=TRUE, quote="", comment.char="#", sep="\t",colClasses=c("character",rep("numeric",10))) 这里指定了第一列为字符型,其余列都为数值型。这样,第一列再读入时依然是 "01"了。 此外,我们...

文章 ggplot中用字符窜传递变量名称列名等

... group=column ) ) } Another option (ggplot2 > 3.0.0) is to use the tidy evaluation pronoun .data to slice the chosen variable/column from the rates.by.groups data frame. library(ggplot2) theme_set(theme_classic(base_size = 14)) # created by @Moody_Mudskipper rates.by.grou...

文章 吐血推荐!来自科研老司机的国庆7天出游指南!

...们, 作为过来人劝你们几句,还是老老实实地 家里游 or 学校游吧。。。 最后, 如果大家一定要出去玩儿的话, 那么到底去哪儿人少又便宜?小众而又彰显学术实力?还能震慑整个朋友圈? 只有一个答案: 哈哈,开...

文章 转录组差异基因为什么要进行聚类分析?

...   聚类使用的为R中的聚类软件包,所针对的数据为union_for_cluster(差异基因的并集),以基因的相对表达水平值log2(ratios) 进行聚类。其采用相应的距离算法,算出每个基因之间的距离,然后通过反复迭代,计算基因之间的相对...

问题 基因家族分析,提取一个物种信息后,再次取提取另一物种cds 用基因的ID,作为序列的ID,总是无法运行,显示脚本不存在,但是LL查看脚本是存在的。

Can't open per1 script ''/work/gene_family/scripts//get_gene_longest_fa.p1'':No such file or directory

文章 uniref50 数据库下载及格式化

...l split_taxid_uniref.pl division.dmp nodes.dmp uniref90.fasta.gz ./ for i in PLN MAM INV PHG PRI ROD SYN UNA VRL VRT ENV BCT ;do diamond makedb --in ${i}_uniref90.fa -d  ${i}_uniref90.fa done