找到约 15 条结果

文章 pROC包绘制ROC曲线

...s100b对outcome进行分析 > outcome=roc(response=aSAH$outcome,predictor = aSAH$s100b)> outcomeCall:roc.default(response = aSAH$outcome, predictor = aSAH$s100b)Data: aSAH$s100b in 72 controls (aSAH$outcome Good) < 41 cases (aSAH$outcome Poor).Area under the curve: ...

文章 blast在linux上的用法

...名称 i：要查询(要研究的物种序列)的序列文件名(Query File)，fasta格式 -e(数学)：期望值(Expectation value)，E值是个统计阈值，缺省值10，意指比对结果中由于随机偶然性产生的匹配结果不大于10,E值越小结果越可靠 -o：查询结...

文章 统计不同样品-生物学重复基因表达量数据合并和绘图统计

...<-read.table("All_gene_fpkm.xls",header=TRUE,comment.char="",sep = "\t",check.names=FALSE,row.names=1)head(myfpkm)group=factor(c(rep(c("DAF2", "DAF5", "DAF11", "DAF16"),each=3),rep(c("GDAF2", "GDAF5", "GDAF11", "GDAF16"),each=3)),levels = c("DAF2","DAF5","DAF11","DAF16","GDAF2","GDAF5","GDAF11","...

问题 请问老师有根据基因ID一次性提取UTR,EXON,INTRON的脚本吗？

...========= IDs that you used as input do       NOT present in gff/gtf file.      Examples of valid IDs is showed as bellowedGuess FileType is:GTFtranscript:AT1G03987.1transcript:AT1G01020.6transcript:AT1G01010.1transcript:AT1G01020.3transcript:AT1G01020.2transcript:AT1G01020.5transcript:AT1...

问题 vcftools过滤除了等位基因选项外，其他参数通通没用

...fter filtering, kept 560 out of 560 Individuals      26 Outputting VCF file...      27 After filtering, kept 118111 out of a possible 118111 Sites      28 Run Time = 121.00 seconds 请问老师可能原因是什么？

文章 python 版本的MCScanX分析物种之间分共线性

...suming --qbed=grape.bed --sbed=peach.bed 20:34:31 [blastfilter] Load BLAST file `grape.peach.last` (total 403868 lines) 20:34:31 [base] Load file `grape.peach.last` 20:34:38 [blastfilter] running the cscore filter (cscore>=0.70) .. 20:34:39 [blastfilter] after filter (294217->31229) .. 20:34:3...

文章 单细胞转录组数据挖掘流程记录-肝细胞癌（HCC）(GSE149614)

...xt.gz \ --nUMI.min 500 \ --nUMI.max 75000 \ --nGene.min 250 \ --mito.gene.pattern "^MT.*-" \ --percent_mito 20 \ --log10GenesPerUMI 0.7\ -o 01.qc -p GSE149614 #单细胞聚类分群分析 Rscript $scripts/seurat_sc_cluster.r --rds 01.qc/GSE149614.afterQC.rds \ -p GSE149614 --reso...

文章 B/T细胞分类：

...张泽民课题组描绘泛癌B细胞全景图） Tian, C., Zhang, Y., Tong, Y. et al. Single-cell RNA sequencing of peripheral blood links cell-type-specific regulation of splicing to autoimmune and inflammatory diseases. Nat Genet 56, 2739–2752 (2024). https://doi.org/10.1038/s41588-024-0...

文章 IF=3.47|非肿瘤-眼部老年黄斑变性ceRNA文章思路

...-miRNA-mRNA network in age-related macular degeneration: From construction to identification. Exp Eye Res. 2021 Feb;203:108427. doi: 10.1016/j.exer.2020.108427. Epub 2020 Dec 29. PMID: 33383027. 延伸阅读 GEO数据库挖掘—WGCNA鉴定骨肉瘤转移相关基因GEO、TCGA多数据库联合挖...

文章 群体遗传进化与GWAS分析动植物文献

...    文章:Whole-genome sequencing of giant pandas provides insights into demographic history and local adaptation·        基础数据:34只熊猫,4.7x覆盖深度 NC-牦牛群体进化-2015·        文章:Yak whole-genome resequencing reveals domestication signatures and prehis...

文章 tar 多线程加速打包速度

...数 bash# 排除单个文件或目录 tar -czvf archive.tar.gz --exclude="file.txt" folder/ tar -czvf archive.tar.gz --exclude="node_modules" project/ # 排除多个文件或目录 tar -czvf archive.tar.gz --exclude="*.log" --exclude="temp" --exclude=".git" folder/ # 使用模式匹配 tar -czvf ...

文章 IF=4.8 | 不同菊花品种萜类合成相关基因的转录组代谢组联合分析

...2月 BMC Plant Biology杂志发表了一篇题为“Integrative transcriptome and metabolome analysis reveals candidate genes related to terpene synthesis in Chrysanthemum × morifolium”的转录组代谢组联合分析文章。” 01 — 研究背景 菊花（Chrysanthemum × morifolium）...

文章 VirtualBox虚拟机网络设置（四种方式）

...口，映射到主机的 80 端口。使用下面的命令。 'C:\Program Files\innotek VirtualBox\VBoxManage.exe' setextradata 'debian' 'VBoxInternal/Devices/pcnet/0/LUN#0/Config/huzhangsheng/Protocol' TCP 'C:\Program Files\innotek VirtualBox\VBoxManage.exe' setextradata 'debian' 'VBoxInternal/De...

文章 TCGA数据库中应该下载哪种表达量数据HTSeq-Counts，HTSeq-FPKM，HTSeq-FPKM-UQ

...alculation normalizes read count by dividing it by the gene length and the total number of reads mapped to protein-coding genes. Upper Quartile FPKM The upper quartile FPKM (FPKM-UQ) is a modified FPKM calculation in which the total protein-coding read count is replaced by the 75th percentile read...

文章 bioperl批量计算蛋白质分子量等电点pI

...:pICalculator;use Data::Dumper; #读入序列 my $in = Bio::SeqIO->new(-file   => "$ARGV[0]",-format => 'Fasta');open OUT,">$ARGV[1]" or die "$!";print OUT "#ID\tlength\tMV(Da)\tpI\n";my $calc = Bio::Tools::pICalculator->new(-places => 2,-pKset => 'EMBOSS');#逐条读取序列...