search for: approx

Hello all, I'm struggling to grasp how I might use lapply() instead of looping to run approx() on a list consisting of multiple arrays - each of equal dimension. But simpler than that, I haven't been able to successfully apply approx() to an array, unless I loop through the third dimension and extract the matrix. See example code below... Any suggestions will be gratefully received. Th...

Hi, I am trying to construct a logistic regression model from my data (104 patients and 25 events). I build a full model consisting of five predictors with the use of penalization by rms package (lrm, pentrace etc) because of events per variable issue. Then, I tried to approximate the full model by step-down technique predicting L from all of the componet variables using ordinary least squares (ols in rms package) as the followings. I would like to know whether I am doing right or not. > library(rms) > plogit <- predict(full.model) > full.ols <- ols(plog...

Hi, I have a object 'zoo': dim(zz) [1] 720 5551 where some columns only have NA's values (representing land data in a sea surface temperature dataset) I find straightforward the use of 'na.approx' for individual columns from the zz matrix, but when applied to the whole matrix: zz.approx<-na.approx(zz) Erro en approx(along[!na], y[!na], along[na], ...) : need at least two non-NA values to interpolate The message is clear, but how do I could skip those 'full-NA's'...

Full_Name: joerg van den hoff Version: 1.9.1 OS: MacOS and SunOS Submission from: (NULL) (149.220.4.88) follow up to ID 7166. something like approx(c(1,2),c(NA,NA),1.5,rule=2) crashes 1.9.1 on both systems (MacOS 10.3.5.: bus error, SunOS 5.9: segmentation fault) even if xout is within given x range (as in example above) where rule=2 seems not be relevant anyway.

I have now committed a version "fulfilling" your wish, partly at least, to R-devel . In the new approx(*, na.rm=FALSE) cases, the result of how NA's are treated does depend on the 4 different extrapolation rules {1, 2, 1:2, 2:1} The main reason was that I kept the low level code in C to do +- what it did before which automatically was using 'rule' to determine these results. The he...

Readers, Am trying to use the function 'approx' to interpolate time series data sets: data1: 01:23:40 5 01:23:45 10 01:23:50 12 01:23:55 7 data2: 01:23:42 01:23:47 01:23:51 01:23:54 The objective is to obtain interpolated values of 'data1' column 2 (5, 10, 12, 7) for the times shown in data2. Tried the following command but recei...

Dear R-devel, The following line crashes R > approx(1, 1, 0, method='const', rule=2, f=0, yleft=NULL, ties='ordered')$y Process R:2 exited abnormally with code 5 at Tue Jul 21 14:18:09 2009 > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status ma...

There is something I do not think is right in the approx() function in base R, with method="constant" and in the presence of NA values. I have 3.6.0, but the behavior seems to be the same in earlier versions. My suggested fix is to add an "na.rm" argument to approx(), as in mean(). If this argument is FALSE, then NA values should be...

There is something I do not think is right in the approx() function in base R, with method="constant" and in the presence of NA values. I have 3.6.0, but the behavior seems to be the same in earlier versions. My suggested fix is to add an "na.rm" argument to approx(), as in mean(). If this argument is FALSE, then NA values should be...

Hi, I am looking for a function that generalizes 'approx' to two (or more) dimensions. The references on the approx help page point toward splines, but a) splines is what I am trying to avoid in the first place and b) splines (except for mgcv splines) seem to be one dimensional. Here is a more detailed account. Using mgcv:gam I fit an additive mode...

>>>>> Robert Almgren >>>>> on Fri, 3 May 2019 15:45:44 -0400 writes [ __ to R-help __ -- here diverted to R-devel on purpose] > There is something I do not think is right in the approx() > function in base R, with method="constant" and in the > presence of NA values. I have 3.6.0, but the behavior > seems to be the same in earlier versions. (of course; the behavior has been unchanged, and "as documented" forever) > My suggested fi...

...79.57 12 C 2012-01-01 13:00:21.847 81.53 13 C 2012-01-01 13:00:21.860 81.53 14 C 2012-01-01 13:00:21.873 NA 15 C 2012-01-01 13:00:43.493 84.69 16 D 2012-01-01 12:01:21.520 24.63 17 D 2012-01-01 12:02:18.880 21.13 I tried the following using na.approx from zoo package interpolatedPrice<-unlist(tapply(portfolio$Price, portfolio$St, na.approx, portfolio$timestamp, na.rm=FALSE)) but keep getting error "Error in na.approx.default(X[[1L]], ...) : x and index must have the same length" I checked the length of every variable in the fo...

Hello, I am looking for a highly efficient matrix version of linear interpolation (like approx). As an example I have data like follows: x<-data.frame(time=1:20, x=(1:20)/10, y=runif(20)) t <- seq(1.5, 15.5 ,by=0.5) # and I found the following solution: nam <- names(x) app <- lapply(x[2:3],approx,x=x$time, xout=t) r <- c(1, 2*(1:(ncol(x)-1))) x.new <- as.data.frame(unlis...

Full_Name: joerg van den hoff Version: 1.9.0 and 1.7.1 OS: MacOS (1.9.0), SunOS (1.7.1) Submission from: (NULL) (149.220.4.88) something like (sure not the originally intended input, but something like this can happen...): approx(c(1,2),c(NA,NA),1.5) crashes R (bus error under MacOS, segmentation fault under SunOS). search of the bug archive did not work. I hope this bug was not reported already regards joerg

...PSTAMP: Steve Christensen INSTDATE: Mar 26 2002 18:15 EMAIL: steve at smc.vnet.net STATUS: completely installed FILES: 25 installed pathnames 2 shared pathnames 4 directories 4 executables 247 blocks used (approx) PKGINST: SMCossh NAME: openssh CATEGORY: application ARCH: sparc VERSION: 3.1p1 BASEDIR: /usr/local VENDOR: The OpenSSH Group PSTAMP: Steve Christensen INSTDATE: Mar 26 2002 18:30 EMAIL: steve at smc.vnet.net STATUS: completely installed F...

...redict(lm1, newdata = data.frame(volume = new_volume)) pred_elev The results just estimated for the 159 rows of the dataframe, not the new volume. > tail(pred_elev) 154 155 156 157 158 159 2254.296 2254.296 2254.296 2254.296 2254.296 2254.296 Also I have used the approx function, but it does not work for the new volume, too. > a = x6[,1] > b = x6[,4] > estimate <- 3,000,000,000 > appro <- approx(b,a, xout = estimate) > appro $x [1] 3e+09 $y [1] NA I do not know why it has happened. Is there any way to do this? Or maybe there is another way...

...1277648894 0.304075 -0.793544 263 25.237911 103 5 1277648898 0.304085 -0.793455 263 33.322525 108 6 1277648902 0.304064 -0.793387 256 40.042988 115 As you can see, the samples have irregular holes in the time column. How can I fill in the missing samples using na.approx? I've tried to creating a blank series with no gaps and combine them, but "merge" just adds columns and "rbind" compains about duplicate indexes. P.S. My GPS still has holes in the data when I turn off "smart recording" :( Thanks, Jason

This is R1.9.0, windows XP I run across the following: It is admittedly silly, but it shouldn't segfault: test <- list(y=1:100, x=rep(NaN, 100)) approx(test, xout=rep(NaN, 10)) Bombs! Kjetil Halvorsen

Hi all, I tried to produce some kriged surfaces with geoR (latest version). The size of the grid should be around 900 x 650 cells (what I find is not a very big grid), and the number of points is around 2500. The command krige.conv stopped after arround 5 min saying it can not allocate a vector with around 1.5 billion units. Sounds reasonable. Is there a workaround? How would I partition the

...d is a file with %%DocumentMedia: special 216 216 0 () () %%BoundingBox: 0 0 216 216 (where of course "216" means 3*72 points = 3 inches), and when this is displayed (using 'gv') I see the pie-chart distinctly offset from the centre of the BoundingBox, with its leftmost point at approx X=84 points from the left of the box (at 0), its rightmost point at X=162 approx 54 points from the right of the box (at 216), so the diameter of the circle is approx 162-84 = 78 points (a bit over 1 inch), and the centre of the circle is approx at 123 points (X) which is 15 points to the right of...