找到约 15 条结果

文章 多分组venn图绘制脚本

...pdf设备dev.off()# 导出结果inter <- get.venn.partitions(venn_list)for (i in 1:nrow(inter)) inter[i, "values"] <- paste(inter[[i, "..values.."]], collapse = "|")inter <- subset(inter, select = -..values..)inter <- subset(inter, select = -..set..)write.table(inter, paste0(output_prefi...

文章 TCGA癌症数据库挖掘发表的大文章CNS级别汇总2019

... Analysis Network continues to study and mine the data. Cancers Selected for Study lists original marker publications by cancer type. Supplemental and associated data files are located in the GDC.2008Comprehensive genomic characterization defines human glioblastoma genes and core pathwaysExit Di...

文章 cBioPortal数据库使用教程

...于该数据库的文章: The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data.Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. 大规模的癌症基因组计划,比如The cancer genome atlas(TCGA) and ...

文章 qiime2 分类器建立 SILVA数据库

...可以直接下载qiime2官方网站i提供的文件:https://docs.qiime2.org/2020.8/data-resources/ wget -c https://data.qiime2.org/2020.8/common/silva-138-99-seqs.qzawget -c https://data.qiime2.org/2020.8/common/silva-138-99-tax.qzaln -s silva-138-99-tax.qza silva-138-ssu-nr99-tax.qzaln -s silva...

文章 外显子比对率太低,参考基因组索引建立refFlat.txt文件有问题:

... -genePredExt saccharomyces_cerevisiae.gtf refFlat.tmp.txtno exons defined for group HRA1, feature noncoding_exon (perhaps try -ignoreGroupsWithoutExons) 改进: gtfToGenePred -genePredExt -ignoreGroupsWithoutExons saccharomyces_cerevisiae.gtf refFlat.tmp.txt

文章 Monocle2|单细胞测序的拟时序分析

...LL] --downsample downsample subset cells numbers for analysis [default None] -g group, --group group the group that coloring the cells by [default seurat_clusters] -d density.group, --density.group density.group ...

问题 (ERR): "--dta" does not exist Exiting now ...为什么hisat2有这个选项却报错

...载的是clean_data所以跳过质控直接做map,修改了代码如下 for i in SRR9113321 SRR9113322 SRR9113323 SRR9113324 SRR9113325 SRR9113326 SRR9113327 SRR9113328 SRR9113329 SRR9113330 SRR9113331 SRR9113332 SRR9113333; do echo "RUN CMD: hisat2 -p 1 --rg-id=${i} --rg SM:${i} --rg LB:${i} --...

问题 老师,您好!全基因组选择,利用自己的SNP数据(芯片SNP)进行基因型数据处理,beagle.27Feb25.75f.jar处理完成后,得到的imputed_snp里面没有SNP标记信息。

....jar   gt=GENO.vcf   out=imputed_snp   nthreads=48 No genetic map is specified: using 1 cM = 1 Mb java.lang.IllegalArgumentException: Duplicate marker: 1 77806149        AX-90568401     C       T         at vcf.Markers.markerSet(Markers.java:127)         at vcf.Marke...

问题 使用xpclr_smooth_line_plot.r绘制时,没有报错,但没有图片生成

...chr.xpclr -F $FAI -f 19226500  -n all.chr.xpclr.smoothline --> 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 Loading required package: ggplot2 chr_name =  1 chr_na...

文章 R 包 osqp安装报错

... Build files cannot be regenerated correctly. ERROR: configuration failed for package ‘osqp’ 解决版本,降低版本: wget https://cran.r-project.org/src/contrib/Archive/osqp/osqp_0.6.0.7.tar.gzR CMD INSTALL osqp_0.6.0.7.tar.gz

文章 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。...

文章 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分析排序后出现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....