Thanks, Eric, for looking into that. The values are below and since I subset the new abcissa is smaller range grid_new.l <- list(abcissa=c(-15.0,-14.),ordinate=y) I am emailing form gmail - don't know why is using html to format when all is in ascii
x [1] -15.20180 -15.01948 -14.86533 -14.73180 -14.61402 -14.50866 -14.41335 [8] -14.32634 -14.24629 -14.17219 y [1] 0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35 0.40 0.45 z [5,] 0.6467642 0.6467642 0.6467642 0.6467642 0.6467642 0.6467642 0.6467642 [6,] 0.5597143 0.5597143 0.5597143 0.5597143 0.5597143 0.5597143 0.5597143 [7,] 0.4854133 0.4854133 0.4854133 0.4854133 0.4854133 0.4854133 0.4854133 [8,] 0.4278326 0.4278326 0.4278326 0.4278326 0.4278326 0.4278326 0.4278326 [9,] 0.3834149 0.3834149 0.3834149 0.3834149 0.3834149 0.3834149 0.3834149 [10,] 0.3433031 0.3433031 0.3433031 0.3433031 0.3433031 0.3433031 0.3433031 [,8] [,9] [,10] [1,] 1.1900951 1.1900951 1.1900951 [2,] 1.0636935 1.0636935 1.0636935 [3,] 0.8927228 0.8927228 0.8927228 [4,] 0.7554456 0.7554456 0.7554456 [5,] 0.6467642 0.6467642 0.6467642 [6,] 0.5597143 0.5597143 0.5597143 [7,] 0.4854133 0.4854133 0.4854133 [8,] 0.4278326 0.4278326 0.4278326 [9,] 0.3834149 0.3834149 0.3834149 [10,] 0.3433031 0.3433031 0.3433031 On Tue, Dec 25, 2018 at 12:45 AM Eric Berger <ericjber...@gmail.com> wrote: > Since you don't provide lambda, rh or qext it is impossible to reproduce > what you are seeing. > Also note that in this mailing list HTML formatted emails are not passed > along. > > > > On Tue, Dec 25, 2018 at 4:13 AM M P <mzp3...@gmail.com> wrote: > >> Hello, >> I used commands below to obtain a surface, can plot it and all looks as >> expected. >> How do I evaluate values at new point. I tried as below but that produces >> errors. >> Thanks for suggestions/help. >> >> x <- log(lambda) >> y <- rh >> z <- qext[,,2] >> >> grid.l <- list(abcissa=x,ordinate=y) >> xg <- make.surface.grid(grid.l) >> out.p <- as.surface(xg,z) >> plot.surface(out.p,type="p") >> >> tried: >> grid_new.l <- list(abcissa=c(-15.0,-10.),ordinate=y) >> xg_new <- make.surface.grid(grid_new.l) >> >> out_new.p <- predict.surface(out.p,xg_new) >> >> results in this prompt: >> predict.surface is now the function predictSurface> >> >> [[alternative HTML version deleted]] >> >> ______________________________________________ >> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see >> 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. >> > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see 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.