Questions tagged [geor]

An R package to conduct geostatistical analysis.

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…
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…
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…
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
1
2 3