An R package to conduct geostatistical analysis.
Questions tagged [geor]
40 questions
5
votes
2 answers
R package geoR won't install on macOS despite available tcl-tk on the system
I've installed R and RStudio via brew. For a project I need a function from the R library geoR which in turn is dependent on tcl-tk. I installed this as well via brew: brew install tcl-tk. However, trying to install geoR still leads to the error…

Sander W. van der Laan
- 790
- 8
- 26
5
votes
1 answer
Calculate variogram of raster data with NAs in R
Summary: I have a raster dataset which contains NA values, and want to calculate a variogram of it, ignoring the NAs. How can I do this?
I have an image which I have loaded into R using the readGDAL function, stored as im. To make this reproducible,…

robintw
- 27,571
- 51
- 138
- 205
4
votes
2 answers
How do I estimate a variogram for data on a globe?
I have data with associated longitudes and latitudes. How do I get a variogram for this data based on the great-circle distances between the points?
This simple example has all the data on the equator:
require(geoR)
long <- seq(-179, 180)
x <-…

pete
- 2,327
- 2
- 15
- 23
4
votes
0 answers
Create heatmap with distribution of attribute values in R (not density heatmap)
Some of you might have seen Beyond "Soda, Pop, or Coke". I am facing a similar problem and would like to create a plot like that. In my case, I have a very large number of geo-coded observations (over 1 million) and a binary attribute x. I would…

user2503795
- 4,035
- 2
- 34
- 49
3
votes
1 answer
Time taken to krige in gstat package in R
The following R program creates an interpolated surface using 470 data points using walker Lake data in gstat package.
source("D:/kriging/allfunctions.r") # Reads in all functions.
source("D:/kriging/panel.gamma0.r") # Reads in…

Chandan
- 764
- 2
- 8
- 21
2
votes
1 answer
rcpp - R session Aborted
library(Rcpp)
library(geoR)
elevationd=as.matrix(data.frame(xcoords=elevation$coords[,1],
ycoords=elevation$coords[,2], elev=elevation$data))
elevationd
cppFunction('void a(NumericMatrix data){
int nr = data.nrow();
int nc = data.ncol();
…

Jiwon Kim
- 139
- 3
- 12
2
votes
0 answers
Cross-validation for kriging in R: how to include the trend while reestimating the variogram using xvalid?
I have a question very specific for the function xvalid (package geoR) in R which is used in spatial statistics only, so I hope it's not too specific for someone to be able to answer. In any case, suggestions for alternative functions/packages are…

Vera
- 29
- 4
2
votes
0 answers
How to get geoR and sp R packages working in AWS Linux EC2 AMI
My end goal is to be able to use the geoR and sp packages in R on an AWS EC2 Linux instance. The first hurdle is getting an older version of R since 3.3 doesn't support TCL/TK (at least that I've found).
After unpacking R-3.2.5.tar.gz and following…

GISUser9
- 125
- 2
- 8
2
votes
1 answer
Why are my predicted values constant?
mydata
mydata_x <- rep(seq(200,600,by=200),times=1,each=4)
mydata_y <- rep(seq(600,1200,by=200),times=4,each=1)
mydata_z <- c(0,0,0,0,0,0,0,529,0,0,0,0)
mydata <- data.frame(x=mydata_x,y=mydata_y,z=mydata_z)
I've taken samples from a virtual field…

braxtonlewis
- 41
- 1
- 8
2
votes
1 answer
covariance matrix
I was wondering if any one could explain to me how the geoR package calculates the covariance function? I mean how you would do it by hand?
library(geoR)
#suppose I have the following coordinates
X = c(60,30,20,40)
Y = c(50,20,50,50)
my_coordinates…

ToNoY
- 1,358
- 2
- 22
- 43
1
vote
0 answers
How do I fit directional variogram on gridded data, and for a given azimuth return the estimate range?
I have a digital elevation model of topographical features. My overall objective is to make a predictive model by regressing on the range of a directional variogram for a given azimuth. Are there any good python libraries which can accommodate this…

F.Dunbar
- 37
- 6
1
vote
0 answers
My variogram code result different from variog() result
I am writing code for producing a variogram. For validating my result, I checked with geoR::variog() but both variograms are different.
I tried to understand the code of variog() to see what happens under the hood but there are so many things…

Prakkash Manohar
- 11
- 2
1
vote
1 answer
Subsetting R lists according to the value of a function in geoR
I'm looking into the behavior of R, and figured out something wierd.
I will use the package geoR as a reference, since it is the package I am working with.
I seem not to grasp how subsetting lists work. I have the following…

Easymode44
- 148
- 9
1
vote
0 answers
Boxcoxfit function in R without intecept
I have a problem with boxcoxfit function.
I simulated some data and now I want to create estimators for regression parameters and parameter in box-cox transformation.
I use a package geoR
I have a matrix X with 2 columns and Y with not-negative…

Pavel Pejřimovský
- 11
- 1
1
vote
2 answers
R can't find function even though I installed its package. .nlmP (from geoR)
I write
install.packages("geoR")
library(geoR)
.nlmP(f,parameters)
and it returns
Error: could not find function ".nlmP"
What to do?
I did what the first answer suggest, then it said
Error in get(".ind.prof.phi", pos = 1) : object '.ind.prof.phi'…

Jaood
- 385
- 1
- 2
- 12