fields package question

classic Classic list List threaded Threaded
6 messages Options
mzp
Reply | Threaded
Open this post in threaded view
|

fields package question

mzp
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]]

______________________________________________
[hidden email] 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.
Reply | Threaded
Open this post in threaded view
|

Re: fields package question

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 <[hidden email]> 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]]
>
> ______________________________________________
> [hidden email] 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]]

______________________________________________
[hidden email] 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.
mzp
Reply | Threaded
Open this post in threaded view
|

Re: fields package question

mzp
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 <[hidden email]> 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 <[hidden email]> 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]]
>>
>> ______________________________________________
>> [hidden email] 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]]

______________________________________________
[hidden email] 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.
mzp
Reply | Threaded
Open this post in threaded view
|

Re: fields package question

mzp
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 <[hidden email]> 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 <[hidden email]>
> 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 <[hidden email]> 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]]
>>>
>>> ______________________________________________
>>> [hidden email] 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]]

______________________________________________
[hidden email] 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.
Reply | Threaded
Open this post in threaded view
|

Re: 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 <[hidden email]> 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 <[hidden email]> 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 <[hidden email]>
>> 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 <[hidden email]> 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]]
>>>>
>>>> ______________________________________________
>>>> [hidden email] 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]]

______________________________________________
[hidden email] 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.
Reply | Threaded
Open this post in threaded view
|

Fwd: fields package question

Eric Berger
Forwarding to the full R-help list.
The protocol is for everyone to see the full set of questions and answers.
This allows for:
    1. more people to contribute, which generally improves the quality of
the answer
    2. people with similar questions can find this thread via search and
see the full discussion

---------- Forwarded message ---------
From: M P <[hidden email]>
Date: Wed, Dec 26, 2018 at 6:32 PM
Subject: Re: [R] fields package question
To: Eric Berger <[hidden email]>


Thanks, dput outputs below. I'd like to evaluate surface e.g. at (-10.5,
0.935)
> dput(x)
structure(c(-10.9251387646973, -10.8977397823155, -10.8197783096724,
-10.7568035566262, -10.701995233594, -10.5966347583582, -10.4868838826637,
-10.4143132015642, -10.2601624697659, -10.1266311291125), .Dim = 10L)
> dput(y)
structure(c(0.899999976158142, 0.910000026226044, 0.920000016689301,
0.930000007152557, 0.939999997615814, 0.949999988079071, 0.959999978542328,
0.970000028610229, 0.980000019073486, 0.990000009536743), .Dim = 10L)
> dput(z)
structure(c(0.0170888844877481, 0.0163554549217224, 0.0143800554797053,
0.0132564017549157, 0.0123229259625077, 0.0108034228906035,
0.00943510234355927,
0.00868543982505798, 0.00763010373339057, 0.00705513963475823,
0.0173761900514364, 0.0166240241378546, 0.0145960496738553,
0.0134453792124987,
0.0124927284196019, 0.0109430542215705, 0.00954837258905172,
0.00878583826124668, 0.00772488862276077, 0.00716213695704937,
0.0178007110953331, 0.0170208644121885, 0.0149147948250175,
0.0137240244075656,
0.0127429272979498, 0.0111485198140144, 0.0097147086635232,
0.00893292296677828,
0.00786342192441225, 0.00731897540390491, 0.0180797930806875,
0.0172817632555962, 0.0151241403073072, 0.0139069128781557,
0.012907050549984,
0.0112831527367234, 0.0098235122859478, 0.00902892742305994,
0.00795362330973148, 0.00742132356390357, 0.0186295621097088,
0.0177957788109779, 0.015536243095994, 0.014266736805439,
0.0132297882810235,
0.0115476427599788, 0.0100369155406952, 0.00921681523323059,
0.00812961626797915, 0.00762135302647948, 0.0190353523939848,
0.0181752592325211, 0.0158403068780899, 0.0145321246236563,
0.0134677225723863,
0.0117424847558141, 0.0101939048618078, 0.00935473665595055,
0.00825834088027477, 0.00776780024170876, 0.0190353523939848,
0.0181752592325211, 0.0158403068780899, 0.0145321246236563,
0.0134677225723863,
0.0117424847558141, 0.0101939048618078, 0.00935473665595055,
0.00825834088027477, 0.00776780024170876, 0.0190353523939848,
0.0181752592325211, 0.0158403068780899, 0.0145321246236563,
0.0134677225723863,
0.0117424847558141, 0.0101939048618078, 0.00935473665595055,
0.00825834088027477, 0.00776780024170876, 0.0190353523939848,
0.0181752592325211, 0.0158403068780899, 0.0145321246236563,
0.0134677225723863,
0.0117424847558141, 0.0101939048618078, 0.00935473665595055,
0.00825834088027477, 0.00776780024170876, 0.0190353523939848,
0.0181752592325211, 0.0158403068780899, 0.0145321246236563,
0.0134677225723863,
0.0117424847558141, 0.0101939048618078, 0.00935473665595055,
0.00825834088027477, 0.00776780024170876), .Dim = c(10L, 10L))

On Wed, Dec 26, 2018 at 1:24 AM Eric Berger <[hidden email]> wrote:
>
> 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 <[hidden email]> 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 <[hidden email]> 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 <[hidden email]>
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 <[hidden email]> 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]]
>>>>>
>>>>> ______________________________________________
>>>>> [hidden email] 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]]

______________________________________________
[hidden email] 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.