In Biojava, and I assume other Bio* projects? Indexes for sequences positions using the standard API are 1-based inclusive rather than the more standard 0-based exclusive on the top end. This means if you are extracting a subsequence up to a certain number of leading bases which match something else, say, you end up having to be careful , as you can't call .subSeqeuence(0,-1) which throws an exception rather than returning an empty sequence (as .substring(0,0) would).
Is there a particular reason for this design choice? It seems confusing as it is counter-intuitive to other programming indexing systems.
•
written
by
Nickengland
The choice of indexing is the source of quite a bit of contention - from a programming perspective 0 based, open end indexing is far more preferable. From a user interface perspective and when communicating with life scientists a 1 based inclusive indexing is absolutely required.
Being off by one is an endemic problem in bioinformatics and has probably already caused tens if not hundreds of millions of dollars in wasted resources (incorrect results etc).
Historically numerical oriented languages such as Fortran and for example R are one based. General programming oriented languages such as C, Java, Perl, Python are zero based. Personally I believe that Bio* projects should use one based indexing as it is the lesser of two evils, it seems the BioJava developers think the same way.
I believe that Bio* projects should use one based indexing as it is the lesser of two evils
Really? I'd say that since the bio* projects are for developers, 0-based indexing makes more sense. Then it's just a matter of converting to 1-based when showing to end-users.
That makes sense.
It is just that the difficulty of train life scientist to use programming concepts while littering the way with time bombs like this makes the task a lot more difficult than it needs to be. Having a zero based system is like starting with a wall instead of a mellow slope. I get occasionally tripped up when switching languages - and it is one of the most devious of errors as many times it remains hidden.
•
written
by
Istvan Albert
Log in to answer this question.
I don't think there really is a reason it chose to use 1-based arrays. Perhaps because gene sequences are 1-based inclusive and they wanted arrays to be reflective of that.
Really? Biopython sequences are definitely zero-based, NOT one-based. This is to follow Python conventions.