Bio.Affy.CelFile 模組
從 Affymetrix CEL 檔案版本 3 和 4 讀取資訊。
- exception Bio.Affy.CelFile.ParserError(*args)
繼承自:
ValueError
Affymetrix 解析器錯誤。
- __init__(*args)
初始化類別。
- class Bio.Affy.CelFile.Record
繼承自:
object
儲存 cel 檔案中的資訊。
使用範例
>>> from Bio.Affy import CelFile >>> with open("Affy/affy_v3_example.CEL") as handle: ... c = CelFile.read(handle) ... >>> print(c.ncols, c.nrows) 5 5 >>> print(c.intensities) [[ 234. 170. 22177. 164. 22104.] [ 188. 188. 21871. 168. 21883.] [ 188. 193. 21455. 198. 21300.] [ 188. 182. 21438. 188. 20945.] [ 193. 20370. 174. 20605. 168.]] >>> print(c.stdevs) [[ 24. 34.5 2669. 19.7 3661.2] [ 29.8 29.8 2795.9 67.9 2792.4] [ 29.8 88.7 2976.5 62. 2914.5] [ 29.8 76.2 2759.5 49.2 2762. ] [ 38.8 2611.8 26.6 2810.7 24.1]] >>> print(c.npix) [[25 25 25 25 25] [25 25 25 25 25] [25 25 25 25 25] [25 25 25 25 25] [25 25 25 25 25]]
- __init__()
初始化類別。
- Bio.Affy.CelFile.read(handle, version=None)
讀取 Affymetrix CEL 檔案並傳回 Record 物件。
支援 CEL 檔案格式版本 3 和 4。如果知道版本,請將 CEL 檔案格式指定為 3 或 4 作為 version 參數。如果未指定版本號碼,解析器將嘗試從檔案內容偵測版本。
此函式傳回的 Record 物件會將 CEL 檔案中的強度儲存在 record.intensities 中。目前,在解析版本 4 CEL 檔案時,不會設定 record.mask 和 record.outliers。
使用範例
>>> from Bio.Affy import CelFile >>> with open("Affy/affy_v3_example.CEL") as handle: ... record = CelFile.read(handle) ... >>> record.version == 3 True >>> print("%i by %i array" % record.intensities.shape) 5 by 5 array
>>> with open("Affy/affy_v4_example.CEL", "rb") as handle: ... record = CelFile.read(handle, version=4) ... >>> record.version == 4 True >>> print("%i by %i array" % record.intensities.shape) 5 by 5 array