找到约 15 条结果

文章 转录组中Count, TPM,FPKM如何计算

...量获得的read count 数,是如何计算成TPM和FPKM的呢?参考《What the FPKM? A review of RNA-Seq expression units》这篇博客, 1. TPM 公式如下: 2. FPKM 公式如下: 3. TPM与FPKM的关系 对应的R 代码如下: countToTpm <- function(counts, effLen) { ...

文章 R 包安装报错:configure: error: gdal-config not found or not executable.

... * removing ‘/share/work/biosoft/R/R-v3.6.1/lib64/R/library/adegenet’ The downloaded source packages are in ‘/tmp/RtmpuYCUHI/downloaded_packages’ Updating HTML index of packages in '.Library' Making 'packages.html' ... done Warning messages: 1: In install.packages("adegenet") : in...

问题 juicer 报错: Failure during chimera handling of

... 老师好,在运行juicer时报错: awk: calling undefined function and input record number 1116, file /home/data/t090519/JF/20.Jui/juicer/work/splits/JF2.fastq.gz.sam source line number 298***! Failure during chimera handling of /home/data/t090519/JF/20.Jui/juicer/work/splits/JF2.fastq.gz...

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

...(description='plot otu tax bar plot')parser$add_argument( "-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]",me...

问题 R包安装

...n of this package for your version of R might be available elsewhere, see the ideas at https://cran.r-project.org/doc/manuals/r-patched/R-admin.html#Installing-packages    请问老师,出现这种情况应该怎么处理?

问题 WGCNA-MAD

...t.char = "",check.names=F) datExpr = datExpr[,-1] # Take a quick look at what is in the data set dim(datExpr) ##################check missing value and filter ################### datExpr = t(datExpr[order(apply(datExpr,1,mad), decreasing = T)[1:5000],]) dim(datExpr) 根据上面的代码,...

文章 基因组Hic挂载AGP文件格式说明

...片段和gap在大片段上的编号。 5. 该段序列的类型(component_type)     常用的是W、N和U。W表示WGS contig;N表示指定大小的gap;U表示不明确长度的gap,一般用100bp长度。 6. 小片段的ID或gap长度(component_id or gap_length)     如果第5列不...

文章 分离NR NT库,快速blast本地比对同源注释基因

...有脚本内存消耗较大,建议确保100G以上的内存空间: #The gi_taxid_nucl.dmp is about 160 MB and contains two columns: the nucleotide's gi and taxid. #The gi_taxid_prot.dmp is about 17 MB and contains two columns: the protein's gi and taxid. #Divisions file (division.dmp): # ...

文章 mapman非模式植物基因注释时遇到的问题

...。 Next,:: about Mercator - PlabiPD在注释网站发现一个大坑,The Mapman desktop application requires a diagram that illustrates the pathway or process to be visualized. Each diagram consists of a pathway image file (.svg) and an extra protein/gene mapping file (.xml). Notice, the diagr...

文章 R 安装github包报错:

...set()` and unset GITHUB_PAT in .Renviron (or elsewhere) if you want to use the more secure git credential store instead. Error: Failed to install 'unknown package' from GitHub:   HTTP error 401.   Bad credentials   Rate limit remaining: 53/60   Rate limit reset at: 2024-11-18 08:36:57 UTC...

文章 GATK报错| this implies a problem with the file locking support in your operating system

... was reached while trying to acquire a lock on file raw.vcf.tmp.idx. Since the GATK uses non-blocking lock acquisition calls that are not supposed to wait, this implies a problem with the file locking support in your operating system. ##### ERROR -----------------------------------------------------...

文章 cnv_dumbbell_plot.r 绘制哑铃图ggplot2包

... -g GENE_SYMBOL, --gene_symbol GENE_SYMBOL Enter the file path that contains the gene symbol[required] -H HEIGHT, --height HEIGHT the height of dumbbell plot[optional,default 6] -W WIDTH, --width WIDTH ...

文章 bedtools批量提取基因组指定位置序列

...  BED/GFF/VCF file of ranges to extract from -fi        -name   Use the name field for the FASTA header        -split  given BED12 fmt., extract and concatenate the sequencesfrom the BED "blocks" (e.g., exons)        -tab    Write output in TAB delimited format.              ...

文章 enrichGSEA_pip.r GSEA富集分析

...fault:10 ]  -n prefix, --prefix prefix                        the output file prefix [optional, default: GSEA ]  -o outdir, --outdir outdir                        output file directory [default cwd]  -H height, --height height                        the height ...

问题 老师好,我在做群里遗传进化算pca时,得到的文件都是0。但是没有报错

...iants loaded from .bim file.2045265 missing IDs set.684 people (0 males, 0 females, 684 ambiguous) loaded from .fam.Ambiguous sex IDs written to plink_pca.nosex .Using up to 23 threads (change this with --threads).Before main variant filters, 684 founders and 0 nonfounders present.Calculating allele...