Search Archives
Search tips
from:Name
Search by author name, e.g. from:Duncan Murdoch
"exact phrase"
Match an exact phrase
word1 word2
Match messages containing both words
Date range
Use the date pickers to filter results to a time period
Use the list dropdown to narrow results to a specific mailing list. Combine from: with other terms to filter by author and content.
Dear List, dear Edzer, is it correct if I use the term "regression kriging" when I run this function?: kriged_value <- gstat::krige(z ~ x + y, [...])@data$var1.pred Or should I call it "kriging with external drift" (or "universal...
Hello Andy, you'll need these: library(rgeos) centroids <- gCentroid(states, byid=TRUE, id = NULL) after you convert your data.frame to SpatialPolygons or SpatialPolygonsDataFrame (package sp). HTH, ?kos 2016.05.20. 5:41 keltez?ssel, Andy Bunn ?rta: > Hello...
Dear Milu, ESRI shapefile format cannot contain field names longer than 10 characters. So you should do one of the following options: - choose other file format when saving the data, - rename the columns of the sf object before saving. HTH...
Dear All, I have a simple question about regression kriging (gstat:::krige()) and variogram fitting (gstat:::fit.variogram()), which I have not found exact answer for (neither in the description of package gstat, nor on the internet). Does it cause...
Dear Jonathan, Function get.snow() is not vectorized, i.e. it does not return a vector of the same length as its input. According to the help of app(), "The function should return a vector or matrix that is divisible...
Dear Ani, Package 'raster' has a function called getData(), which can download SRTM tiles. Then you can simply plot() them. Or you can use package 'elevatr'. But I'm sure there are many other alternative solutions. HTH, ?kos Bede-Fazekas...
Dear Nina, you should make prediction on the evaluation dataset (Note that since kriging is an exact interpolation method, predicting on the original dataset does not make sense): prediction <- krige(dependent_variable ~ predictor, locations = original_dataset, newdata = evaluation_dataset) Then...
Dear Diego, Yes, you calculated the spatial autocorrelation of the predicted values in this way, which is not too informative and intresting. Autocorrelation of the residuals is what you need. If you subtract the predictions from the observations (using the...
Hello, try newplot <- ggplot(data = GEV[1:4, ], aes(x = RL45)) + geom_density(color = "midnightblue") + xlab("Location (mm/day") + ggtitle("Global Location under RCP4.5") + xlim(300, 350) HTH, ?kos ----------- ?kos Bede-Fazekas Centre for Ecological Research, Hungary 2022.10...
Dear Patrick, Package 'spatstat' can generate point pattern from several point processes, e.g. from uniform Poisson process (aka. Complete Spatial Randomness, CSR), general Poisson cluster process and Simple Sequential Inhibition process, and can also generate independent uniform random points...
Dear Patrick, please take a look at sf::st_layers(dsn), this gives you the name of, and basic information on, the layers found in the specific file/folder (dsn). HTH, ?kos Bede-Fazekas 2021.04.13. 8:35 keltez...
Dear Enoch, Without a reproducible example it is hard to guess what you data looks like and what causes the unexpected results. Anyway, if you are sure that raster object called "Trying" has perfect CRS and resolution, then this may...
Dear Manuel, There are some broken links in the WorldClim website. If you find one, you can write to info at worldclim.org. Although I have negative experiences, you should try this way. Anyway, the global, non-tiled raster can...
Dear Thiago, a reproducible example and a desired output would help us to understand and solve the problem. The description was not clear, at least, for me. But it seems to be a simple apply/lapply/aggregate/by task that...
Dear Cristabel, function focal() of package raster is what you are searching for. focal(x, w = matrix(1/9, ncol = 3, nrow = 3), fun = sum) or focal(x, w = matrix(1, ncol = 3, nrow = 3), fun = mean) HTH, ?kos Bede...
Dear Manuel, technically, the result of st_intersection(x, y), where both x and y are POLYGONs, can be POINT, LINESTRING, POLGYON and GEOMETRY as well. The result is GEOMETRY if the type of the different features is not the...
Dear Cristi?n, there migt be several packages that can compute confusion matrix (aka. contingency table): ??contingency One of them that I can recommend you is package 'ROCR' and functions ROCR::prediction() and ROCR::performance(). ## Creates a ROCR::prediction object...
Dear Erin, I think that projectRaster(from = import_raster, to = r) or crop(x = import_raster, y = r) is what you need. crop() simply crop x with the extent of y, while projectRaster(from, to) create a raster from from...
Dear Laura, Below are two examples that you may find relevant: https://github.com/TomNash/shinySDM https://github.com/wallaceEcoMod/wallace Have a nice week, ?kos ____________ ?kos Bede-Fazekas Centre for Ecological Research, Hungary 2024. 10. 15. 20:55 keltez...
Dear Milu, I think that you are looking for as.data.frame(x, xy = TRUE). HTH, ?kos Bede-Fazekas Hungarian Academy of Sciences 2018.08.14. 17:03 keltez?ssel, Miluji Sb ?rta: > Dear all, > > I have the following rasterbrick...
Can't find what you're looking for? Try searching with Google .