找到约 15 条结果

文章 NCBI的SRA Lite和SRA Normalized数据有什么区别

...ssing overall read quality, setting a per-read quality flag (Read_Filter), and removing base quality scores from the file. In the resulting files, all reads have a Read_Filter flag with value pass or reject. Importantly, it is still possible to produce fastq formatted files from SRA Lite format usin...

问题 edger差异分析,如何查询分组差异的结果

...LSE]##归一化TMM y <- calcNormFactors(y,method = 'TMM')#差异分析 #design <- model.matrix(~groups) ###计算离散度,估计离散值(Dispersion) ##估算方法1 #y <- estimateDisp(y,design)####一次性计算common dispersion和tagwise disperson ##估算方法2, y <- estimate...

文章 iTOL对进化树添加外部标签

...e node branch, or outside the tree#lines starting with a hash are comments and ignored during parsing#=================================================================##                    MANDATORY SETTINGS                           ##=========================================...

文章 MEGAN 分类学树状图 饼图说明

 Each circle of the MEGAN tree represents a taxon in the NCBI taxonomy and is labeled by its name and the number of reads that were assigned to the taxon and not to a subtaxon. The size of the circle is scaled logarithmically to represent the number of reads.  图中圆圈大小代表该物...

文章 R语言基础入门—数据结构(2)

...ingsAsFactors = FALSE)# Add the "dept" coulmn. emp.data$dept <- c("IT","Operations","IT","HR","Finance") v <- emp.dataprint(v) 当我们执行上面的代码,它产生以下结果 - emp_id emp_name salary start_date dept 1 1 Rick 623.30 2012-01-01 IT ...

文章 GEO数据芯片数据挖掘-在线工具分析GDS数据

...sis Tools 可以查看数据整体的分布情况,例如点击Experiment design and value distribution,可以出现数据box分布图,如下: 单击图片位置,可以查看整个GDS所有样品的数据分布情况,同时查看样品到处理信息(如下图)。 总结 针对...

文章 NCBI批量下载数据,省时又省力

...   os.mkdir(args.out_dir)    dout=os.path.abspath(args.out_dir)output_handle = open(dout+'/'+args.name+'.%s'%args.rettype, "w")Entrez.email = "huangls@biomics.com.cn"     # Always tell NCBI who you are#handle = Entrez.efetch(db="nucleotide", id="EU490707", rettype="gb", retmode="text")#print(h...

文章 lasso_cox.r lasso_cox 模型优化分析

...alysis optional arguments: -h, --help show this help message and exit -i data, --data data input data file path[required] -t time, --time time set suvival time column name [required] -e event, --event event set event column name [required] -v variate [...

文章 R安装units包报错解决方法

...ource* package ‘units’ ...** package ‘units’ successfully unpacked and MD5 sums checked** using staged installationconfigure: units: 0.8-2checking whether the C++ compiler works... yeschecking for C++ compiler default output file name... a.outchecking for suffix of executables...checking whe...

问题 RNAseq有参转录组数据自主分析课程,在建立基因组索引的时候出现如下问题

...务器上运行,最后获取gene length时候出错: get gene length and gene.bed from gtf: RUN CMD: python /home/zpp/work/my_rnaseq/scripts/get_gene_length_from_gtf.py -g Homo_sapiens.GRCh38.99.chromosome.22.gtf -p gene_length Traceback (most recent call last):   File "/home/zpp/work/my_rn...

文章 R 包RIdeogram安装报错

...ackage ‘grImport2’ ...** package ‘grImport2’ successfully unpacked and MD5 sums checked** using staged installation** R** inst** byte-compile and prepare package for lazy loading** help*** installing help indices** building package indices** testing if installed package can be loaded from te...

文章 芯片差异分析分组矩阵

...数完成,获得的分组矩阵结构类似下方显示结果: > design           A BGSM1626001 1 0GSM1626002 1 0GSM1626003 1 0GSM1626010 0 1GSM1626011 0 1GSM1626012 0 1attr(,"assign")[1] 1 1attr(,"contrasts")attr(,"contrasts")$`description`[1] "contr.treatment" ...

文章 保姆式教程!手把手教你使用Simplot,绘制高大上的序列相似性图片

...很美观!!! 参考文献: 《Monsavirus in monkey rectal swab and throat swab specimens in China: Proposal for Posaliviridae as a new family in Picornavirales》; Han Z, Song Y, Xiao J, Zhao X, Lu H, Zhang K, Jia S, Zhou J, Li J, Si F, Sun Q, Zhu S, Wang D, Yan D, Xu W, Fu X, Zhang Y. ...

文章 ensembl和NCBI基因组下载,基因序列下载查看

...ensemblgenomes.org/ 2.NCBI     是National Centerfor Biotechnology Information的缩写,指美国国立生物技术信息中心。NCBI的全面和强大,相信大家都深有感触,NCBI在参考基因组信息展示上同样表现出色。 地址:https://www.ncbi.nlm.nih.gov/ 详...

问题 WGCNA软阈值计算报错

老师,您好!我在WGCNA第二步,network construction中软阈值计算的过程中遇到了问题,暂时不知道如何解决,想请教下各位老师!!