Dear UseR,

Based on exonmap, probeset.to.gene returns the gene (ensembl identifier) corrseponding to the submitted probeset identifier. Likewise, gene.to.probeset returns all probeset identifiers belonging to the submitted ensembl gene identifier. Hence the following 2 examples should result in a TRUE.

	"2376039" %in% gene.to.probeset(probeset.to.gene("2376039"))
	"2376044" %in% gene.to.probeset(probeset.to.gene("2376044"))

However, this is not the case. So, if I understand it correctly, these 2 probesets are linked to a certain gene but, if for that gene, one asks for the corresponding probesets, these 2 are not retrieved. Has this to do with the annotation quality of these two probesets? The first probeset contains only one probe, the second only three (instead of the standard 4). The other probesets for the same gene work fine. Or do I need the EST data, xmapExtras(), in addition to xmapDatabase()?

Thanks in advance for your help!

Kind regards,
An


	[[alternative HTML version deleted]]

