用docker，运行rstudio-server，可以尝试rocker/tidyverse
好物分享——R语言版本的bedtools不要重复造轮子 ——“哈喽，能请教你一下怎样在R server中调用Linux下的工具？”——“你为什么要这么做呢？”——“我想把Linux工具输出的结果直接在R中分析，不想...
发简信
IP属地：勃兰登堡
-
-
老师发表学术论文，需要学会哪些Word技巧?
老教师除了使用PPT备课， 还会高频率使用word进行论文编写、科研项目组编等等。 大家都知道论文存在以下特点： 章节多、页码多、图表多、文献多、编号多 因此如果不掌握一些W...
-
解决复杂问题的思维模型
"一个人能同时保有全然相反的两种观念，还能正常行事，是第一流智慧的标志。" 到底什么才是深度思考？ 我们都认为提升深度思考能力在于见多识广，知识积累多，自然水到渠成，有句话不...
-
Linux小技巧汇总对于 Mac pro 来说，连接服务器时经常出现一种情况是 ：ssh登录没问题，但是超过 5 分钟不操作就会出现问题，没法输入，不得不重启终端 或者 iTerm 。解决方法是...
-
R语言ggplot2做箱线图的时候如何添加表示平均值的线箱线图展示的就是分位数，中间的线表示的是中位数，也就是50%分位数，如果非要在箱线图上画上表示平均值的线段也是可以实现的，今天介绍一下实现代码 示例数据集我们用R语言的内置数...
-
-
可能会用到的R语言包~2aplot Y叔开发，拼图解决坐标对齐的问题 extrafont 解决ggplot2绘图导出pdf文件中文字体乱码的问题 对应的B站介绍视频https://www.bilib...
-
感谢作者的分享
生成geno2go那一步有点慢，用这个来生成geno2go的表会快很多
gos_list <- function(x){
the_gos <- str_split(x[2], ",", simplify = FALSE)[[1]]
df_temp <- data.frame(GID = rep(x[1], length(the_gos)),
GO = the_gos,
EVIDENCE = rep("IEA", length(the_gos)))
return(df_temp)
}
gene2gol <- apply(as.matrix(gos),1,gos_list)
gene2gol_df <- do.call(rbind.data.frame, gene2gol)生信 | 构建物种的OrgDb
写在前面 本文主要参考生信小白2018[https://www.jianshu.com/p/bb4281e6604e]与多啦A梦的时光机[https://www.jiansh...
-
args <- commandArgs(trailingOnly=TRUE)
table <- read.table(args[1], header=T)
windowCenter=(table$windowStart+table$windowEnd)/2
pdf(args[2], height=7, width=7)
plot(windowCenter, table$D, type="l", xlab="Position", ylab="D (gray) / fD (black)", ylim=c(0,1), main="Chromosome 5 (altfas, neocan, telvit)", col="gray")
lines(windowCenter, table$f_d)
dev.off()
改了一下滑窗绘制的代码，终端使用时用Rscript调用起来方便一点生信学习笔记：使用SNP data做基因渗入分析 (4)
继续上一周的内容，上次讲到了如何使用--tree Dsuite的功能。要使用--treeDsuite的功能，我们显然需要一个树文件。一般树文件可以通过SNP数据画取进化树来获...
-
作图
-
收藏数量遗传学和生物统计学知识
-
遗传学
