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

Question

Я хотел бы построить переменное количество точек по мере увеличения моего размера выборки. Однако по какой-то причине функция «вариограмма» только каждый раз отображает 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 очков. Кто-нибудь знает, почему? Я не видел это по умолчанию.

Answer 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', так как линия, соединяющая образцы вариограммы, предлагает больше, чем есть, действительно.