找到约 15 条结果

文章 使用DNA序列做主成分分析(PCA)——R语言adegenet包

...的数据,则输入准确的分组信息 a <- read.table(file="pop.txt")                 #1列,对应每条序列的分组信息df.pca.scores$population <- factor(a$V1, levels=c("group1", "group2","group3",                                                ...

文章 KEGG GO 富集显著性图 名字太长,柱状图上文字解决

... "KEGG_PENTOSE_AND_GLUCURONATE_INTERCONVERSIONS",     "KEGG_METABOLISM_OF_XENOBIOTICS_BY_CYTOCHROME_P450",    "KEGG_PORPHYRIN_AND_CHLOROPHYLL_METABOLISM",    "KEGG_NEUROACTIVE_LIGAND_RECEPTOR_INTERACTION"  ),  pvalue = c(0.000396688, 0.000246396, 0.000508482, 0.005582234, 0.001678246))# 数...

问题 你好,打扰了,想请教一下我是mac版R studio,安装WGCNA总报错,不太看得懂哪里出了问题,麻烦指导一下,急,谢谢!

...of C/C++/Fortran:   ‘WGCNA’ Do you want to attempt to install these from sources? (Yes/no/cancel) Yes installing the source package ‘WGCNA’ 试开URL’https://cran.rstudio.com/src/contrib/WGCNA_1.69.tar.gz' Content type 'application/x-gzip' length 1168289 bytes (1.1 MB) ===========...

文章 单细胞转录组数据挖掘流程记录-CRC直肠癌(GSE178341)

...ftp.ncbi.nlm.nih.gov/geo/series/GSE178nnn/GSE178341/suppl/GSE178341%5Fcrc10x%5Ffull%5Fc295v4%5Fsubmit.h5" -O GSE178341_crc10x_full_c295v4_submit.h5 wget -c "https://ftp.ncbi.nlm.nih.gov/geo/series/GSE178nnn/GSE178341/suppl/GSE178341%5Fcrc10x%5Ffull%5Fc295v4%5Fsubmit%5Fcluster.csv.gz" -O GSE178341_c...

文章 10X的空间转录组数据报错

Error in .local(obj, ...) : cannot derive coordinates from non-numeric matrix #9240 coordinates <- Seurat::Read10X_Coordinates("tissue_positions_list.csv", filter.matrix = TRUE) is.numeric(coordinates) head(coordinates) tissue_positions_list.csv 这个文件不应该有表头的; ...

文章 核苷酸替换模型-GTR是所有模型中考虑参数最多最复杂的模型,其他所有模型是GTR模型的特例

...ations of parameters for DNA site substitution. The models are listed here from the least complex to the most parameter rich. Jukes-Cantor (JC, nst=1): equal base frequencies, all substitutions equally likely (PAUP* rate classification: aaaaaa, PAML: aaaaaa) (Jukes and Cantor 1969) Felsenstein 198...

文章 geo_gene_exp_download.r GEO数据自动下载与整理

...lette, --palette palette                        A palette name from RColorbrewer [default Accent]  -G gpl, --gpl gpl     GPL file for annotation [default None]  -s Gene Symbol, --symbol Gene Symbol                        Gene Symbol column name [default Gene Symbol]  -x...

文章 高分文章1

...背景和功能 文献:The mutational constraint spectrum quantified from variation in 141,456 humans 发表时间:2020.05 简介:对来自141,456个体中得到的125,748个全外显子组序列数据和15,708个全基因组序列数据中预测的功能缺失变异进行整合,揭示...

问题 误: 无法分配大小为7.4 Mb的矢量 Killed after 9 m

... 载入程辑包:‘BiocGenerics’ The following objects are masked from ‘package:parallel’:     clusterApply, clusterApplyLB, clusterCall, clusterEvalQ, clusterExport,     clusterMap, parApply, parCapply, parLapply, parLapplyLB, parRapply,     parSapply, parSapplyLB The fol...

文章 单细胞转录组数据挖掘流程记录-结直肠癌 CRC(GSE178341)

...ftp.ncbi.nlm.nih.gov/geo/series/GSE178nnn/GSE178341/suppl/GSE178341%5Fcrc10x%5Ffull%5Fc295v4%5Fsubmit.h5" -O GSE178341_crc10x_full_c295v4_submit.h5 wget -c "https://ftp.ncbi.nlm.nih.gov/geo/series/GSE178nnn/GSE178341/suppl/GSE178341%5Fcrc10x%5Ffull%5Fc295v4%5Fsubmit%5Fcluster.csv.gz" -O GSE178341_c...

文章 绘制venn图,2-7维

...er R version 4.2.3genes <- paste("gene",1:1000,sep="")set.seed(20210419)x <- list(A=sample(genes,300),          B=sample(genes,525),          C=sample(genes,440),          D=sample(genes,440),          E=sample(genes,440),          F=sample(genes,150),          G=...

问题 UPARSE to pick OTUs in QIIME

...o they use the same algorithm as UCLUST in otu picking, exept the function of low abundance clusters filter and chimera detection in "usearch"  and "usearch61"?  And what is more, what is the difference between  "usearch"  and "usearch61", just different version? Is there any methods using the a...

文章 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.  图中圆圈大小代表该物...

文章 微生物多样性分类柱状图按照分组进行排序脚本升级;tax_bar_plot.r

...ment( "-i", "--input", type="character",required=T,help="input the abuance of tax in each sample,[required]",metavar="filepath")parser$add_argument( "-a", "--xlab", type="character",required=F, default="sample",help="input xlab [default sample]",metavar="xlab")parser$add_argument( "-t", "--title", t...

文章 awk筛选差异基因

...选条件:log2fc的绝对值大于1,FDR<0.05 命令: awk 'BEGIN{OFS=FS="\t"}{if(FNR==1) print $0; else {abs_log2fc=($2<0?$2*(-1):$2);if(abs_log2fc>=1 && $3<0.05) print $0;}}' de_gene.txt $2代表的是第二列的log2fc,你的文件的log2fc在哪一列就用$第几列...