Actually, let's set it grid_new.l <- list(abcissa=c(-15.0,-14.5),ordinate=y) to avoid out of bounds
On Tue, Dec 25, 2018 at 4:41 PM M P <mzp3...@gmail.com> wrote: > 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.
