Re: [R] Help with error in "if then" statement

From: Erik Iverson <iverson_at_biostat.wisc.edu>
Date: Fri, 11 Jul 2008 15:30:50 -0500

Hello -

Andrew Rominger wrote:
> Dear list,
>
> I'm afraid this is a mundane question. Here's the background: I've
> produced a function which allows me to sequentially measure angles and
> distances from a specified reference point to a curve defined by
> empirical data points while (i.e. using a while loop) the angle being
> measured is within certain bounds. Because the curve is circular I need
> to pars the data into basically an "upper" curve and "lower" curve. I
> tried to do this with an if statement, specifically:
>
> ycrit<-subset(data,subset=data$x==min(data$x)

You might consider not using 'data' as a variable name (it is the name of a function in the utils package) , and clearer code might be

ycrit <- subset(data, x == min(x))

But what does this return? It can definitely return a data.frame as in the following example:

t1 <- data.frame(a = rep(2,2))
t1
class(subset(t1, a == min(a)))

You probably do not want to them subsequently use this returned subset in an 'if' clause. Also, consider cases where there are NAs in the data:

t2 <- data.frame(a = c(NA,rep(2,2)))
t2
subset(t2, a == min(a))

This code returns a data.frame with 0 rows.

I don't know what your specific problem is, but hopefully this might lead you in the right direction.

> y.max<-length*sin(angle)+y.ref #length, angle and y.ref are given
> if(y.max<ycrit){
> ........
> }
>
> Now the problem: The while loop works for 4 iterations, until I get the
> error message:
>
> "Error in if (y.max < ycrit) { : missing value where TRUE/FALSE needed"
>
> When I forced the function to print y.max and ycrit for each iteration
> of the while loop, it returns finite real numbers, including ycrit =
> 153.5 and y.max = 245.16 for the step which returns the error message.
>
> Any ideas about what's going on here--why does R "think" that
> 245.16<153.5 is "missing," or is anything other than TRUE/FALSE? Am I
> using "if" incorrectly? In which case would it be more appropriate to
> perhaps create subsets of the data points based on < or > ycrit?
>
> Thanks in advance for any guidance--
> Andy
>



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. Received on Fri 11 Jul 2008 - 20:35:23 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 11 Jul 2008 - 21:32:24 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.

list of date sections of archive