Skip to content

fields package question

5 messages · Eric Berger, m p

Original
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.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>
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 <mzp3769 at gmail.com> wrote:

            

      
      
      
    

        

      
      
    

        

  
  


  


      

    

    

      
      

        


  

        

      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 -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:

        

            

      
      
      
    

        

      
      
    

            

      
      
      
    

        

  
  


  


      

    

    

      
      

        


  

        

      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 at gmail.com> wrote:

        

            

      
      
      
    

        

      
      
    

            

      
      
      
    

        

  
  


  


      

    

    

      
      

        


  

        

      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 at gmail.com> wrote: