[BioC] restrict transcriptDB object to "known" genes

[Stefanie]

Hi,

I just managed to do it the following way:

I get the transcript IDs of all known genes by the following command:

ensembl = useDataset("hsapiens_gene_ensembl",mart=ensembl)

res = getBM(attributes = "ensembl_transcript_id", filters = "status", values =
"known", mart = ensembl)

humanDb = makeTranscriptDbFromUCSC(genome = "hg19", tablename =
"ensGene",transcript_ids = as.character(res[,1]))

tx = transcriptsBy(humanDb)


Works! :)