Sometimes the following function call causes a database exception:> gene.seq <- getSequence (id=gene.map[,"ensembl_transcript_id"], type="ensembl_transcript_id",+ seqType="3utr", mart=hmart) I understand the above function must be called by try to capture the eventual error. WHat is not clear to me is how to realize that an error has occurred. The on-line documentation mentions an invisible object of class "try-error". How shall I test whetehr such object has been created or not ? I assume it is created whenever an error does occur ....? Thank you for your help. Maura tutti i telefonini TIM! [[alternative HTML version deleted]]
Hi, On Aug 3, 2009, at 11:53 AM, <mauede at alice.it> wrote:> Sometimes the following function call causes a database exception: > >> gene.seq <- getSequence (id=gene.map[,"ensembl_transcript_id"], >> type="ensembl_transcript_id", > + seqType="3utr", mart=hmart) > > I understand the above function must be called by try to capture the > eventual error. > WHat is not clear to me is how to realize that an error has > occurred. The on-line documentation > mentions an invisible object of class "try-error". > How shall I test whetehr such object has been created or not ? > I assume it is created whenever an error does occur ....?By "invisible" object, the documentation is saying that it's returning an object via the "invisible" function (see ?invisible). The object is no longer invisible if the function is called and the return value is assigned to a variable: which you are doing in your case. In your case, you just have to check to see if the object returned by getSequence is a (or inherits from) the "try-error" class. The last line of the help in ?try actually shows you how to do that. Here's a simpler example: R> m <- matrix(1:100, 10, 10) R> a <- try(m[,5]) # no error R> a [1] 41 42 43 44 45 46 47 48 49 50 R> b <- try(m[,11]) # error here b [1] "Error in try(m[, 11]) : subscript out of bounds\n" attr(,"class") To programmatically test if b is an error: R> is(b, 'try-error') [1] TRUE R> > inherits(b, 'try-error') # used in last line of ?try [1] TRUE -steve -- Steve Lianoglou Graduate Student: Computational Systems Biology | Memorial Sloan-Kettering Cancer Center | Weill Medical College of Cornell University Contact Info: http://cbio.mskcc.org/~lianos/contact
gotError <- try(gene.seq <- getSequence
(id=gene.map[,"ensembl_transcript_id"],
type="ensembl_transcript_id",
seqType="3utr", mart=hmart))
if (inherits(gotError, 'try-error')){.....error occurred, do
recovery........}
On Mon, Aug 3, 2009 at 11:53 AM, <mauede at alice.it>
wrote:> Sometimes the following function call causes a database exception:
>
>> ?gene.seq <- getSequence
(id=gene.map[,"ensembl_transcript_id"],
type="ensembl_transcript_id",
> + ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?seqType="3utr", mart=hmart)
>
> I understand the above function must be called by try to capture the
eventual error.
> WHat is not clear to me is how to realize that an error has occurred. The
on-line documentation
> mentions an invisible object of class "try-error".
> How shall I test whetehr such object has been created or not ?
> I assume it is created whenever an error does occur ....?
>
> Thank you for your help.
> Maura
>
>
> tutti i telefonini TIM!
>
>
> ? ? ? ?[[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
--
Jim Holtman
Cincinnati, OH
+1 513 646 9390
What is the problem that you are trying to solve?
Put your function in try() and the test it's class
gene.seq <- try(getSequence (id=gene.map[,"ensembl_transcript_id"],
type="ensembl_transcript_id", seqType="3utr", mart=hmart))
if(class(gene.seq) == "try-error"){
#code to run when an error occurs
} else {
#code to run when no error occurs
}
HTH,
Thierry
------------------------------------------------------------------------
----
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature
and Forest
Cel biometrie, methodologie en kwaliteitszorg / Section biometrics,
methodology and quality assurance
Gaverstraat 4
9500 Geraardsbergen
Belgium
tel. + 32 54/436 185
Thierry.Onkelinx at inbo.be
www.inbo.be
To call in the statistician after the experiment is done may be no more
than asking him to perform a post-mortem examination: he may be able to
say what the experiment died of.
~ Sir Ronald Aylmer Fisher
The plural of anecdote is not data.
~ Roger Brinner
The combination of some data and an aching desire for an answer does not
ensure that a reasonable answer can be extracted from a given body of
data.
~ John Tukey
-----Oorspronkelijk bericht-----
Van: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org]
Namens mauede at alice.it
Verzonden: maandag 3 augustus 2009 17:54
Aan: r-help at stat.math.ethz.ch
Onderwerp: [R] How to catch an error using try
Sometimes the following function call causes a database exception:
> gene.seq <- getSequence
(id=gene.map[,"ensembl_transcript_id"],
> type="ensembl_transcript_id",
+ seqType="3utr", mart=hmart)
I understand the above function must be called by try to capture the
eventual error.
WHat is not clear to me is how to realize that an error has occurred.
The on-line documentation mentions an invisible object of class
"try-error".
How shall I test whetehr such object has been created or not ?
I assume it is created whenever an error does occur ....?
Thank you for your help.
Maura
tutti i telefonini TIM!
[[alternative HTML version deleted]]
______________________________________________
R-help at r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Dit bericht en eventuele bijlagen geven enkel de visie van de schrijver weer
en binden het INBO onder geen enkel beding, zolang dit bericht niet bevestigd is
door een geldig ondertekend document. The views expressed in this message
and any annex are purely those of the writer and may not be regarded as stating
an official position of INBO, as long as the message is not confirmed by a duly
signed document.
Seemingly Similar Threads
- Conflict command getSequence {biomaRt} and getSequence {seqinr} !!
- Problem with biomaRt::getSequence.
- why biomaRt cannot extract 3UTR sequences for 1941 ENSGxxxxx ?
- Cannot install pakages from Bioconductor besides the default installation
- Fortran-90 and R