fields package question
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 <mzp3769 at 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 <mzp3769 at 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 <ericjberger at 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 <mzp3769 at 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 at 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.