ggplotで黄土の表面をプロットする方法

Question

私はこのコードを持っています。私は私のデータフレームの黄土色のサーフェスを作成します。

library(gstat) 
    library(sp) 
    x<-c(0,55,105,165,270,65,130,155,155,225,250,295, 
     30,100,110,135,160,190,230,300,30,70,105,170, 
     210,245,300,0,85,175,300,15,60,90,90,140,210, 
     260,270,295,5,55,55,90,100,140,190,255,285,270) 
    y<-c(305,310,305,310,310,260,255,265,285,280,250, 
     260,210,240,225,225,225,230,210,215,160,190, 
     190,175,160,160,170,120,135,115,110,85,90,90, 
     55,55,90,85,50,50,25,30,5,35,15,0,40,20,5,150) 
    z<-c(870,793,755,690,800,800,730,728,710,780,804, 
     855,813,762,765,740,765,760,790,820,855,812, 
     773,812,827,805,840,890,820,873,875,873,865, 
     841,862,908,855,850,882,910,940,915,890,880, 
     870,880,960,890,860,830) 

    dati<-data.frame(x,y,z) 

    x.range <- as.numeric(c(min(x), max(x))) 
    y.range <- as.numeric(c(min(y), max(y))) 

    meuse.loess <- loess(z ~ x * y, dati, degree=2, span = 0.25, 
         normalize=F) 
    meuse.mar <- list(x = seq(from = x.range[1], to = x.range[2], by = 1), y = seq(from = y.range[1], 
                        to = y.range[2], by = 1)) 
    meuse.lo <- predict(meuse.loess, newdata=expand.grid(meuse.mar), se=TRUE) 

は今、私はggplot2機能をmeuse.lo[[1]]をプロットしたい...しかし、私は、Y（グリッドの座標）とZ（補間値）の列のxでのデータフレームにmeuse.lo[[1]]を変換する方法を知りません。ありがとう。

Answer 1

あなたの問題と簡単な解決策はここにあるあなたがloess()のために新しいデータを生成するためにgrid.expand()を使用する場合loess()が行列を返すことです。

これは?loess.predictのヘルプに記載されている：

NEWDATAは（要求された場合とSEの）、予測をexpand.gridするための呼び出しの結果であった場合には、適切な寸法の配列となります。

ここでもgrid.expand()を使用して新しいデータを計算できますが、この関数は強制的にデータフレームを返して属性を削除します。 ?grid.expandから

：

KEEP.OUT.ATTRS： "out.attrs" 属性を示す論理（下記参照）が計算されて返さなければなりません。

ので、これを試してみてください。

nd <- expand.grid(meuse.mar, KEEP.OUT.ATTRS = FALSE) 
meuse.lo <- predict(meuse.loess, newdata=nd, se=TRUE) 

# Add the fitted data to the `nd` object 
nd$z <- meuse.lo$fit 

library(ggplot2) 
ggplot(nd, aes(x, y, col = z)) + 
    geom_tile() + 
    coord_fixed() 

結果：

Answer 2

ggplot2は、おそらく3次元グラフには最適ではありません。しかし、ここでrgl

library(rgl) 
plot3d(x, y, z, type="s", size=0.75, lit=FALSE,col="red") 
surface3d(meuse.mar[[1]], meuse.mar[[2]], meuse.lo[[1]], 
     alpha=0.4, front="lines", back="lines")