search for: 1e

...I tried the following codes given below. But it produced the error messages.Can any one tell me the code in 'drc' package to find out the EC (Effective Concentration) values and Confidence Intervals. The data is given below (also attached a ".csv" file). dat conc Length 1 1e-04 28.55 2 1e-04 14.47 3 1e-04 29.93 4 1e-04 31.66 5 1e-04 10.78 6 1e-04 16.83 7 1e-04 10.52 8 1e-04 31.44 9 1e-04 16.56 10 1e-04 10.30 11 1e-04 20.08 12 1e-04 14.04 13 1e-04 19.93 14 1e-04 14.94 15 1e-04 38.68 16 1e-04 37.92 17 1e-04 3.54 18 1e-04 13.58...

I am getting theses values when I disply my ncdf file 1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 [227,] -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 [228,] -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 [229,] -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 -1e+34 [230,] -1e+34 -1...

...e in order that we collaboratively propose an implementation suiting all needs. Thanks, Karl Forner Annex: I ran the benchmarks on a linux Intel(R) Xeon(R) with 2 cpus of 4 cores each ( CPU E5520 @ 2.27GHz). type threads n error time time_per_chunk 1 lecuyer 1 1e+07 2.105472e-04 1.538 0.00153800 2 lecuyer 1 1e+08 4.441492e-05 15.265 0.00152650 3 lecuyer 1 1e+09 2.026819e-05 153.209 0.00153209 4 lecuyer 2 1e+07 3.182633e-04 0.821 0.00082100 5 lecuyer 2 1e+08 7.375036e-05 7.751 0.00077510 6...

...ead.cat <- cut(df$ead, breaks=c(0, 100000, 500000, 1000000, 2000000, 5000000 , 10000000, 100000000) ) df          df_sorted <- df[order(df$rating),]      # the output is as given below. > df_sorted    rating         ead                     ead.cat 1       A          169229.93        (1e+05,5e+05] 3       A         5877794.25        (5e+06,1e+07] 6       A            21000.00               (0,1e+05] 12      A          150000.00       (1e+05,5e+05] 13      A          167490.00       (1e+05,5e+05] 18      A             9100.00               (0,1e+05] 23      A             3000.00   ...

Hi, This problems has bothered me for the lase couple of hours. > 1e-100==0 [1] FALSE > (1-1e-100)==1 [1] TRUE How can I tell R that 1-1e-100 does not equal to 1, actually, I found out that > (1-1e-16)==1 [1] FALSE > (1-1e-17)==1 [1] TRUE The reason I care about this is that I was try to use qnorm() in my code, for example, > qnorm(1e-100) [1] -21.2...

dear all, here's a couple of questions that puzzled me in these last hours: ##### issue 1 qnorm(1-10e-100)!=qnorm(10e-100) qnorm(1-1e-10) == -qnorm(1e-10) # turns on to be FALSE. Ok I'm not a computer scientist but, # but I had a look at the R inferno so I write: all.equal(qnorm(1-1e-10) , -qnorm(1e-10)) # which turns TRUE, as one would expect, but all.equal(qnorm(1-1e-100) , -qnorm(1e-100)) # turns FALSE: Indeed # qnor...

Dear mailing list readers! Using R and the plot function I stumbled over this little issue: plot(x=NULL, y=NULL, xlim=range(1:10), ylim=range(1e-9:1), log="y") produces an plot empty plot, where the y-axis is in the specified range. Changing ylim to ylim=range(1e-9:1e-3) creates an y-axis in the range of 1e-12 to 1e-6. This appearance of the y-axis is a bit unexpected for me. Could anyone point me to a FAQ entry why t...

Full_Name: Version: 2.7.1 (2008-06-23) OS: windows vista Submission from: (NULL) (216.82.144.137) Hello, I have been working with various probability distributions in R, and it seems the gamma distribution is inaccurate for some inputs. For example, qgamma(1e-100, 5e-101, lower.tail=FALSE) gives: 1.0. However, it seems this is incorrect; I think the correct answer should be 0.082372029620717283. When I check these numbers using pgamma, I get: pgamma(1,5e-101, lower.tail=FALSE) = 9.1969860292859463e-102 and pgamma(0.082372029620717283,5e-101, lower.tail...

...t;------------------------" [1] "file /users/huxxx249/MEGANv2.04/Output/latlong_spec2_WOCONV_150s_2004187.ncf has 31 variables:" [1] "int TFLAG[DATE-TIME,VAR,TSTEP] Longname:TFLAG Missval:NA" [1] "float ISOPRENE[COL,ROW,LAY,TSTEP] Longname:ISOPRENE Missval:1e+30" [1] "float TRP1[COL,ROW,LAY,TSTEP] Longname:TRP1 Missval:1e+30" [1] "float MEOH[COL,ROW,LAY,TSTEP] Longname:MEOH Missval:1e+30" [1] "float ACET[COL,ROW,LAY,TSTEP] Longname:ACET Missval:1e+30" [1] "float CH4[COL,ROW,L...

Rsync 2.5.5; transfer via ssh; sparc-sun-solaris2.5 (Ultra 1); gcc 2.8.1. Every time rsync on this box causes boundary violation. I attached two strange examples. I have a few rsync's corefiles (100MB and above) and I can examine them. I will be grateful for any help. -Mirek <-------------------------------------------------------------------------> received 966761 names done

...mes take 1000x longer than taking the same number of samples from a slightly larger population. For my application, this resulted in a runtime of several hours rather than a few seconds. Looking into it, I saw that sample.int() is hardcoded to switch algorithms when the population is larger than 1e+7, and I was straddling this boundary: sample.int <- function(n, size = n, replace = FALSE, prob = NULL) { if (!replace && is.null(prob) && n > 1e7 && size <= n/2) .Internal(sample2(n, size)) else .Internal(sample(n, size, replace, prob)) } do_sam...

...ge.mode(seq(from=1L,to=2.5)) => integer --------------------------------- On the other hand, according to ? seq Currently, the default method returns a result of type '"integer"' if 'from' is (numerically equal to an) integer This is not correct: > seq(1e16,1e16+10) [1] 1e+16 1e+16 1e+16 1e+16 1e+16 1e+16 1e+16 1e+16 1e+16 1e+16 1e+16 > seq(1e16,1e16+10)-1e16 [1] 0 0 2 4 4 4 6 8 8 8 10 # not enough precision to represent all the distinct integers, OK > storage.mode(seq(1e16,1e16+10)) [1] "double" In this case, the ? :...

Hello All I am looking for high precision values for the normal distribution in the tail,(1e-10 and 1 - 1e-10) as the R package that I am using sets any number which is out of this range to these values and then calls the qnorm and qt function. What I have noticed is that the qnorm implementation in R is not symmetric when looking at the tails. This is quite surprising to me, as it is wel...

Full_Name: Morten Welinder Version: 1.6.1 OS: Solaris/sparc Submission from: (NULL) (65.213.85.144) qbeta(0.1, 1e-8, 0.5, TRUE, FALSE) seems to hang for me.

...ear R-helpers, I need some help to produce a set of contour plots that I am trying to make in order to compare surfaces between the levels of a factor. For example: library(lattice) g <- expand.grid(x = 60:100, y = 1:25, ti = c("a","b","c")) g$z <-with(g, (-1e-4*x-1e-3*y-1e-5*x*y)*(ti=="a") + (1e-2*x-1e-3*y-1e-4*x*y)*(ti=="b") + (1e-3*x-1e-3*y-1.e-5*x*y)*(ti=="c") ) contourplot(z~ x * y|ti, data = g, cuts=20, pretty=T, screen = list(z = 30, x = -60)) As you can see in...

Full_Name: Mike Meredith Version: 2.3.0 OS: WinXP SP2 Submission from: (NULL) (210.195.228.29) Using dt() with a non-centrality parameter and near-zero values for 'x' results in erratic output. Try this: tst <- c(1e-12, 1e-13, 1e-14, 1e-15, 1e-16, 1e-17, 0) dt(tst,16,1) I get: 0.2381019 0.2385462 0.2296557 0.1851817 0.6288373 3.8163916 (!!) 0.2382217 The 0.238 values are okay, the others nonsense, and they cause confusing spikes on plots of dt() vs 'x' if 'x' happens to include tiny values....

Hi, I am trying to extract subset of data from my original data frame based on some condition. For example : (mydf -original data frame, submydf - subset dada frame) >submydf = subset(mydf, a > 1 & b <= a), here column a contains values ranging from 0.01 to 100000. I want to extract only those matching condition 1 i.e a > . But when i execute this command it is

...T[ind.r], Rat[ind.r], sdrat[ind.r], log = TRUE) dB.r <- dnorm(y[, wholeranged], Dat[, wholeranged], sddat[, wholeranged], log = TRUE) } else { dA.r <- 0 dB.r <- 0 } dA.p[dA.p < -1e+08] <- -1e+08 dA.p[dA.p > +1e+08] <- +1e+08 dA.r[dA.r < -1e+08] <- -1e+08 dA.r[dA.r > +1e+08] <- +1e+08 drat <- sum(dA.p, dA.r, na.rm = TRUE) if (unif) drat <- 0 if (any(RAT < minRat | RAT > maxRat, na.rm = T...

...ice is a vector of 789 asset prices of the following form, > mode(price) [1] "numeric" > price[1:10] [1] 1.1315 1.1314 1.1313 1.1314 1.1315 1.1316 1.1315 1.1317 [9] 1.1316 1.1316 Then the returns are defined to be the differences, > ret <- diff(price) > ret[1:10] [1] -1e-04 -1e-04 1e-04 1e-04 1e-04 -1e-04 2e-04 -1e-04 [9] 0e+00 0e+00 Here R recognizes the first two entries and also the second ones to be equal (I did not found any two equal entries, where R didnt recognized this relation), > ret[1]==ret[2] [1] TRUE > ret[3]==ret[4] [1] TRUE But after...

...7-06-22 19:49 GMT+02:00 Uwe Ligges <ligges at statistik.tu-dortmund.de>: > On 22.06.2017 17:11, Bernd Funovits wrote: >> >> Hello, >> >> I experienced some unexpected behaviour while determining the rank of matrices (sometimes 1x1 matrices): >> base::qr(matrix(1e-20))$rank returns 1 (incorrect) >> base::qr(diag(c(1, 1e-20)))$rank returns 2 (incorrect) > > > Why do you believe this is incorrect? > > 0 != 1e-20 > > and 1e-20 is well representable without significant rounding errors given > > > .Machine$double.xmin > [1]...