????????去年寫的是二倍體的SSR數(shù)據(jù)轉(zhuǎn)0，1。代碼倒也不是不能跑單倍型，但是會存在多統(tǒng)計一個空值的情況，因此需要在過濾數(shù)據(jù)的時候加上過濾空值，這樣就可以正常統(tǒng)計。

????? ? 目的：本文的代碼是直接針對單倍體數(shù)據(jù)進(jìn)行處理，將單倍體的SSR分子標(biāo)記數(shù)據(jù)轉(zhuǎn)為的0，1數(shù)據(jù)矩陣。

? ? ????解決思路：使用R對數(shù)據(jù)進(jìn)行處理。原始數(shù)據(jù)讀取→分析分子標(biāo)記等位變異數(shù)n→將其轉(zhuǎn)化為0，1數(shù)據(jù)矩陣。

一、加載需要的R包及讀取數(shù)據(jù)

二、提取每一列的數(shù)據(jù)，進(jìn)行位點統(tǒng)計并轉(zhuǎn)化

這里輸出的文件在R中是：

txt格式是：

三、最后的數(shù)據(jù)轉(zhuǎn)化

最后是需要手動調(diào)節(jié)的部分，把樣品信息粘貼到最左側(cè)。

R中界面有一個這樣的統(tǒng)計，這里就是統(tǒng)計了每個引物的變異數(shù)。對應(yīng)到輸出的文件就是V1-V7表示的是SSR107的等位變異。V8-V26表示的是SSR125的等位變異。

【最復(fù)雜的代碼，做最簡單的事?！?/span>