定閱號

你好，請問一下崇决，這個bedgraph文件第二列和第三列是起始位置和終止位置嗎材诽？ 哪個代表C，哪個代表G呢恒傻？
以及這個bedgraph是0-based還是1-based呀脸侥？如果是0-based的，是不是說chr1 10496 10497 94 16 1中：
chr1：10497是C甲基化位點呢盈厘？

其次睁枕，我有點疑問的地方是：
chr1 10524 10525 88 15 2
chr1 10525 10526 86 13 2

這里第一行的chr1:10524 -10525如果代表10525處是一個C堿基，
那么第二行的chr1:10525-10526也代表10526處是一個C堿基，加上C后面會跟一個G譬重，那么從10525-10527不就變成了CCG了嗎拒逮？ 是可以這樣理解嗎？這樣還滿足第一行的CG結構嗎臀规？

還是說這里不區(qū)分正負鏈滩援，chr1:10524 -10525-10526是正鏈xCGx，然后chr1:10524-10525 -10526-10527是負鏈xGCx

這樣在做callDML的時候塔嬉，pos:10525需要考慮正負鏈嗎玩徊？ 即把
chr1 10524 10525 88 15 2
chr1 10525 10526 86 13 2
這兩行的beta信息考慮進去當成一個位點呢，但是我看up主你的代碼pos = T1[,3]，只考慮了后面這列，所以是不是當成10525和10526兩個位點進行計算了读第？

這里我感到有些疑惑，困擾了好久了畔塔，請up主指教

R語言 -- 尋找差異甲基化區(qū)域（DMR）-- DSS 包

最好的文檔其實還是官方文檔。鸯屿。澈吨。http://bioconductor.org/packages/release/bioc/vignettes/DSS/inst/doc/DS...

生信擺渡

9371 8 10

定閱號

沒有看到GQ婶恼，只有GERMQ桑阶，這個是GQ嘛？

基因突變檢測之Mutect2

寫在前面的話 今天是2021年1月3號勾邦，年底一番瞎忙活蚣录，停更了好長一段時間的簡書。人生中有個很重要伯樂導師告訴我眷篇，在這個事多的工作環(huán)境下萎河，沒人愿意看你太多的廢話。直奔重點才是...

靜小沐

24878 6 41

定閱號

接著上一次內(nèi)容繼續(xù)介紹铅歼，BBTools、BBMap中的一些實用的小工具换可。 BBMap Read Merger 合并雙末端（PE）reads,在預期中這些reads有重疊的位置...

定閱號

想問一下椎椰，就是一個qseqid對應多個taxid的這種，應該怎么處理呢沾鳄？ 應該給它歸為哪個物種呢慨飘？

宏基因組之物種注釋（基于nr庫）

昨天下午搗鼓了一下宏基因組物種注釋過程(基于nr庫)，現(xiàn)在將整個流程記錄一下。軟件需求：blast瓤的，diamond休弃，taxonkit(安裝自行百度) 構建細菌子庫 blast...

凱凱何_Boy

14608 14 10

定閱號

樓主討論的這個是不是目標序列是同一條染色體的情況啊，如果目標序列是不同物種的不同序列圈膏，而這個指標都一樣的情況該怎么辦呢塔猾？

我是想看自己的序列與哪個物種最相似，來判斷有沒有目標菌種存在稽坤。這種情況怎么辦呀丈甸？ 請陳老師指教

「JCVI」如何篩選得到最優(yōu)blast比對結果？

JCVI尿褪，包含了太多的功能睦擂，但是我感覺好像又沒有一個特別好的說明文檔（小聲bb，感謝唐老師的開發(fā)的好用工具） blast比對 未過濾的blast比對結果杖玲，所使用參數(shù)是：-ou...

陳有樸

2524 1 11

定閱號

你好顿仇，請問，如果只想返回最好的一個序列結果摆马，比如query3和query1只保留一條最好的臼闻。有沒有參數(shù)設置啊，還是只能自己寫代碼篩選今膊。

如果用一致性些阅，覆蓋度，比對得分進行排序后斑唬，依舊還是有很多序列的排名一樣市埋，比如100%匹配，得分都是278. 那么應該再怎么進行篩選呢恕刘？ 我是想看我的序列有多少條比對到了某個細菌上缤谎。我希望一個序列只有一個物種結果，這樣就方便后續(xù)統(tǒng)計褐着，這種應該怎么辦呀坷澡？請洲更老師指教

如何讓BLAST返回最優(yōu)的一個搜索結果，看看你沒有有進坑

大部分時候含蓉，我們都是看著別人的教程频敛，然后嘗試處理自己的數(shù)據(jù)，結果跑完了馅扣，如果和預期相符合就不會懷疑這個工具有啥問題斟赚。如果你要學習生物信息學，那么有一個信條一定要記住差油，不要盲目...

xuzhougeng

11177 4 29

定閱號

你好拗军，我想請問任洞，我star比對完的bam文件，我想提取沒有比對上的序列发侵，但是samtools view -f 4那個命令不行交掏，不知道怎么提取了。samtools flagstat unmapped_reads.bam說我百分百比對上了刃鳄，但是實際上我這個樣本的比對率只有16%盅弛，log文件里提示，短序列比對失敗的有40%多铲汪。

比對軟件STAR的使用

在之前的學習和練習里熊尉，比對這一步我使用過bowtie2（DNA比對）和hisat2（RNA-seq比對），現(xiàn)在學習另一個很火的軟件：STAR掌腰。STAR能夠發(fā)現(xiàn)非典型拼接和嵌合...

生信start_site

43922 19 50

定閱號

你好呀狰住，之前的推送看不到了，有其他平臺的齿梁，比如公眾號嗎催植？如何篩選感興趣的GO和KEGG進行繪圖~這個推送

KEGG數(shù)據(jù)庫 | 通路結果分類展示

在做完KEGG通路富集之后，紛繁復雜的通路映入眼簾勺择，略微有些雜亂無章创南，缺乏一個有序的排列，因此省核，本推送依據(jù)KEGG一級類目稿辙，將不同的通路歸納總結，分類排列气忠，使結果更加規(guī)整邻储，幫...

生命數(shù)據(jù)科學

1578 4 10

定閱號

您好，想咨詢一下旧噪，在做交叉驗證之前吨娜，需要先挑選好較少的特征嗎？即最終的模型特征淘钟？然后再進行交叉驗證宦赠？

56-caret包學習：模型訓練與調(diào)優(yōu)

1、模型訓練與參數(shù)優(yōu)化 在進行建模時米母，需對模型的參數(shù)進行優(yōu)化勾扭，在caret包中其主要函數(shù)是train。一旦定義了模型和調(diào)優(yōu)參數(shù)值铁瞒，就應該指定重采樣的類型妙色。目前，k折交叉驗證重...

wonphen

4544 4 11

定閱號

老師您好精拟，我有一個問題燎斩，在選到13作為最佳超參數(shù)后，是不是需要回到所有訓練集中蜂绎，再使用這個超參數(shù)后訓練一個分類模型栅表。由于超參數(shù)是13，是不是就是意味著保留13個特征對于這個模型的分類效果最好呀师枣。這個時候是不是選擇重要性前13個特征作為最終的特征怪瓶，然后預測測試集呢？因為我不會python践美，然后這個k鄰近模型也不懂洗贰，用的是R語言和隨機森林的理解方式理解的，不知道有沒有理解錯誤陨倡，請老師指正

單一驗證敛滋、k折交叉驗證（特例：留一法（LOOCV））、交叉驗證確定最優(yōu)超參數(shù)

一兴革、單一驗證 1绎晃、單一訓練集和測試集 最簡單的樣本集劃分就是只有訓練集和測試集，而沒有驗證集杂曲，因此無法利用驗證集反過來對模型參數(shù)進行調(diào)整庶艾。只能先給定一組超參數(shù)C，然后訓練得到...

馬爾代夫Maldives

11061 2 6

定閱號

數(shù)據(jù)在人工智能技術里是非常重要的！本篇文章將詳細給大家介紹3種數(shù)據(jù)集：訓練集棚饵、驗證集煤裙、測試集。同時還會介紹如何更合理的講數(shù)據(jù)劃分為3種數(shù)據(jù)集蟹地。最后給大家介紹一種充分利用有限數(shù)...

定閱號

1998 -- 2002 | 0歲 -- 4歲 話說我也不是從出生就開始記事兒的呀积暖，前幾年關于我的事情只能靠我家人跟我講述的了，今年過年要好好采訪我的家人關于我幼兒時對...

定閱號

今天在看碩士與博士的區(qū)別時看到下面一段話： 對博士而言夺刑，最大的挑戰(zhàn)其實是克服自我厭惡，然而在讀碩士時很少會有這種情緒的體現(xiàn)分别。太多焦慮的博士遍愿，并非沒有學術能力，而是無法克服這種...

定閱號

您好耘斩，請問一下沼填，您下載mixcr的時候，使用的時候括授，它有要許可證嗎坞笙？我這邊下載了3的版本和4的版本都需要有許可證岩饼，不知道怎么回事

mixcr3.0?軟件使用

MiXCR是一個通用的框架，處理從原始測序數(shù)據(jù)到定量clonotype的免疫組數(shù)據(jù)薛夜。 支持PE和SE的reads籍茧，考慮測序質(zhì)量，糾正PCR錯誤梯澜，還能識別生殖系高突變寞冯。 支持部...

11的霧

3253 4 4