search for: nan

...------------ Illustration -------------------- R : Copyright 2003, The R Development Core Team Version 1.8.0 (2003-10-08) > A <-matrix(0, 5, 5) > i <- c() > a <- A[i, ] ; crossprod(a) [,1] [,2] [,3] [,4] [,5] [1,] 6.578187e-313 NaN NaN NaN NaN [2,] NaN 1.273197e-313 NaN 1.485397e-313 NaN [3,] NaN 4.243992e-313 2.121996e-314 NaN NaN [4,] NaN 1.697597e-313 NaN 4.880590e-313 NaN [5,] 5.941588e-313 NaN NaN 1.697597e-313 NaN &g...

COD.EST VARIABLE AÑO DIA ENERO FEBRERO MARZO ABRIL MAYO JUNIO JULIO AGOSTO SEPBRE OCTUBRE NOVBRE DICBRE 1 2501001 PT 1986 1 NaN NaN NaN NaN 2.0 .0 .0 9.0 NaN NaN NaN NaN 2 2501001 PT 1986 2 NaN NaN NaN NaN .0 11.0 19.0 .0 NaN NaN NaN NaN 3 2501001 PT 1986 3 NaN NaN NaN NaN 11.0 3.0 .0 .0 NaN NaN NaN NaN 4 2501001 PT 1986...

Hi all... I built a matrix binding vectors with rbind, and have something like this: [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] CLS 3.877328 4.087636 4.72089 4.038361 3.402942 2.786285 2.671222 3.276419 ORD NaN NaN NaN NaN 5.770780 5.901113 11.888054 7.934823 FAM NaN NaN NaN NaN NaN 3.699455 4.551196 2.885390 GEN NaN NaN NaN NaN NaN 3.967411 4.390296 2.885390 SPP NaN NaN NaN NaN NaN NaN NaN 2.885390 Then I...

I'm having trouble grokking complex NaN's. This first set examples using complex(re=NaN,im=NaN) give what I expect > Re(complex(re=NaN, im=NaN)) [1] NaN > Im(complex(re=NaN, im=NaN)) [1] NaN > Arg(complex(re=NaN, im=NaN)) [1] NaN > Mod(complex(re=NaN, im=NaN)) [1] NaN > abs(complex(re=NaN, im=NaN))...

Dear all, I am struggling to understand this. What happens when you raise a negative value to a power and the result is a very large number? B [1] 47.73092 > -51^B [1] -3.190824e+81 # seems fine # now this: > x <- seq(-51,-49,length=100) > x^B [1] NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN NaN <snip> > is.numeric(x^B) [1] TRUE > is.real(x^B) [1] TRUE > is.infinite(x^B) [1] FALSE FALSE FALSE FALSE FALSE I am lost, I checked the R mailing help, but could not find anything directly. I loaded package Brobdingna...

...0.5008249 0.4960220 [50] 0.4960220 > #print(eW$vector) > > > #library(MASS) > #n = nrow(C) > #y = runif(n) > > #lm1 = lm.gls(y~1, W = C) > #summary(lm1) > > eW <-eigen(C, symmetric = TRUE) > l_eW <- eW$values > # the thirdteens eigenvector contrains NaN > print(eW$values) [1] 4.5730646 4.5730646 3.3066738 3.3066738 3.3066738 3.3066738 2.3935268 [8] 2.3935268 1.9367508 1.9367508 1.9347787 1.9347787 1.4276845 1.4276845 [15] 1.4276845 1.4276845 0.9858318 0.9858318 0.9858318 0.9858318 0.9123115 [22] 0.9123115 0.7945283 0.7945283 0.7880493 0.78804...

...uage. I was perfoming some calculation on fractional exponential and I found a strange behaviour. I do not know if it is really a bug, but I would expect a different answer from R. I was trying the following : x <- seq(-3,3, by =0.1) n <- 2.2 y <- exp(-x^n) well, the y vector contains (NaN for all negative value of x) but if you ask for single value calculation like y <- exp(-(-3)^2.2) or y <- exp(-(-2.9)^2.2) the answer is correct. It seem it does not make the calculation in vector form. I got the same behaviour (NaN) in a for loop > for(i in 1:length(x)) y[i]=exp(...

Dear R What does NaN mean? I recently did a correlation on a batch of data for some reason it didn't like one column cor(sleep,use="complete.obs") BodyWt BrainWt SlowSleep ParaSleep TotalSleep BodyWt 1.00000000 0.95584875 -0.3936373 -0.07488845 -0.3428373 BrainWt 0.9558...

...read a file with observations taken every 30 minutes: > aur2009=read.csv(paste(datadir,"AUR_ECPP_2009.csv",sep="/"),sep=";",stringsAsFactors=F) > aur2009[1:3,1:5] Date.Time E_filled E_filled_flag LE_filled LE_filled_flag 1 1/1/2009 0:00 0 NaN 5.86 NaN 2 1/1/2009 0:30 0 NaN 5.05 NaN 3 1/1/2009 1:00 0 NaN 5.56 NaN > delme = strptime(aur2009[,1], "%m/%d/%Y %H:%M") > aur2009[,1]=as.POSIXct(delme) Date.Time E_filled E_filled_flag LE...

Hi Oleg, On 01/09/14 18:46, Oleg Ranevskyy wrote: > Hi Duncan, > > I looked through the IEEE standard and here is what I found: > > *6.2 Operations with NaNs* > /"For an operation with quiet NaN inputs, other than maximum and minimum > operations, if a floating-point result is to be delivered the result shall be a > quiet NaN which should be one of the input NaNs"/. > > *6.2.3 NaN propagation* > /"An operation that prop...

...#39;t know whether this is a bug, or a problem with the way I >>> built R 2.0.1 (under Mac OS 10.3 on a G5), or something else. Can >>> anyone else confirm (or otherwise) that this happens in their R >>> 2.0.1 on Mac OS X? >>> >>> > d<-matrix(NaN,3,3) >>>> d >>> [,1] [,2] [,3] >>> [1,] NaN NaN NaN >>> [2,] NaN NaN NaN >>> [3,] NaN NaN NaN >>>> solve(d) >>> Error in solve.default(d) : Lapack routine dgesv: system is exactly >>> singular >>&g...

...;>#print(eW$vector) >> >> >>#library(MASS) >>#n = nrow(C) >>#y = runif(n) >> >>#lm1 = lm.gls(y~1, W = C) >>#summary(lm1) >> >>eW <-eigen(C, symmetric = TRUE) >>l_eW <- eW$values >># the thirdteens eigenvector contrains NaN >>print(eW$values) >> >> > [1] 4.5730646 4.5730646 3.3066738 3.3066738 3.3066738 3.3066738 2.3935268 > [8] 2.3935268 1.9367508 1.9367508 1.9347787 1.9347787 1.4276845 1.4276845 >[15] 1.4276845 1.4276845 0.9858318 0.9858318 0.9858318 0.9858318 0.9123115 >[22] 0.9123...

...FALSE) Autocorrelations of series 'X', by lag -4 -3 -2 -1 0 1 2 3 4 0.333 -0.500 0.667 -0.833 1.000 -0.833 0.667 -0.500 0.333 > ccf(x, y, plot = FALSE) Autocorrelations of series 'X', by lag -4 -3 -2 -1 0 1 2 3 4 NaN NaN NaN NaN NaN NaN NaN NaN NaN > y <- c(1, 1, 1, 1, 1, 1) > ccf(x, y, plot = FALSE) Autocorrelations of series 'X', by lag -4 -3 -2 -1 0 1 2 3 4 NaN NaN NaN NaN NaN NaN NaN NaN NaN i don't see why the result from ccf() would be NaN if the elements of y are a...

Hello. I have for example 4 or more dataframe which like such this example: date value 2006-11 0.4577647 2006-12 NaN 2006-10 0.1577647 2006-11 0.3577647 2006-12 NaN 2007-01 NaN 2007-02 NaN 2007-03 0.2956429 2007-01 0.3677647 2007-02 NaN They have the same length. I need merge by date undefined number of dataframes. End result output will look like such example date name1...

...vision code diff --git a/libswfdec/swfdec_script.c b/libswfdec/swfdec_script.c index d96ec0f..a9f16e0 100644 --- a/libswfdec/swfdec_script.c +++ b/libswfdec/swfdec_script.c @@ -283,6 +283,10 @@ swfdec_value_to_number_7 (JSContext *cx, if (JSVAL_IS_OBJECT (val)) { *d = *cx->runtime->jsNaN; return JS_TRUE; + } else if (JSVAL_IS_STRING (val) && + JS_GetStringLength (JSVAL_TO_STRING (val)) == 0) { + *d = *cx->runtime->jsNaN; + return JS_TRUE; } else { return JS_ValueToNumber (cx, val, d); } @@ -973,20 +977,17 @@ swfdec_action_binary (JSContext...

As far as I know, LLVM does not try very hard to guarantee constant folded NaN payloads that match exactly what the target would generate. —Owen > On Sep 16, 2014, at 10:30 AM, Oleg Ranevskyy <llvm.mail.list at gmail.com> wrote: > > Hi Duncan, > > I reread everything we've discussed so far and would like to pay closer attention to the the ARM'...

Hi, Thank you for all your helpful comments. To sum up, below is the list of correct folding examples for fadd: (1) fadd %x, -0.0 -> %x (2) fadd undef, undef -> undef (3) fadd %x, undef -> NaN (undef is a NaN which is propagated) Looking through the code I found the "NoNaNs" flag accessed through an instance of the FastMathFlags class. (2) and (3) should probably depend on it. If the flag is set, (2) and (3) cannot be folded as there are no NaNs and we are not guaranteed t...

This is some interesting: > -8^(1/3) [1] -2 > x=(-8:8) > y=x^(1/3) > y [1] NaN NaN NaN NaN NaN NaN NaN NaN 0.000000 1.000000 [11] 1.259921 1.442250 1.587401 1.709976 1.817121 1.912931 2.000000 So, can anybody explain this?! (Why is x[1]^(1/3)=y[1]=NaN, but -8^(1/3)=-2?) Thx!!! [[alternative HTML version deleted]]

...Thank you for all your helpful comments. >> >> To sum up, below is the list of correct folding examples for fadd: >> (1) fadd %x, -0.0 -> %x >> (2) fadd undef, undef -> undef >> (3) fadd %x, undef -> NaN (undef is a NaN which >> is propagated) >> >> Looking through the code I found the "NoNaNs" flag accessed through >> an instance >> of the FastMathFlags class. >> (2) and (3) should probably depend on it. >> If the flag is set, (2) and (3) cann...

Hi Oleg, On 01/09/14 15:42, Oleg Ranevskyy wrote: > Hi, > > Thank you for your comment, Owen. > My LLVM expertise is certainly not enough to make such decisions yet. > Duncan, do you have any comments on this or do you know anyone else who can > decide about preserving NaN payloads? my take is that the first thing to do is to see what the IEEE standard says about NaNs. Consider for example "fadd x, -0.0". Does the standard specify the exact NaN bit pattern produced as output when a particular NaN x is input? Or does it just say that the output is a N...