找到约 15 条结果

文章 叶绿体基因做跨物种系统发育分析

... Zhang Z, He Z, Xu S, et al. Transcriptome analyses provide insights into the phylogeny and adaptive evolution of the mangrove fern genus Acrostichum[J]. Scientific reports, 2016, 6: 35634.

文章 linux命令之AWK

...上数字/[0-9][0-9]*/    一个或一个以上数字FILENAME 文件名OFS      输出字段分隔符， 默认也是空格，可以改为制表符等ORS        输出的记录分隔符，默认为换行符,即处理结果也是一行一行输出到屏幕-F'[:#/]'   定义三个分隔符...

问题 我共线比对结果就这么点是为什么

...LUE: 1e-05# MAX GAPS: 25############### Statistics ################ Number of collinear genes: 12, Percentage: 0.05# Number of all genes: 23995############################################ Alignment 0: score=277.0 e_value=0 N=6 NW_014685782.1&NW_014685782.1 minus  0-  0: XP_014782229.1 XP_014...

文章 R语言基础入门—数据结构（2）

...R 语言的数据框可以使用 data.frame() 函数来创建。 # Create the data frame. emp.data <- data.frame( emp_id = c (1:5), emp_name = c("Rick","Dan","Michelle","Ryan","Gary"), salary = c(623.3,515.2,611.0,729.0,843.25), start_date = as.Date(c("2012-01-01", "2013-09-23", "20...

文章 R及R包的更新、移除、卸载

...例： detach(""package:R包名"", unload=TRUE)#detach: remove it from the search() path of available R objects; unload=T: to attempt to unload the namespace when a package is being detached. 如果要移除所有已加载的包，可参考https://stackoverflow.com/questions/7505547/detach-all-pa...

问题 asv分类报错

... set. Current permissions are '0o41755.' This most likely means something other than QIIME 2 created the directory '/home/wy/rawdata/am-test/data/my_ampliseq/tmp/qiime2' or QIIME 2 failed between creating '/home/wy/rawdata/am-test/data/my_ampliseq/tmp/qiime2' and setting permissions on it.The above ...

文章 linux 硬盘分区 格式化 挂载

...inux 2.23.2).Changes will remain in memory only, until you decide to write them.Be careful before using the write command.Device does not contain a recognized partition tableBuilding a new DOS disklabel with disk identifier 0xeb311dde.Command (m for help): m #打印帮助Command action   a  ...

文章 SCTransform ,Error: useNames = NA is defunct

...ixed slope and excluding poisson genes.Variance stabilizing transformation of count matrix of size 20445 by 2713Model formula is y ~ log_umiGet Negative Binomial regression parameters per geneUsing 2000 genes, 2713 cellsError: useNames = NA is defunct. Instead, specify either useNames = TRUE or useN...

问题 重测序

...bam \ failed with ret: 32512, going to retry. sh: -c: line 0: unexpected EOF while looking for matching `"' sh: -c: line 1: syntax error: unexpected end of file sh: -O: command not found warning, cmd: sh: -c: line 1: syntax error: unexpected end of file warning, cmd:   echo "RUN CMD: gatk --j...

文章 GEO、TCGA多数据库联合挖掘胰腺导管腺癌预后关键基因

...ancreatic carcinoma identified by co-expression analysis[J]. International Journal of Biological Sciences, 2018, 14(2):124-136. 相关课程： TCGA-基因差异表达分析、 TCGA-生存分析 、 GEO芯片数据挖掘、GEO芯片数据标准化、WGCNA-加权基因共表达网...

文章 宏病毒组分析一区

阿尔兹海默症病毒组——A metagenomic study of gut viral markers in amyloid-positive Alzheimer's disease patients 推荐语：越来越多的证据表明，病毒参与了阿尔茨海默症（AD）的发展和治疗。然而，对AD患者肠道病毒组的研究仍然存在显著的研...

文章 awk筛选差异基因

...选条件：log2fc的绝对值大于1，FDR<0.05 命令： awk 'BEGIN{OFS=FS="\t"}{if(FNR==1) print $0; else {abs_log2fc=($2<0?$2*(-1):$2);if(abs_log2fc>=1 && $3<0.05) print $0;}}' de_gene.txt $2代表的是第二列的log2fc,你的文件的log2fc在哪一列就用$第几列...

文章 samtools使用

...ls都有哪些命令 $samtools Program: samtools (Tools for alignments in the SAM format) Version: 1.9 (using htslib 1.9) Usage: samtools <command> [options] Commands: -- Indexing dict create a sequence dictionary file faidx index/extract FASTA fqidx ...

问题 hmmer搜索结构域

...提取 #结构域位置 grep -v '#' WRKY_hmmerOut_final.txt|awk 'BEGIN{OFS="\t"}$10==1 {print $1,$18,$19 }' >domain1_final.bed grep -v '#' WRKY_hmmerOut_final.txt|awk 'BEGIN{OFS="\t"}$10==2 {print $1,$18,$19 }' >domain2_final.bed #截取序列 seqtk subseq  ../01.data_prepare/Mlu.gene.pe...

文章 总结：生化检测实验样本值不准的原因分析

...high light stress tolerance in Arabidopsis thaliana [J]. Journal of pineal research, 2018. 这个J Pineal Res期刊，2018年IF=11.613，也是很牛了！！ 学霸姐姐大致看了一下，这篇文章主要是讲拟南芥中褪黑激素的诱导及其在高光胁迫下作用，文献报道了...