????R剛?cè)腴T的時(shí)候,能夠正確讀取單個(gè)文件就覺得小有成就,隨著時(shí)間的積累,單一文件地讀取已經(jīng)不能滿足需求了,此時(shí),批量地做就是解放雙手地過程。 使用for循環(huán)把下載地TCG...
????R剛?cè)腴T的時(shí)候,能夠正確讀取單個(gè)文件就覺得小有成就,隨著時(shí)間的積累,單一文件地讀取已經(jīng)不能滿足需求了,此時(shí),批量地做就是解放雙手地過程。 使用for循環(huán)把下載地TCG...
1. Kaplan-Meier生存分析 Kaplan Meier是一種單因素生存分析,它可用于研究1個(gè)因素對(duì)于生存時(shí)間的影響。參考:生存分析之Kaplan-Meier曲線都告...
三個(gè)月前R 4.0版本就已經(jīng)正式上線了,之前怕麻煩就一直沒更新,而現(xiàn)在由于要用到新版本的R包,才想起順便一起更新了R。 目前在官網(wǎng)上的最新版R是2020-06-22的4.0....
list轉(zhuǎn)化為data.frame ??在數(shù)據(jù)處理過程中會(huì)遇到各式各樣的要求,格式轉(zhuǎn)換、數(shù)據(jù)結(jié)構(gòu)重構(gòu)等需求可以說是家常便飯了,所以R語言本身也提供了很多數(shù)據(jù)格式轉(zhuǎn)換的函數(shù)或者...
一、4Dtv(四倍簡(jiǎn)并位點(diǎn)顛換率)的定義 4DTV stands for four-fold synonymous (degenerative) third-codon tr...
DepthOfCoverage會(huì)輸出7個(gè)文本結(jié)果。其中一個(gè)是按照interval上的每個(gè)堿基,輸出一行統(tǒng)計(jì)信息,所以會(huì)比較慢: .DepthOfCoverage.txt結(jié)果 ...
上一篇講了GenVisR 基因組數(shù)據(jù)可視化實(shí)戰(zhàn) (一)[http://www.itdecent.cn/p/bb695b053246]是GenVisR畫waterfall 圖...
上次我們整理到bwa比對(duì)后得到bam文件,下一步我們要通過GATK流程從bam文件中call variant。一、使用GATK前須知事項(xiàng): 對(duì)GATK的測(cè)試主要使用的是人類全...
首先是加載相關(guān)的包,mutate主要屬于dplyr包里,這里我們統(tǒng)一使用tidyverse包。tidyverse包中含有各種數(shù)據(jù)整理以及畫圖的包,如下加載tidyverse包...