order_mm_tapered: Maximin ordering using a tapered covariance matrix

In katzfuss-group/NPvecchia: Large Scale Covariance Estimation

Maximin ordering using a tapered covariance matrix

Description

This function is a wrapper for order_maximin_dist to take in data and locations to calculate the ordering based on a tapered sample covariance matrix. It provides an argument to control the amount of tapering.

Usage

order_mm_tapered(locs, datum, tapering_range = 0.4)

Arguments

locs	matrix of locations of points (to match input argument of fields::rdist)
datum	Data where column i corresponds to observations at the i'th point of locs
tapering_range	Percentage of the maximum distance for Exponential tapering, which defaults to 0.4 * the maximum distance.

Value

Maximin ordering of the points

Examples

n <- 100
d <- 2
N <- 50
locationss <- matrix(runif(n * d), nc = d)
dataa <- matrix(rnorm(n * N), nr = N)
order_mm <- order_mm_tapered(locationss, dataa)
#reorder locations/data
locationss <- locationss[order_mm, ]
dataa <- dataa[, order_mm]

katzfuss-group/NPvecchia documentation built on April 15, 2022, 2:23 a.m.