找到约 15 条结果

文章 perl 中 grep 的使用

...perl里面的一个列表筛选函数,功能还是不错的。 grep有2种表达方式: 1. grep  BLOCK  LIST 2. grep  EXPR,  LIST BLOCK表示一个code块,通常用{}表示;EXPR表示一个表达式,通常是正则表达式。LIST是要匹配的列表。 示例: 统计匹配表达式...

问题 在NCBI网站上,同一个SAMN具有多个SRR文件,获得的fastq文件合并成一个文件

在NCBI网站上,同一个SAMN具有多个SRR文件,如何将这些下载下来的多个fastq文件合并成一个文件,并且最后获得read_1和read_2文件

问题 电脑内存不能升级,运行Rstudio内存一直上升,请问还有其他方法解决这个问题吗。

...ple_type <- "Primary solid Tumor"    > legacy <- FALSE > # 工作目录如果不存在,则创建目录 > if( !file.exists(work_dir) ){ +   if( !dir.create(work_dir, showWarnings = FALSE, recursive = TRUE) ){ +     stop(paste("dir.create failed: outdir=",work_dir,sep=""))...

文章 VirtualBox 中练习安装linux系统centos7

...好了:下载地址:https://www.virtualbox.org/wiki/Downloads 1.2 下载centos7系统iso文件: 海量的高通量测序数据,windows是处理不了的。因此果断选择linux服务器,开源免费而且生物数据分析软件也免费。这里选择CentOS版本的linux作为服...

问题 重复序列鉴定

...题,具体如下:运行这个脚本:perl $scriptsdir/ProtExcluder1.2/ProtExcluder.pl ${lib}_blast_results.txt ${lib},会显示ProtExcluder1.2/这个文件夹里的多个脚本有权限问题,可是我并没有运行有权限问题的脚本。其次,这些脚本为什么会存在权限问...

问题 基因家族分析提取启动子序列出现问题,望老师解答

...列时,biolinux系统运行脚本失败。 基因位置信息如下 XP_016539240.1    NC_029984.1    140133542    140145037    +XP_016549783.1    NC_029988.1    2729089    2737473    +XP_016550150.1    NC_029988.1    148442883    148468095    +XP_016553914.1 ...

文章 剑桥女大学生跳机自杀:学术研究遇到瓶颈到底怎么解决?

...豁然开朗。 3,换一个简单的方式 面对高强度的研究工作,遇到失败是正常的, 你可以在碰到特别棘手的事情时,同时去做一些更简单、更有可能得出积极结果的事情, 比如一个小小的实验、或者一次小长跑等等。 即使...

文章 单细胞转录组数据挖掘流程记录-头颈癌(HNSCC)(GSE181919)

...CD4+FOXP3+调节性T细胞在白斑中扩增并在胶原刺激下表达LAIR2,这可能为肿瘤进展提供有利的微环境。总之,基质和免疫细胞相互作用在 HNSCC 的逐步进展中发挥着关键作用,我们已经确定了 HNSCC 进展过程中可能的治疗干预措施的步...

问题 老师,您好,我在运行gatk的HaplotypeCaller时主脚本通过 SLURM 调度系统提交作业,使用 GNU Parallel 并行处理多个样本并调用 HaplotypeCaller 脚本,目前运行了约六七个小时,我的工作日志.err和.out都为空这是正常的吗?

文章 相关系数小结

...,不知道出自哪里,暂且一看: 0.1以下 没有关系 0.10-0.29 微弱关系 0.30-0.49 中等关系 0.50以上 强关联 pearson相关系数使用时需要注意的一些情况: XY联合分布符合正太分布 对异常点敏感 只有X和Y完全线性相关的时候,pearson系数...

文章 IT项目管理经验分享, 如何做好IT项目管理

...样资源情况下,良好的项目管理可以提高项目完成质量和工作效率、避免延期或者缩短项目完成时间。今天聊一些实战经验,略偏重 IT 开发项目,也是我们基于在线和远程办公的一点思考。 远程办公改变了组织和团队的形式...

问题 hmmer网站搜索蛋白序列中的结构域,得到几个domain,但是序列范围重叠。

...hits found是3即得到3个Thioredoxin domain: thioredoxin、thioredoxin_2、thioredoxin_9,但是序列范围重叠,72-174、85-172、64-173,是什么情况?实际这个蛋白序列中只有一个thioredoxin domain吗? 用NCBI-CDD搜索的结果类似,也是给出multi-dom,但是序列...

问题 转录组分析时 jaspar_2018_human.txt 文件是怎样下载的

问题 Rscript $scriptdir/gwas_manhattan_plot.r -i pvalue_$m.txt -F $FAI -T Sugar_$m -n Sugar_${m}_manhattan -c 1.3e-6绘制曼哈顿图时,输入这段代码后,报错,显示颜色数量不足,怎么解决

文章 picrust2用法

...iime/pick_de_novo_otus/rep_set/qiime_rep_set.fasta rep_set.farm -rf picrust2_outpicrust2_pipeline.py -s rep_set.fa -i $workdir/5.pick_otu_qiime/pick_de_novo_otus/otu_table.biom -o picrust2_out --processes 1  --in_traits COG,EC,KO,PFAM,TIGRFAMadd_descriptions.py -i picrust2_out/EC_metagenome_out/pre...