找到约 15 条结果

文章 blast报错

...复制粘贴官网的代码并执行: makeblastdb –in mydb.fsa –dbtype nucl –parse_seqids 总会有如下错误提示: USAGE  makeblastdb [-h] [-help] [-in input_file] [-input_type type]    -dbtype molecule_type [-title database_title] [-parse_seqids]    [-hash_index] [-mask_data m...

问题 您好!我是R语言的新手,在用hchinamap程序绘制地图时,出现cannot coerce type 'closure' to vector of type 'character',出现这类错误怎样解决呢?请指教,谢谢!

步骤如下: install.packages('hchinamap') download.file('https://mdniceczx.oss-cn-beijing.aliyuncs.com/chinadf.rda', file.path(dir, 'chinadf.rda')) load(file.path(dir, 'chinadf.rda'), verbose = TRUE)chinadf

问题 RNA 速率分析Pseudotime vs latent time

...enerates a metric known as latent time. I am having a very had time seeing the difference between latent time and standard pseudotime analysis (usually carried out by slingshot, monocle 3, or PAGA) and whether one of them is redundant or not.

文章 人类癌症数据库突变注释

cosmic70: the latest COSMIC database with somatic mutations from cancer and the frequency of occurence in each subtype of cancer. For more updated cosmic, see instructions below on how to make them.icgc21: International Cancer Genome Consortium version 21 mutations.nci60: NCI-60 human tumor cell lin...

问题 转录组KEGG富集分析出现问题。

... normal_vs_tumor --pvalueCutoff 1 --ann.db org.Hs.eg.db --organism hsa --idtype ENSEMBL --totype ENTREZID 输入标黄代码后,出现如下错误: [root@bdc6beb897b9  16:59:07 /work/zhaoran/6.enrich]#Rscript $scriptdir/enrichKEGG_pip.R --deg.file $workdir/5.deg/normal_vs_tumor.DEG.final.tsv ...

问题 circos绘制基因组共线性图绘图命令运行问题

...rence/best_practices   The debugging facility is helpful to figure out what's happening under the   hood       http://www.circos.ca/documentation/tutorials/configuration/debugging   If you're still stumped, get support in the Circos Google Group.       http://groups.google.com/...

文章 Monocle2|单细胞测序的拟时序分析

...单细胞测序文章《Single-cell profiling of tumor heterogeneity and the microenvironment in advanced non-small cell lung cancer》(IF=14.92)就多次采用拟时序分析,探讨不同细胞亚型之间的分化关系。 拟时序分析 为了方便大家能高效快捷的对单细胞数...

问题 GO注释:blast了nr数据库并通过blast2go获得了GO注释号;也通过blast了swiss-prot数据库获得了GO注释号。这两个结果需要合并吗?

...espectively, were subjected to GO analysis based on sequence alignment and the results were used in GO enrichment analysis. In total, 261 (C. exilicauda) and 197 (M. martensii) CDSs were assigned with 1562 (C. exilicauda) and 1390 (M. martensii) GO terms of knownfunction. The percentages of CDSs wit...

文章 multi_cox.r 基因表达量做多因素cox分析

...e [predict.time ...]                        Time point to draw the ROC curve [default 365 1095                        1825]  -c cut.score, --cut.score cut.score                        set cut score value to divide high and low risk groups                ...

文章 hclust_analysis.r 转录组数据层次聚类分析

...  输入基因表达矩阵文件,或者免疫侵润矩阵文件 cell_type TCGA-B7-A5TK-01A-12R-A36D-31 TCGA-BR-7959-01A-11R-2343-13 TCGA-IN-8462-01A-11R-2343-13 TCGA-BR-A4CR-01A-11R-A24K-31 TCGA-CG-4443-01A-01R-1157-13 aDC 0.612131 0.452721 0.434065 0.352635 0.268974 B ce...

文章 compare_stat_boxplot.r metadata中变量分组T检验比较并绘制box图

...p name from metadata to subset [default NULL]  -G geom, --geom geom  set type of plot: boxplot, violin, splitviolin                        [default=boxplot]  -p palette [palette ...], --palette palette [palette ...]                        fill color palette [default ('#e41...

文章 全基因组加倍/复制介绍(whole genome duplication, WGD)

...on,Horticultural Plant Journal,20192. Van de Peer Y, Maere S, Meyer A. The evolutionary significance of ancient genome duplications[J]. Nature Reviews Genetics, 2009, 10(10): 725-732.3. Badouin, H., Gouzy, J., Grassa, C.J., Murat, F., Staton, S.E., Cottret, L., Lelandais-Briere, C., Owens, G.L.,...

文章 ANNOVAR人类各个数据库变异注释结果表格说明

...DERATE or LOW impact).GENEGene name (usually HUGO)GENEIDGene ID)FEATUREThe type of feature is in the next field (e.g. transcript, motif, miRNA, etc.)FEATUREIDTranscript ID (preferably using version number), Motif ID, miRNA, ChipSeq peak, Histone mark, depending on the annotation.BIOTYPEDescription o...

文章 bioperl处理fasta和fastq序列截取反向互补翻译等等

...   -format => 'fasta');while ( my $seqobj = $in->next_seq() ) { # the human read-able id of the sequence my $id=$seqobj->id();  # string of sequence my $seq=$seqobj->seq(); # a description of the sequence my $desc=$seqobj->desc();  # one of 'dna','rna','protein'  https://www.bi...

文章 maf_oncoplot.r突变注释文件分类可视化maftools

...message and exit -i maffile, --maf maffile input the maf file[required] -m metadata, --meta metadata input metadata file path[required] -g group, --group group input group id in metadata file to classification ...