找到约 15 条结果

文章 空间转录组数据读取

...者从网络上下到的原始数据,fq文件,那我们如何使用r来分析这些数据呢?很多同学做单细胞,对于矩阵的读取没啥问题,但是空间转录组多了一个图片,那如何处理呢 fq数据 fq数据直接用肯定行,需要借助软件spaceranger进行预处...

问题 镜像安装,镜像安装

老师,您好!我尝试跟着教程执行docker pull omicsclass/pan-gene-family:v1.0命令行时根本下了镜像,请问一下您可以发一份本地版本镜像文件吗

文章 多分组venn图绘制脚本

绘制多分组venn图时发现好多包都支持,只能用原始的包,脚本如下 #!/share/work/biosoft/R/latest/bin/Rscriptargs <- commandArgs(trailingOnly = TRUE)library(venn)library(VennDiagram)library(purrr)# 读取命令行参数input_file <- args[1]output_prefix <- args[2]#...

文章 OTUtable标准化方法

...ormalize_table.htmlmkdir normalized_otu_css/normalize_table.py -i $workdir/5.pick_otu/pick_de_novo_otus/otu_table_clean.biom -a CSS -o normalized_otu_css/CSS_normalized_otu_table.biombiom convert -i  normalized_otu_css/CSS_normalized_otu_table.biom -o  normalized_otu_css/CSS_normalized_otu_table1....

文章 打破部门壁垒,提升跨部门协作效率,让沟通更加有效

...些工作在部门内沟通完全没问题,但是跨部门沟通就问题断;一项工作分配好由一个部门主导另一个部门协作,但因配合当,迟迟无法交付;跨部门工作经常出现多个部门重复做了同一份工作;工作出现问题,各部门相互推...

文章 进化树构建的方法及原理详解

...。有根树反映了树上基因或蛋白质进化的时间顺序,通过分析有根树的树枝的长度,可以了解同的基因或蛋白质以什么方式和速率进化。而无根树只反映分类单元之间的距离,而涉及谁是谁的祖先问题。做有根树需要指定外类...

问题 Launch_PASA_pipeline.pl脚本报错连上mysql?

...ing to MySQL db: Specie_name_pasa -*** Running PASA pipeine: * [Wed May 15 16:19:38 2024] Running CMD: /share/work/biosoft/PASApipeline/PASApipeline-v2.5.2/scripts/create_mysql_cdnaassembly_db.dbi -c alignAssembly.config -S '/share/work/biosoft/PASApipeline/PASApipeline-v2.5.2/schema/cdna_alignmen...

文章 O2PLS-DA分析

两组学数据的正交偏最小二乘法判别分析(2-way orthogonal Partial Least Square with Discriminant Analysis ) O2PLS 模型用于两个数据组间的整合分析,包括系统生物学组学间关联、分子调控机制-表型间关联等各种大数据组的内在联系都可通...

文章 单细胞转录组数据挖掘流程记录-非小細胞肺癌(NSCLC)(GSE127465)

...促进这些群体的转化研究,我们还对小鼠中的 TIM 进行了分析。在比较同物种的 TIM 时,我们发现树突状细胞和单核细胞之间的群体结构几乎完全一致;中性粒细胞亚群保守;巨噬细胞之间存在物种差异。相比之下,患者血液...

文章 GATK SPARK 环境安装 加速GATK

...run-example SparkPi # 测试Pi 计算 3.设置 spark java环境; 版本要太高,23 会报错 cp conf/spark-env.sh.template conf/spark-env.sh #添加下面的行export SCALA_HOME=/share/biosoft/scala/bin/export JAVA_HOME=/share/biosoft/java/jdk-17.0.9/export SPARK_HOME=/share/bioso...

文章 从原始测序数据开始分析单细胞数据

#umi_tools whitelist --stdin data/12-22-3_L2_1.fq.gz \ #                    --bc-pattern="(?P<discard_1>CTCTTTCCCT)(?P<cell_1>.{10})(?P<discard_2>ACGACGCTCTTCGAGTGATTGCTTGTGACGCCTT)(?P<cell_2>.{8})(?P<umi_1>.{12})T{4}.*" \ #                    --se...

文章 RagTag:同源锚定延伸

...affolding过程中,可能会根据同参数或图谱数据产生多个版本的基因组组装结果,该步骤可以将多个结果根据权重进行最终组装结果的生成。 ragtag.py merge <asm.fa> <scf1.agp> <scf2.agp> [...]  <asm.fasta>           assem...

问题 reseqv2.0镜像annovar报错

...现报错错误信息如下 [root@8fbb0720ada0  09:47:30 /work/my_reseq/5.var_ann]# table_annovar.pl $workdir/4.snp_indel/var_qc/pyb4.clean.vcf.gz $refdir      -buildver unknown -out $workdir/5.var_ann/snp      -remove -protocol refGene -operation g -nastring . -vcfinput NOTICE: the --polis...

文章 deseq_analysis.r 差异基因分析DESeq2

...k 输入基因表达矩阵文件,为fpkm 或者 tpm文件 用于相关性分析: ID TCGA-B7-A5TK-01A-12R-A36D-31 TCGA-BR-7959-01A-11R-2343-13 TCGA-IN-8462-01A-11R-2343-13 TCGA-BR-A4CR-01A-11R-A24K-31 TCGA-CG-4443-01A-01R-1157-13 TSPAN6 59.65411 32.83064 40.7596 39.84131 53....

文章 零假设与备择假设

...检查人有病’作为原假设就恰当了。 假设检验是统计分析的一种重要方法,正确的理解假设检验过程中由于更换原假设与备择假设所引发的结论的差异,对于我们理解假设检验的思想以及掌握其方法都是十分重要的。只有更...