Bio.PopGen.GenePop 套件

子模組

模組內容

用於處理 GenePop 的程式碼。

請參閱 http://wbiomed.curtin.edu.au/genepop/,此格式在此處有說明文件:http://wbiomed.curtin.edu.au/genepop/help_input.html

類別:Record 保留 GenePop 資料。

函式:read 將 GenePop 記錄(檔案)解析為 Record 物件。

部分靈感來自 MedLine 程式碼。

Bio.PopGen.GenePop.get_indiv(line)

提取該行上個體資訊的詳細資料。

Bio.PopGen.GenePop.read(handle)

解析包含 GenePop 檔案的 handle。

handle 是一個類似檔案的物件,其中包含 GenePop 記錄。

class Bio.PopGen.GenePop.Record

基底類別:object

保留來自 GenePop 記錄的資訊。

成員

  • marker_len 標記長度(每個等位基因 2 或 3 位數字代碼)。

  • comment_line 註解行。

  • loci_list 基因座名稱的清單。

  • pop_list 族群名稱的清單。

  • populations 族群資料的清單。

在大多數 genepop 檔案中,族群名稱並不可信。強烈建議依索引參照族群。

populations 每個族群都有一個元素。每個元素本身都是個體清單,每個個體都是一對,由個體名稱和等位基因清單組成(每個標記 2 個,或單倍體 1 個):範例

[
    [
        ('Ind1', [(1,2),    (3,3), (200,201)],
        ('Ind2', [(2,None), (3,3), (None,None)],
    ],
    [
        ('Other1', [(1,1),  (4,3), (200,200)],
    ]
]
__init__()

初始化類別。

__str__()

傳回(重建)GenePop 文字表示。

split_in_pops(pop_names)

將 GP 記錄分割成每個條目 1 個族群的字典。

給定一個具有 n 個族群和 m 個基因座的記錄,傳回記錄字典(鍵為 pop_name),其中每個項目都是具有單一族群和 m 個基因座的記錄。

引數:- pop_names - 族群名稱

split_in_loci(gp)

將 GP 記錄分割成每個條目 1 個基因座的字典。

給定一個具有 n 個族群和 m 個基因座的記錄,傳回記錄字典(鍵為基因座名稱),其中每個項目都是具有單一基因座和 n 個族群的記錄。

remove_population(pos)

移除族群(按位置)。

remove_locus_by_position(pos)

按位置移除基因座。

remove_locus_by_name(name)

按名稱移除基因座。