2012-09-26 10 views
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 

コード:

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

に:ggplot2での作業

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

答えて

8

を、あなたはあなたの例のデータを使用して設定して等高線図を作成することができますgive

enter image description here

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

aes(x$x="X", y$y="Y", z$z="Zn") 
+0

ああ、私は間違った機能を使っていました。これは機能しますが、色で線の間を埋めることは可能ですか? –

+0

試してみてください: 'g + stat_density2d(aes(塗りつぶし= ..濃度..)、contour = FALSE、geom =" tile ")' – csgillespie

+0

ありがとう! 2つを組み合わせると、間違いなく私が探しているプロットのタイプです。しかし、密度の鍵は亜鉛の価値を与えるものではなく、e-08の値を与えます。 –

関連する問題