找到约 15 条结果

文章 blast序列比对中E值是什么意思?

...。E值接近零或为零时,具本上就是完全匹配了 E-value - this tells you the probability of getting this score and homology by chance (from random sequence), that is, the lower (the higher the number behind e-), the more reliable the result is. As I have mentioned above, it's not only...

文章 nomogram_multi_cox.r 多因素cox分析构建模型并绘制列线图nomogram

...                        set event column name must 0 or 1 code format                        [required]  -v variate [variate ...], --variate variate [variate ...]                        variate for cox analysis [required]  -P predict.time [predict.time ...], --p...

文章 外源插入流程报错运行不了

...: You may be importing Biopython from inside the source tree. This is bad practice and might lead to downstream issues. In particular, you might encounter ImportErrors due to missing compiled C extensions. We recommend that you try running your code from out...

文章 python版本MCScanX分析方法

...d C.baccatum The genomes of C. annuum, C.chinenese and C.baccatum were retrieved from PGP (http://peppergenome.snu.ac.kr/download.php).The C. annuum genome was compared to C.chinenese and C.baccatum  genomes using the MCScan toolkit (V1.1)[1]. To call synteny blocks, we performed all-against-al...

文章 常见的植物数据库(更新-Ing)

1、NCBI:https://www.ncbi.nlm.nih.gov/ 2、Uniprot:https://www.uniprot.org/ 3、KEGG:https://www.kegg.jp/ 4、InterPro:https://www.ebi.ac.uk/interpro/ 5、DDBJ:https://www.ddbj.nig.ac.jp/index-e.html 6、 Phytozome :https://phytozome.jgi.doe.gov/pz/portal.html# 7、http://plants.ensembl.o...

问题 NGS数据分析时,使用ANNOVAR的table_annovar.pl报错。

...ne with 1609 transcripts (including 43 without coding sequence annotation) for 1609 unique genesNOTICE: Processing next batch with 6421 unique variants in 6421 input linesNOTICE: Reading FASTA sequences from HPdb/HP_refGeneMrna.fa ... Done with 220 sequences 但是使用table_annovar.pl进行注释...

问题 基因家族分析circos时显示基因基因数目超出最大值怎么办

...       /home/manager/gene_family/circos/genome.txt   but the maximum is currently set at [25000]. To increase this number change   max_links in etc/housekeeping.conf. Keep in mind that drawing that many links   may create an image that is too busy and uninterpretable. 请问出现如...

问题 GWAS筛选阈值设置

For SNPs, we set the thresholds as −log10(P value) = 6, 6, and 4 for litter size, numbers of horns and nipples, respectively. 请问,这种自选的阈值,怎么选择呢?为什么有的是6有的是4呢?

文章 MFP - Molecular Functional Portrait

...ing_profile_metrics_plotfrom portraits.utils import read_gene_sets, ssgsea_formula, median_scale# Example script# Read signaturesgmt = read_gene_sets('signatures.gmt')  # GMT format like in MSIGdb# Read expressionsexp = pd.read_csv('expression.tsv', sep='\t', index_col=0)  # log2+1 transformed; Ge...

文章 topGO报错

...     pval.info <- function(whichNodes) {          ret.val <- format.pval(termsP.value[whichNodes], digits = 3, eps = 1e-20)          names(ret.val) <- whichNodes          return(ret.val)        }        library("DBI")        .getTermsDefinition <- function(wh...

问题 blastp提示ensemble中提取的序列文件不对

...ulate ungapped Karlin-Altschul parameters due to an invalid query sequence or its translation. Please verify the query sequence(s) and/or filtering options 已去除序列中的终止子符号。 e.g: >ENST00000620518 ENSG00000081692  JMJD4 1:227732392-227735273(-) MDRETRALADSHFRGLGVDVPGV...

问题 WGCNA-MAD

...############ input data########################## # The following setting is important, do not omit. options(stringsAsFactors = FALSE) #读取数据,删除第一列GeneSymbol datExpr = read.table(expr_file,header=T,row.names=1,comment.char = "",check.names=F) datExpr = datExpr[,-1] # Take a ...

文章 RagTag:同源锚定延伸

...   build a FASTA file from an AGP file      agpcheck        check for valid AGP file format      asmstats        assembly statistics      splitasm        split an assembly at gaps      delta2paf       delta to PAF file conversion      paf2delta       PAF to delta ...

文章 ANNOVAR琳琅满目的人类变异信息注释数据库

...位置(chromEnd)、名称(name)、染色标识(gieStain)‘1000g2015aug’ for alternative allele frequency in the 1000 Genomes Project 。1000基因组项目(和ExAV 外显子集合联合一样,是公开、开放的数据库)里面供选择的等位基因频率信息。‘exac03’for the var...

问题 R语言分析GEO数据,您好老师,19年参加培训的时候就是用的R3.6.1,自己运行出现如下报错

...0709cdf' Installation path not writeable, unable to update packages: boot, foreign, KernSmooth, lattice, Matrix, mgcv, nlme, nnet, survival Old packages: 'BH', 'bibtex', 'BiocParallel', 'bit', 'blob', 'broom', 'callr', 'caret', 'caTools', 'checkmate', 'chron', 'cli', 'clusterProfiler', 'data.table...