[R] priority of operators in the FOR ( ) statement
Gabor Grothendieck
ggrothendieck at gmail.com
Tue Aug 23 13:28:46 CEST 2005
On 8/23/05, Ravi.Vishnu at outokumpu.com <Ravi.Vishnu at outokumpu.com> wrote:
> Dear All,
> I spent an entire evening in debugging a small, fairly simple program in R
> - without success. It was my Guru in Bayesian Analysis, Thomas Fridtjof,
> who was able to diagonose the problem. He said that it took a long time
> for him also to locate the problem.
> This program illustrates in some ways the shortcomings of the error
> messages that R responds with. In this case, it was quite misleading and
> directs attention to a location far removed the actual problem statement.
> Without any more introductory comments, let me directly discuss the
> essential problem. I am enclosing the entire program after a brief
> discussion.
>
> The problem arises from the following statement (nr is an integer
> constant) :
> for ( i in 1:nr-1) {.......}
> The unexpected problem (at least for me) is that R reads the above
> statement as (i in (1:nr)-1) {.....}. This makes i be initially as zero
> which leads to an error because the for loop in R starts from 1. The
> problem is easily fixed by writing the for loop as ( i in 1:(nr-1))
> {.......}. This would be an easy problem to fix if R directly indicates
> what the problem is. Instead, it gives mystifying error messages which are
> totally misleading. For example, to the program given below, I got the
> following error message (these point to commands elsewhere in the program)
> :
> Error in if ((x >= 0) & (x < s2)) return(x/2) else if ((x >= s2) & (x < :
>
> missing value where TRUE/FALSE needed
>
> I would like clarifications on the following points :
> 1. I am just curious to know if the priority of operators in the for
> statement ( the colon before the minus operator, for example) is a
> deliberate design decision. I have tested Matlab and found that it
> interprets my original statement correctly without an extra paranthesis.
?Syntax gives the operator precedence.
Also, note that : is probably best not used in functions since it does
not handle boundary conditions properly. If n were 0 then 1:n
results in two iterations corresonding to 1 and 0 but what you
really wanted was likely no iterations at all. To do that you need
seq(length = n) rather than ":".
Also I have found expressions like 0:1/10 handy to generate
0, .1, .2, ..., 1 and that works with the current precedence.
More information about the R-help
mailing list