找到约 15 条结果

问题 在RStudio上用setwd()命令未能修改数据资料路径

...是运行不成功,显示: >  setwd("D:/20251124_data") Error in setwd("D:/20251124_data") : cannot change working directory >  setwd("D:/20251124data") Error in setwd("D:/20251124data") : cannot change working directory > 为什么? 用dir("D:/"),显示如下: >...

文章 将fasta序列多行变单行——使用 awk 或 seqkit 或 python

...快) awk '{if($0~/>/) name=$0 ;else seq[name]=seq[name]$0;}END{for(i in seq) print i"\n"seq[i]}}' test.fa>seq1ATTAAAGGTTTATACCTTCCAATCTGTGTGGCTGTCACTCACACGAGTAACTCGTCTATC>seq2ATTAAAGGTTTATACCTTCCAATCTGTGTGGCTGTCACTC 方法二:seqkit工具 将多行序列转换为一行序列 seqki...

文章 Cufflinks,Stringtie 合并转录本之后,如何筛选新转录本?

...a reference transcript4eSingle exon transfrag overlapping a reference exon and at least 10 bp of a reference intron, indicating a possible pre-mRNA fragment.5iA transfrag falling entirely within a reference intron6oGeneric exonic overlap with a reference transcript7pPossible polymerase run-on fragme...

文章 2020年-基因家族分析文章怎么发!

... Liu Y, Cui J, Zhou X, et al. Genome-wide identification, characterization and expression analysis of the TLP gene family in melon (Cucumis melo L.)[J]. Genomics, 2020. 重点福利来了! 组学大讲堂专注于生物信息培训,已有20000+的线上网课学员和1000+的线下培训学员...

问题 老师您好,在KAKS分析时duplication_gene.out导出为空怎么解决呢?

在运行perl ../script/KAKS_SHAIXUAN.pl -in1 JmjC_cds_confirmed.fa.fai -in2 JmjC_cds_confirmed_blast.out -out duplication_gene.out这一步后结果为空,但中途也没有报错。 以下是生成的文件,还想麻烦老师看看是什么问题。

文章 uniref50 数据库下载及格式化

...t_taxid_uniref.pl division.dmp nodes.dmp uniref90.fasta.gz ./ for i in PLN MAM INV PHG PRI ROD SYN UNA VRL VRT ENV BCT ;do diamond makedb --in ${i}_uniref90.fa -d  ${i}_uniref90.fa done

问题 在进行KAKS计算时出现不报错,但结果文件是空的现象

...or -i 2.axt -o kaks2.txt<br>Method(s): MA<br>Genetic code: 1-Standard Code<br>Please wait while reading sequences and calculating...<br>manager@bl8vbox[dupl-id1]  这是无结果文件输出的显示! manager@bl8vbox[dupl-id1] /home/manager6/KaKs_Calculator2.0/KaKs_Calcu...

文章 总结:生化检测实验样本值不准的原因分析

...居然有11分的文献引用! Lee H Y, Back K. Melatonin induction and its role in high light stress tolerance in Arabidopsis thaliana [J]. Journal of pineal research, 2018. 这个J Pineal Res期刊,2018年IF=11.613,也是很牛了!! 学霸姐姐大致看了一下,这篇文章主...

文章 R语言绘制多分组差异分析结果图

...atial transcriptomics analysis of esophageal squamous precancerous lesions and their progression to esophageal cancer. Nat Commun 14, 4779 (2023).

文章 微生物多样性alpha分析看不明白?看这里!

...用来估计群落中含有OTU 数目的指数,同样由Chao提出(Chao and Yang, 1993),是生态学中估计物种总数的常用指数之一。默认将序列量10以下的OTU都计算在内,从而估计群落中实际存在的物种数。ACE指数越大,表明群落的丰富度越高。 ...

文章 Rocky Linux安装 时报错: DNF error

Rocky Linux安装 时报错: DNF error : Error in POSTTRANS scriptlet in rpm package kernel-core 解决办法: 重新 安装 ,在手动分区时,将/boot设置为300M以上。 此外,我们在网易云课堂上有各种教学视频,有兴趣可以了解一下:...

文章 seurat_sc_cluster.r 脚本使用帮助:

...mloptional arguments:  -h, --help            show this help message and exit  -i count, --count count                        input read count file , csv or tsv format[ optional]  --sep sep             the field separator character for count file; Values            ...

文章 做微生物研究必懂的OTU table相关知识

...格,后续的分析就很好展开了。 参考文献:[1] T. Magoc and S. Salzberg. FLASH: Fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27:21 (2011), 2957-63.[2] Patrick D. Schloss,et al. Introducing mothur: Open-Source, Platform-Independent, Community-...

文章 TCGA数据下载—TCGAbiolinks包参数详解

...access to an unmodified copy of data that was previously stored in CGHub and in the TCGA Data Portal hosted by the TCGA Data Coordinating Center (DCC), in which uses as references GRCh37 (hg19) and GRCh36 (hg18). GDC harmonized database: data available was harmonized against GRCh38 (hg38) using G...

问题 16S 差异比较使用compare_errorbar.r的时候报错

...表我的分组信息这些都是准确的 这个代码不报错 for i in {2..6};do     Rscript $scriptdir/meta_stat.r -i ../05.export_data/taxa_summary/absolute_abundance/feature_table_tax_L$i.txt \ -m $metadata -t Group \ -a WSD_Lig_PPI -b WSD_PPI \ -o meta_stat --name WSD_Ligatu...