CRAN Package Check Results for Package gstat

Last updated on 2025-12-14 12:50:00 CET.

Flavor	Version	T_install	T_check	T_total	Status
r-devel-linux-x86_64-debian-clang	2.1-4	20.15	185.80	205.95	OK
r-devel-linux-x86_64-debian-gcc	2.1-4	15.07	131.79	146.86	OK
r-devel-linux-x86_64-fedora-clang	2.1-4	34.00	241.07	275.07	ERROR
r-devel-linux-x86_64-fedora-gcc	2.1-4	37.00	225.02	262.02	ERROR
r-devel-windows-x86_64	2.1-4	37.00	243.00	280.00	OK
r-patched-linux-x86_64	2.1-4	21.91	176.80	198.71	OK
r-release-linux-x86_64	2.1-4	19.99	177.66	197.65	OK
r-release-macos-arm64	2.1-4				OK
r-release-macos-x86_64	2.1-4	15.00	205.00	220.00	OK
r-release-windows-x86_64	2.1-4	35.00	238.00	273.00	OK
r-oldrel-macos-arm64	2.1-4				OK
r-oldrel-macos-x86_64	2.1-4	15.00	189.00	204.00	OK
r-oldrel-windows-x86_64	2.1-4	45.00	289.00	334.00	OK

Check Details

Version: 2.1-4
Check: tests
Result: ERROR Running ‘allier.R’ Comparing ‘allier.Rout’ to ‘allier.Rout.save’ ... OK Running ‘blockkr.R’ Comparing ‘blockkr.Rout’ to ‘blockkr.Rout.save’ ... OK Running ‘covtable.R’ Comparing ‘covtable.Rout’ to ‘covtable.Rout.save’ ... OK Running ‘cv.R’ Comparing ‘cv.Rout’ to ‘cv.Rout.save’ ... OK Running ‘cv3d.R’ [3s/12s] Comparing ‘cv3d.Rout’ to ‘cv3d.Rout.save’ ... OK Running ‘fit.R’ Comparing ‘fit.Rout’ to ‘fit.Rout.save’ ... OK Running ‘krige0.R’ [4s/12s] Comparing ‘krige0.Rout’ to ‘krige0.Rout.save’ ... OK Running ‘line.R’ Comparing ‘line.Rout’ to ‘line.Rout.save’ ... OK Running ‘merge.R’ Comparing ‘merge.Rout’ to ‘merge.Rout.save’ ... OK Running ‘na.action.R’ Comparing ‘na.action.Rout’ to ‘na.action.Rout.save’ ... OK Running ‘rings.R’ Comparing ‘rings.Rout’ to ‘rings.Rout.save’ ... OK Running ‘sim.R’ Comparing ‘sim.Rout’ to ‘sim.Rout.save’ ... OK Running ‘stars.R’ [52s/76s] Comparing ‘stars.Rout’ to ‘stars.Rout.save’ ... OK Running ‘variogram.R’ Comparing ‘variogram.Rout’ to ‘variogram.Rout.save’ ... OK Running ‘vdist.R’ Comparing ‘vdist.Rout’ to ‘vdist.Rout.save’ ... OK Running ‘windst.R’ [0m/87m] Running the tests in ‘tests/windst.R’ failed. Complete output: > suppressPackageStartupMessages(library(sp)) > suppressPackageStartupMessages(library(spacetime)) > suppressPackageStartupMessages(library(gstat)) > suppressPackageStartupMessages(library(stars)) > > data(wind) > wind.loc$y = as.numeric(char2dms(as.character(wind.loc[["Latitude"]]))) > wind.loc$x = as.numeric(char2dms(as.character(wind.loc[["Longitude"]]))) > coordinates(wind.loc) = ~x+y > proj4string(wind.loc) = "+proj=longlat +datum=WGS84 +ellps=WGS84" > > wind$time = ISOdate(wind$year+1900, wind$month, wind$day) > wind$jday = as.numeric(format(wind$time, '%j')) > stations = 4:15 > windsqrt = sqrt(0.5148 * wind[stations]) # knots -> m/s > Jday = 1:366 > daymeans = colMeans( + sapply(split(windsqrt - colMeans(windsqrt), wind$jday), colMeans)) > meanwind = lowess(daymeans ~ Jday, f = 0.1)$y[wind$jday] > velocities = apply(windsqrt, 2, function(x) { x - meanwind }) > # match order of columns in wind to Code in wind.loc; > # convert to utm zone 29, to be able to do interpolation in > # proper Euclidian (projected) space: > pts = coordinates(wind.loc[match(names(wind[4:15]), wind.loc$Code),]) > pts = SpatialPoints(pts) > if (require(sp, quietly = TRUE) && require(maps, quietly = TRUE)) { + proj4string(pts) = "+proj=longlat +datum=WGS84 +ellps=WGS84" + utm29 = "+proj=utm +zone=29 +datum=WGS84 +ellps=WGS84" + pts = as(st_transform(st_as_sfc(pts), utm29), "Spatial") + # note the t() in: + w = STFDF(pts, wind$time, data.frame(values = as.vector(t(velocities)))) + + library(mapdata) + mp = map("worldHires", xlim = c(-11,-5.4), ylim = c(51,55.5), plot=FALSE) + sf = st_transform(st_as_sf(mp, fill = FALSE), utm29) + m = as(sf, "Spatial") + + # setup grid + grd = SpatialPixels(SpatialPoints(makegrid(m, n = 300)), + proj4string = m@proj4string) + # grd$t = rep(1, nrow(grd)) + #coordinates(grd) = ~x1+x2 + #gridded(grd)=TRUE + + # select april 1961: + w = w[, "1961-04"] + + covfn = function(x, y = x) { + du = spDists(coordinates(x), coordinates(y)) + t1 = as.numeric(index(x)) # time in seconds + t2 = as.numeric(index(y)) # time in seconds + dt = abs(outer(t1, t2, "-")) + # separable, product covariance model: + 0.6 * exp(-du/750000) * exp(-dt / (1.5 * 3600 * 24)) + } + + n = 10 + tgrd = seq(min(index(w)), max(index(w)), length=n) + pred = krige0(sqrt(values)~1, w, STF(grd, tgrd), covfn) + layout = list(list("sp.points", pts, first=F, cex=.5), + list("sp.lines", m, col='grey')) + wind.pr0 = STFDF(grd, tgrd, data.frame(var1.pred = pred)) + + v = vgmST("separable", + space = vgm(1, "Exp", 750000), + time = vgm(1, "Exp", 1.5 * 3600 * 24), + sill = 0.6) + wind.ST = krigeST(sqrt(values)~1, w, STF(grd, tgrd), v) + + all.equal(wind.pr0, wind.ST) + + # stars: + df = data.frame(a = rep(NA, 324*10)) + s = STF(grd, tgrd) + newd = addAttrToGeom(s, df) + wind.sta = krigeST(sqrt(values)~1, st_as_stars(w), st_as_stars(newd), v) + # 1 + plot(stars::st_as_stars(wind.ST), breaks = "equal", col = sf.colors()) + # 2 + stplot(wind.ST) + # 3 + plot(wind.sta, breaks = "equal", col = sf.colors()) + st_as_stars(wind.ST)[[1]][1:3,1:3,1] + (wind.sta)[[1]][1:3,1:3,1] + st_bbox(wind.sta) + bbox(wind.ST) + all.equal(wind.sta, stars::st_as_stars(wind.ST), check.attributes = FALSE) + + # 4: roundtrip wind.sta->STFDF->stars + rt = stars::st_as_stars(as(wind.sta, "STFDF")) + plot(rt, breaks = "equal", col = sf.colors()) + # 5: + stplot(as(wind.sta, "STFDF")) + st_bbox(rt) + + # 6: + stplot(as(st_as_stars(wind.ST), "STFDF")) + } OMP: Warning #96: Cannot form a team with 24 threads, using 2 instead. OMP: Hint Consider unsetting KMP_DEVICE_THREAD_LIMIT (KMP_ALL_THREADS), KMP_TEAMS_THREAD_LIMIT, and OMP_THREAD_LIMIT (if any are set). Flavor: r-devel-linux-x86_64-fedora-clang