r - R で 3D 曲面プロットを作成するにはどうすればよいですか?

Question

多重線形回帰を使用しています。

ax*y 相互作用の結果である z 出力があります。

これらのデータの表面プロットを作成したいのですが、うまくいきませんでした。コマンドーのワイヤーフレームと persp を試しましたが、コマンドーをまだ設計できていないようです。

美しい scatterplot3d を作成しましたが、同じデータの表面プロットを作成するにはどうすればよいですか?

Answer 1

線形回帰で 3D サーフェスを当てはめたと言うとき、何をしているのかわかっていると思います。persp期待されるのは、 x 値と y 値のグリッドとwireframe、それらの各グリッド ポイントでの予測された z 高さです。を使用してこれを生成できますexpand.grid。これが指標となる例です。

preddf <- expand.grid(x=seq(xmin, xmax, len=51),
                      y=seq(ymin, ymax, len=51))

preddf$z <- predict(model, preddf)

persp(preddf)

、 、およびを予測子の範囲に置き換えxmin、51 をグリッドの目的のサイズ/密度に置き換えます。xmaxyminymax

Answer 2

別の代替手段expand.gridは、使用するouterことです。これはより高速です。

@HongOoi データの使用:

z <- outer(x,y, function(x,y) 
  predict(mod,data.frame(x=x,y=y)))
persp(x,y,z)