how to get all iterations if I meet NaN?

Patrizio Frederic · 2009-03-27T21:59:33Z

2009/3/27 huiming song : > hi, everybody, please help me with this question: > > If I want to do iteration for 1000 times, however, for the 500th iteration, > there is NaN appears. Then the iteration will stop. If I don't want the stop > and want the all the 1000 iterations be done. What shall I do? > > > suppose I have x[1:1000] and z[1:1000],I want to do some calculation for all > x[1] to x[1000]. > > z=rep(0,1000) > for (i in 1:1000){ > ?z[i]=sin(1/x[i]) > } > > if x[90

Patrizio Frederic

Fri, Mar 27, 2009 2:59 PM

2009/3/27 huiming song <huimings at gmail.com>:

not sure I properly understood. Consider:

 x = seq(-pi,pi,length=1001)
 z = sin(1/x)
# Warning message:
# In sin(1/x) : NaNs produced
x[500:502]; z[500:502]
# [1] -0.006283185  0.000000000  0.006283185
# [1] -0.8754095        NaN  0.8754095

one NaN and one warning have been created, the remaining 1000
calculations has been executed.

lim(x->0)sin(1/x) not exists so sin(1/0) is not a number nan

z1 = z[!is.nan(z)]
x1 = x[!is.nan(z)]

# x and z without the z's nan position

hope that help

Patrizio

Thread (8 messages)

huiming song how to get all iterations if I meet NaN? Mar 27 Uwe Ligges how to get all iterations if I meet NaN? Mar 27 Patrizio Frederic how to get all iterations if I meet NaN? Mar 27 Nash how to get all iterations if I meet NaN? Mar 28 Uwe Ligges how to get all iterations if I meet NaN? Mar 28 Nash how to get all iterations if I meet NaN? Mar 29 Peter Dalgaard how to get all iterations if I meet NaN? Mar 29 Nash how to get all iterations if I meet NaN? Mar 29