找到约 15 条结果

文章 BWA比对及Samtools提取mapped reads

...: 1:Only READ1 is set. 2:Only READ2 is set. O:Either both READ1 and READ2 are set;or neither is set. 对于PE测序reads,同时指定-1 R1.fq-2 R2.fq-s singleton.fq 时,samtools会将flag1和flag2配对的序列分别输出到-1-2指定的文件,对于无法匹配的flag 1/2,...

文章 cox模型评估-决策曲线分析(DCA)

...shed 2019 Oct 4.3. Van Calster B, Wynants L, Verbeek JFM, et al. Reporting and Interpreting Decision Curve Analysis: A Guide for Investigators.Eur Urol. 2018;74(6):796-804.4. Yan P, Li JW, Mo LG, Huang QR. A nomogram combining inflammatory markers and clinical factors predicts survival in patients w...

文章 ConsensusClusterPlus 包对基因表达数据进行一致性聚类

...t(1 - fm), method = finalLinkage) # # Cluster method : average # Number of objects: 128 results[[2]][["consensusClass"]][1:5] # 01005 01010 03002 04006 04007 # 1 1 1 1 1 计算聚类一致性 (cluster-consensus) 和样品一致性 (item-consensus) icl <- ca...

文章 单细胞转录组数据挖掘流程记录-BRCA乳腺癌(GSE161529)

...10GenesPerUMI 0.7 \ -o 01.qc-$i -p GSE161529 --metadata.col.name Sample type --metadata.value $i ${i%%-*} done Rscript $scripts/merge_seurat_obj.r -i $(ls 01.qc-*/GSE161529.afterQC.rds) \ -o 02.merge -p GSE161529 Rscript $scripts/seurat_sc_cluster.r --rds 02.merge/GSE161529.rds \ -p GSE1615...

文章 picrust2用法

...ID中注释信息,避免错误:no ASV ids overlap between input FASTA and sequence abundance table get_fa_by_id.pl $workdir/5.pick_otu_qiime/pick_de_novo_otus/otu_table_clean.txt $workdir/5.pick_otu_qiime/pick_de_novo_otus/rep_set/qiime_rep_set.fasta rep_set.farm -rf picrust2_outpicrust2_pipeli...

文章 给一个fastq ID列表,输出指定ID对应的双端fastq文件:

.../);$keep{$tmp[0]}=1;}close(IN);my$i=0;while ( my $obj1=$fq1->next_seq() and my $obj2=$fq2->next_seq() ) {my$id1=$obj1->id;my$id2=$obj2->id;#print "$id1,$id2\n";#die;if( exists $keep{$id1} or exists $keep{$id2}){$fqo1->write_seq($obj1);$fqo2->write_seq($obj2);}}$fq1->close();$fq2...

文章 biolinux虚拟机不能为虚拟机打开一个新任务,或者一直黑屏解决办法

...拟电脑 Bio-Linux-8.0.7 打开一个新任务.VT-x is disabled in the BIOS for all CPU modes (VERR_VMX_MSR_ALL_VMX_DISABLED).返回 代码: E_FAIL (0x80004005)组件: ConsoleWrap界面: IConsole {872da645-4a9b-1727-bee2-5585105b9eed} 报错截图: 解决办法: 当...

文章 非肿瘤文章集锦

...RNA-mRNA调控网络 Maternal smoking-induced lung injuries in dams and offspring via inflammatory cytokines  https://pubmed.ncbi.nlm.nih.gov/33989842/  https://pubmed.ncbi.nlm.nih.gov/34306596/   自噬相关心血管疾病

文章 挖掘公共基因大数据的时代已经到来

...数据的挖掘典范,当属TCGA数据库。 TCGA一A座金矿 TCGA(The Cancer Genome Atlas )数据库,目前收录了33种癌症,11000多列病人的多个组学的数据以及临床表型信息。组学数据包括基因组的变异,转录组,miRNA, 甲基化,蛋白质组等。...

文章 dcaps标记设计软件——SNP_Primer_Pipeline

...CGGGAGAGGAAATTACTGCCTTGGTGCTG 输出结果: index   product_size    type    start   end     diff_number     3'differall     length  Tm      GCcontent       any 3'       end_stability   hairpin primer_seq      ReverseComplement       penalty compl_any       compl_e...

文章 R语言基础入门—运算符!

...结果是布尔值。 运算符描述例& 它被称为元素逻辑AND运算符。它将第一向量的每个元素与第二向量的相应元素组合,并且如果两个元素都为TRUE,则给出输出TRUE。 v <- c(3,1,TRUE,2+3i)t <- c(4,1,FALSE,2+3i)print(v&t)#它产生...

问题 qiime2 python "matplotlib" _tkinter.TclError: couldn't connect to display "localhost:12.0" 报错

...te-packages/qiime2/sdk/action.py", line 245, in bound_callable    output_types, provenance)  File "/share/work/biosoft/conda/Miniconda3.8/envs/qiime2-2020.8/lib/python3.6/site-packages/qiime2/sdk/action.py", line 452, in _callable_executor_    ret_val = self._callable(output_dir=temp_dir, **vie...

文章 ValueError: invalid literal for int() with base 10: ''

...如下: /biosoft/miniconda/miniconda2/bin/python -m jcvi.graphics.karyotype  --format=pdf  --figsize=15x5 mcscan_seqid mcscan_layout 报错信息如下: markers', 'networkx.algorithms.community.__future__', 'networkx.generators.atlas', 'networkx.algorithms.connectivity.edge_kcomponents', 'p...

文章 python特殊属性和方法的详解与运用

...'__main__' >>> Class1.__bases__ # 类型所继承的基类 (<type 'object'>,) >>> Class1.__dict__ # 类型字典,存储所有类型成员信息。 <dictproxy object at 0x00D3AD70> >>> Class1().__class__ # 类型 <class '__main__.Class1'> >>...

文章 rocks mysql数据库多实例数据库配置

...sql.sock -e "show variables like 'server_id';"Warning: Using a password on the command line interface can be insecure. +---------------+-------+ | Variable_name | Value | +---------------+-------+ | server_id     | 7     | +---------------+-------+ https://www.cnblogs.com/lon...