-1
我試圖在xy.plots上面顯示一些變種圖模型。 panel.plots除了我想添加到對應的子圖的行列表之外。用R格子包畫出線條列表
require(gstat)
require(sp)
data(meuse)
names(meuse)
#make directional variograms
b<-variogram(log(zinc)~1, meuse, alpha = c(0, 45, 90, 135))
#split the variogram data by direction
a<-lapply(1:length(unique(b$dir.hor)),
function(i) subset(b, grepl(unique(b$dir.hor)[[i]], b$dir.hor)))
#get the model fit parameters for each directional variogram
a<-lapply(1:length(unique(b$dir.hor)),
function(i) fit.variogram(a[[i]], vgm(0.5,"Exp", 1200, 0.5)))
#generate model data for the directional variograms
a<-lapply(1:length(a), function(i)
variogramLine(a[[i]], maxdist=1500))
#plot
require(lattice)
plot(b, ylim=c(0,1.2), xlim=c(0,1500), cex=1.5,
panel = function(x, y, ...) {panel.xyplot(x, y, ...)
panel.abline(v=1000, lwd=1, lty=3, col=2)
panel.lines(a[[i]], lwd=2, lty=2, col=4)
})
在劇情的最後一個參數不起作用,因爲它是一個列表(見陰謀顯示錯誤),不知道如何行的列表應該格繪製。幫助將不勝感激!
http://stackoverflow.com/help/mcve –
如@哈克-R所提到的,請提供一個工作示例。 'meuse'數據位於'library(sp)'FWIW中。 –
@RomanLuštrik完成 – ToNoY