用docker,運(yùn)行rstudio-server,可以嘗試rocker/tidyverse
好物分享——R語言版本的bedtools不要重復(fù)造輪子 ——“哈嘍,能請教你一下怎樣在R server中調(diào)用Linux下的工具?”——“你為什么要這么做呢?”——“我想把Linux工具輸出的結(jié)果直接在R中分析,不想...
用docker,運(yùn)行rstudio-server,可以嘗試rocker/tidyverse
好物分享——R語言版本的bedtools不要重復(fù)造輪子 ——“哈嘍,能請教你一下怎樣在R server中調(diào)用Linux下的工具?”——“你為什么要這么做呢?”——“我想把Linux工具輸出的結(jié)果直接在R中分析,不想...
老教師除了使用PPT備課, 還會(huì)高頻率使用word進(jìn)行論文編寫、科研項(xiàng)目組編等等。 大家都知道論文存在以下特點(diǎn): 章節(jié)多、頁碼多、圖表多、文獻(xiàn)多、編號(hào)多 因此如果不掌握一些W...
"一個(gè)人能同時(shí)保有全然相反的兩種觀念,還能正常行事,是第一流智慧的標(biāo)志。" 到底什么才是深度思考? 我們都認(rèn)為提升深度思考能力在于見多識(shí)廣,知識(shí)積累多,自然水到渠成,有句話不...
對于 Mac pro 來說,連接服務(wù)器時(shí)經(jīng)常出現(xiàn)一種情況是 :ssh登錄沒問題,但是超過 5 分鐘不操作就會(huì)出現(xiàn)問題,沒法輸入,不得不重啟終端 或者 iTerm 。解決方法是...
箱線圖展示的就是分位數(shù),中間的線表示的是中位數(shù),也就是50%分位數(shù),如果非要在箱線圖上畫上表示平均值的線段也是可以實(shí)現(xiàn)的,今天介紹一下實(shí)現(xiàn)代碼 示例數(shù)據(jù)集我們用R語言的內(nèi)置數(shù)...
aplot Y叔開發(fā),拼圖解決坐標(biāo)對齊的問題 extrafont 解決ggplot2繪圖導(dǎo)出pdf文件中文字體亂碼的問題 對應(yīng)的B站介紹視頻https://www.bilib...
感謝作者的分享
生成geno2go那一步有點(diǎn)慢,用這個(gè)來生成geno2go的表會(huì)快很多
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)
生信 | 構(gòu)建物種的OrgDb寫在前面 本文主要參考生信小白2018[http://www.itdecent.cn/p/bb4281e6604e]與多啦A夢的時(shí)光機(jī)[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()
改了一下滑窗繪制的代碼,終端使用時(shí)用Rscript調(diào)用起來方便一點(diǎn)
生信學(xué)習(xí)筆記:使用SNP data做基因滲入分析 (4)繼續(xù)上一周的內(nèi)容,上次講到了如何使用--tree Dsuite的功能。要使用--treeDsuite的功能,我們顯然需要一個(gè)樹文件。一般樹文件可以通過SNP數(shù)據(jù)畫取進(jìn)化樹來獲...