[R] fields package question

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