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.
>>
>
[[alternative HTML version deleted]]
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. >>> >>[[alternative HTML version deleted]]
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. >>>> >>>[[alternative HTML version deleted]]
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 <mzp3769 at gmail.com>
Date: Wed, Dec 26, 2018 at 6:32 PM
Subject: Re: [R] fields package question
To: Eric Berger <ericjberger at gmail.com>
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 <ericjberger at gmail.com>
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 <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.
[[alternative HTML version deleted]]