Tạo biểu đồ đường bao đầy bằng cách sử dụng dữ liệu trong danh sách

Question

Tôi có tập dữ liệu bao gồm 3 cột trong tệp .csv. Hai cột đầu tiên là tọa độ bản đồ và phần thứ ba là tỷ lệ phần trăm kẽm được tìm thấy trong lỗ khoan tại các tọa độ bản đồ tương ứng. Tôi muốn tạo một bản đồ đường đồng mức để thể hiện sự thay đổi nồng độ Zn với khoảng cách. Tất cả các ví dụ về mã tôi đã có thể tìm thấy dữ liệu sử dụng dưới dạng ma trận, trong khi tôi có trong danh sách. Tôi đã thử một số cách khác nhau để vẽ sơ đồ này mà tôi đã trình bày bên dưới, phần lớn các phương thức cung cấp cho tôi thông báo lỗi dọc theo dòng "đối tượng x không tìm thấy" mà tôi nghĩ là làm với bố cục dữ liệu của tôi. Có ai biết cách để làm điều này không? Tôi đã thêm một tập dữ liệu tương tự để khai thác bên dưới. Nhờ sự giúp đỡ trước. Holly

dữ liệu thiết lập:

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 

Code:

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

Answer 1

Làm việc với ggplot2, bạn có thể tạo một âm mưu đường viền với dữ liệu ví dụ của bạn thiết lập sử dụng:

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

để tặng

Bạn có một số vấn đề với mã số ggplot2. Đặc biệt, nơi bạn đặt thẩm mỹ bạn có:

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