In your printing out of z your copy-paste effort seems to have missed columns 1-7 of rows 1-4. It would be better if you would provide the data by using the dput() function, as in: dput(y) dput(z)
then copy-paste the output from that. On Wed, Dec 26, 2018 at 1:50 AM M P <mzp3...@gmail.com> wrote: > 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.
