Thursday, May 01, 2008

Generating a Affymatrix "tpmap" file for TiMAT

R Code that uses the Bioconductor module Affxparser.

Note: The code has a bug and in normal circumstances returns an error

> writeTpmap(filename="Dm_tiling2_MR_v01.tpmap", bpmaplist=bpmapr)
> Error in as.vector(y) : argument "y" is missing, with no default

The "y" can be fixed as follows

> y <- readBpmap("Dm_tiling2_MR_v01.bpmap")
> writeTpmap(filename="Dm_tiling2_MR_v01.tpmap", bpmaplist='y', verbose = as.integer(0))

This seems to work but doesn't print anything into the tpmap file due to some "sequence skipping".

I will try the Affymatrix windows based program next.

No comments: