Re: [R] fields package question

2018-12-26 Thread Eric Berger
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 wrote: > Actually, let'

Re: [R] fields package question

2018-12-25 Thread M P
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 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,-

Re: [R] fields package question

2018-12-25 Thread M P
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

Re: [R] fields package question

2018-12-24 Thread Eric Berger
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 wrote: > Hello, > I used commands below to obtain a surface, can plot it and all loo

[R] fields package question

2018-12-24 Thread M P
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.g