用docker，運行rstudio-server骤素，可以嘗試rocker/tidyverse
好物分享——R語言版本的bedtools不要重復造輪子 ——“哈嘍家淤，能請教你一下怎樣在R server中調(diào)用Linux下的工具异剥？”——“你為什么要這么做呢？”——“我想把Linux工具輸出的結(jié)果直接在R中分析絮重，不想...
發(fā)簡信
IP屬地：勃蘭登堡
-
-
老師發(fā)表學術(shù)論文冤寿，需要學會哪些Word技巧?
老教師除了使用PPT備課， 還會高頻率使用word進行論文編寫青伤、科研項目組編等等疚沐。 大家都知道論文存在以下特點： 章節(jié)多、頁碼多潮模、圖表多、文獻多痴施、編號多 因此如果不掌握一些W...
-
解決復雜問題的思維模型
"一個人能同時保有全然相反的兩種觀念，還能正常行事羡玛，是第一流智慧的標志别智。" 到底什么才是深度思考？ 我們都認為提升深度思考能力在于見多識廣稼稿，知識積累多薄榛，自然水到渠成，有句話不...
-
Linux小技巧匯總對于 Mac pro 來說，連接服務器時經(jīng)常出現(xiàn)一種情況是 ：ssh登錄沒問題谋右，但是超過 5 分鐘不操作就會出現(xiàn)問題硬猫，沒法輸入，不得不重啟終端 或者 iTerm 倚评。解決方法是...
-
R語言ggplot2做箱線圖的時候如何添加表示平均值的線箱線圖展示的就是分位數(shù)浦徊，中間的線表示的是中位數(shù)，也就是50%分位數(shù)天梧，如果非要在箱線圖上畫上表示平均值的線段也是可以實現(xiàn)的盔性，今天介紹一下實現(xiàn)代碼 示例數(shù)據(jù)集我們用R語言的內(nèi)置數(shù)...
-
-
可能會用到的R語言包~2aplot Y叔開發(fā)冕香，拼圖解決坐標對齊的問題 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)生信 | 構(gòu)建物種的OrgDb
寫在前面 本文主要參考生信小白2018[http://www.reibang.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調(diào)用起來方便一點生信學習筆記：使用SNP data做基因滲入分析 (4)
繼續(xù)上一周的內(nèi)容弊予，上次講到了如何使用--tree Dsuite的功能。要使用--treeDsuite的功能开财，我們顯然需要一個樹文件汉柒。一般樹文件可以通過SNP數(shù)據(jù)畫取進化樹來獲...
-
作圖
-
收藏數(shù)量遺傳學和生物統(tǒng)計學知識
-
遺傳學
