[BioC] Problem with biomaRt::getSequence.
Tanvir Ahamed [guest]
guest at bioconductor.org
Tue May 7 12:54:17 CEST 2013
Hi,
I can run the code some days ago . But cant run now.
Problem 1: Output is ok
ensembl = useDataset("hsapiens_gene_ensembl",mart=ensembl)
utr5 = getSequence(chromosome=3, start=185514033, end=185535839, type="entrezgene",seqType="5utr", mart=ensembl)
Output :
5utr entrezgene
Sequence unavailable 10644
GGAGCGCCGGGTACCGGGCCGGGGGAGCCGCGGGCTCTCGGGGAAGAGACGG 10644
3 GGGGGGCGGAGGAGGAGGAGAGACGAGGGCAGCGGAGGAGGCGAGGAGCGCCGGGTACCGGGCCGGGGGAGCCGCGGGCTCTCGGGGAAGAGACGG 10644
CGGAGGAGGCGAGGAGCGCCGGGTACCGGGCCGGGGGAGCCGCGGGCTCTCGGGGAAGAGACGG 10644
No UTR is annotated for this transcript 10644
Problem 2:Problem is here
protein = getSequence(id=c(100, 5728),type="entrezgene",seqType="peptide", mart=ensembl)
Error in getBM(c(seqType, type), filters = type, values = id, mart = mart, :
Query ERROR: caught BioMart::Exception::Database: Error during query execution: Can't create/write to file '/mnt/ephemeral0/mysqltmp/#sql_40a_0.MYI' (Errcode: 2)
I need help please.
-- output of sessionInfo():
R version 3.0.0 (2013-04-03)
Platform: x86_64-w64-mingw32/x64 (64-bit)
locale:
[1] LC_COLLATE=Bangla_Bangladesh.1252 LC_CTYPE=Bangla_Bangladesh.1252 LC_MONETARY=Bangla_Bangladesh.1252 LC_NUMERIC=C
[5] LC_TIME=Bangla_Bangladesh.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] seqinr_3.0-7 biomaRt_2.16.0
loaded via a namespace (and not attached):
[1] RCurl_1.95-4.1 tools_3.0.0 XML_3.96-1.1
--
Sent via the guest posting facility at bioconductor.org.
More information about the Bioconductor
mailing list