I am processing the raw data on the CMAP where there are a lot of row data using the .CEL format. I know I can read the value using the R software. But ,what I want is to use python to process .CEL data. I have learned that there is a package called biopython which can process CEL data. Could anyone know the detail of how to process .CEL data using python?The following is my code to process .CEL data using python. But there is something wrong.
from Bio.Affy import CelFile
with open('AGENT_p_NCLE_RNA6_HG-U133_Plus_2_A01_436578.CEL') as handle:
c = CelFile.read(handle)
The result is as the followings:
<Bio.Affy.CelFile.Record object at 0x02534730>
What is the wrong with my code? And using R, the CDF is used however in python it is not used,why?
It would be nice of you to answer my problem.