Почему вариограмма всегда отображает 15 точек в R?

Я хотел бы построить переменное количество точек по мере увеличения размера моей выборки. Однако по какой-то причине функция "вариограмма" каждый раз отображает только 15 точек.

Я проверил, чтобы убедиться, что размер данных, которые я передаю, "вариограмма" правильно менялась - так и было.

library(gstat)
library(RandomFields)
library(lattice)
library(latticeExtra)
mod <- RMexp(var=1, scale=5) + RMtrend(mean=3)
# theoretical mean 3 

# (x,y) coordinates for simulation grid
x <- seq(0,50,by=0.5)
y <- seq(0,0,by=0.5)
xx <- rep(x, times=length(y))
yy <- rep(y, each=length(x))

zz <- RFsimulate(mod, x=xx, y=yy,spConform=FALSE)

field <- data.frame(x=xx,y=yy,z=zz)


d <- sample(zz,10)
g <- gstat(formula=z~1, locations=~x+y, data=raw.dat)

# N=10:
n10 <- sample(1:length(field[[1]]),10,replace=F)

#g <- gstat(formula=z~1, locations=~x+y, data=raw.dat) 
f10 = field[n10,]
g10 <- gstat(formula=z~1, locations=~x+y, data=f10)


raw.vgm <- variogram(g10) # create method of class "gstatVariogram"
plot(raw.vgm,main='Variogram of Raw Data for N = 10',type='b') # plot method     for class "gstatVariogram"

# N=25:
n25 <- sample(1:length(field[[1]]),25,replace=F)

#g <- gstat(formula=z~1, locations=~x+y, data=raw.dat) 
f25 = field[n25,]
g25 <- gstat(formula=z~1, locations=~x+y, data=f25)

# f25 имеет длину 25 - я проверил #

raw.vgm <- variogram(g25) # create method of class "gstatVariogram"
plot(raw.vgm,main='Variogram of Raw Data for N = 25',type='b') # plot method for class "gstatVariogram"

Обе сырые вариограммы показывают только 15 точек. Кто-нибудь знает почему? Я не видел это по умолчанию.

1 ответ

Решение

?variogram дает по умолчанию для аргумента width Значение cutoff/15, что вызывает по умолчанию 15 баллов. Если вы делаете значение для width меньше, вы увидите больше очков. Пытаться

raw.vgm <- variogram(g25, width = .5)
plot(raw.vgm,main='Variogram of Raw Data for N = 25')

чтобы получить больше очков, попробуйте изменить cutoff если вы чувствуете себя предприимчивым Я бы не рекомендовал type='b' поскольку линия, соединяющая точки выборки вариограммы, на самом деле предполагает больше, чем есть.

Другие вопросы по тегам