[R] fields package question

[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 <mzp3769 using 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 using 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 using 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 using 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]]