0

網膜の密度細胞に異なる密度点を使用して地形図を作成しようとしています。の例volcano3dは本当に素晴らしく、それは私にアイデアを与えたので、私はこの例に従って私の質問を言い表すと思います。私の質問は、x、y座標がプロットされた島の等高線がある場合geom_path、パス内にのみ等高線マップを作成することは可能ですか？私はそれをやろうとしましたが、rは、問題は私のx軸とy軸が連続していないことであると教えてくれます。これは私のデータの最初の例です。xとyは座標であり、freqはz軸である必要があります。データ全体を入れていないのは残念ですが、等高線には3000ポイントを超えており、内部に約250のサンプリングがあるため、長いデータセットになります。

p <- ggplot(data = contour.xy, aes(x = x, y = y)) + 
  geom_path(data = contour.xy, aes(x = x, y = y)) + 
  geom_path(data = opticnerve.xy, aes(x = x, y = y)) + 
  geom_contour(data = counts, aes(x = "x", y = "y", z = "freq")) + 
  coord_equal()

画像を入れてみましたが、新規ユーザーなのでスパム対策規制のため写真を添付できません。

私は本当に助けやアドバイスに感謝します

みんなに感謝とメリークリスマス

4

2 に答える 2

1

最後に、いくつかの調査を行った後、私が探していたものを正確に実行するのに役立つパッケージを見つけました。これはエイドリアン・バデリーによって書かれたspatstatであり、彼はスクリプトの一部をとても親切に手伝ってくれます。基本的には、不規則な形状のウィンドウを作成することです。この場合は私の「島」であり、分析する必要のあるすべてのポイントがウィンドウの内側にあります。ウィンドウを作成するために、xベクトルをxpとして、yベクトルをypとして使用して次のコードを使用しました。map<-owin（poly = list（x = xp、y = yp））plot（map）

次に、すべてのポイントをポイントパターン（ppp）として作成できます

本当にあなたのすべての助けに感謝しますが、彼のサポートのためにエイドリアン・バデリーに特に感謝します

于 2011-12-22T09:25:46.170 に答える
0

これがvolcano3d例です（再現可能な問題を提供するため。コードはパッケージexample(direct.labels)から取得されます）。directlabels

volcano3d <- melt(volcano)
colnames(volcano3d) <- c("x", "y", "z")

次のように、データの等高線図（この場合はマウンガファウ）を描くことができます。

(p <- ggplot(volcano3d, aes(x, y, z = z)) + 
  geom_contour()
)

aes文字列ではなく式を受け取るサンプルコードにも問題があります。だからあなたはどちらかを持っている必要があります

aes(x = x, y = y, z = freq)

また

aes_string(x = "x", y = "y", z = "freq")

x通常、との引数名を明示的に指定する必要はないyので、実際aes(x, y, z = freq)にはそうします。

于 2011-12-21T11:22:36.730 に答える