5

.csv ファイルに 3 つの列で構成されるデータ セットがあります。最初の 2 列はマップ座標で、3 番目は、対応するマップ座標のボアホールで見つかった亜鉛のパーセンテージです。Zn 濃度が距離によって変化することを示す等高線図を作成したいと思います。私が見つけることができたコードの例はすべて、マトリックスの形式でデータを使用していますが、私のものはリストになっています。以下に示すように、これをプロットするいくつかの異なる方法を試しましたが、ほとんどの方法で、データのレイアウトに関係していると思われる「オブジェクト x が見つかりません」という行に沿ってエラー メッセージが表示されます。誰もこれを行う方法を知っていますか? 同様のデータセットを以下に追加しました。事前に助けてくれてありがとう。ホリー

データセット:

Statsrep <- structure(list(X = c(156000L, 156010L, 156421L, 156450L, 156500L, 156700L, 158420L, 158646L, 158970L, 159050L, 159050L, 159130L, 159155L), Y = c(143630, 143980, 147260, 145000, 146000, 142800, 146700, 145207, 147170, 145200, 144800, 147815, 145890), Zn = c(2, 8, 4, 0, 3, 0, 2, 7, 12, 0, 4, 19, 0)), .Names = c("X", "Y", "Zn"), row.names = c(1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L), class = "data.frame")
Statsrep

コード：

library(ggplot2)
Grade <- read.csv(file="filename.csv", header=TRUE, sep=",")
ggplot(Grade, aes(x$x="X", y$y="Y", z$z="Zn")) +
stat_contour()

library(lattice)
Grade <- read.csv(file="filename.csv", header=TRUE, sep=",")
levelplot(Grade ~x*y, data = Zn,
xlab = "Eastings", ylab = "Northings",
col.regions = terrain.colours)

Grade <- read.csv(file="filename.csv", header=TRUE, sep=",")
x$x <- X
y$y <- Y
z$z <- Zn
filled.contour(x$x, y$y, z$z, color = terrain.colours,
xlab = "Eastings", ylab = "Northings"),
plot.axes = {axis(1, seq(156000, 165000, by=1000)); axis(2, seq(142000, 150000, by=1000))},
key.title = title(main="Zn content\n(percent)"),
key.axes= axis(4, seq(0, 20, by = 2)))
4

1 に答える 1

8

ggplot2を使用すると、次のコマンドを使用して、サンプルデータセットを使用して等高線図を作成できます。

ggplot(Statsrep, aes(x=X, y=Y, z=Zn)) + 
    geom_density2d()

与えるために

ここに画像の説明を入力してください

ggplot2コードにいくつか問題がありました。特に、あなたが持っていた美学を設定する場所：

aes(x$x="X", y$y="Y", z$z="Zn")
于 2012-09-26T12:20:19.773 に答える