Display geographic maps of ancestry coefficients

This function displays interpolated values of ancestry coefficients on geographic maps.

# S3 method for tess3Q
plot(x, coord, method = "map.max", resolution = c(300,
  300), window = NULL, background = TRUE, map.polygon = NULL,
  raster.filename = NULL, interpolation.model = FieldsKrigModel(10),
  col.palette = CreatePalette(), ...)

Arguments

x	an object of class tess3Q containing an ancestry coefficient matrix computed from tess3.
coord	a numeric matrix of dimension \(n\) times 2 where \(n\) is the number of individuals. The matrix must contain (Longitude, Latitude) coordinates for all individuals.
method	a character string "map.max" or "map.all". If "map.all", interpolating surfaces are displayed for the ancestry coefficients of all populations. If "map.max" the union of interpolating surfaces is displayed.
resolution	an integer vector resolution = c(rx,ry) for the resolution of the grid used to computed the interpolating surface. rx and ry are resolution numbers for the x-axis and y-axis respectively.
window	an integer vector for the plotting window, such that window = c(xmin, xmax, ymin, ymax) contains the window's min and max coordinates.
background	if TRUE a raster file is used as a stencil to render only raster pixel on earth.
map.polygon	an object of class sp::SpatialPolygonsDataFrame used to crop the interpolating surface. If NULL, the function getMap is used.
raster.filename	a raster file name used to compute the background stencil. This is an alternative method to crop the interpolating surfaces. The default method uses map.polygon.
interpolation.model	an interpolation model used to compute the interpolating surface. Interpolation models can use the functions FieldsTpsModel or FieldsKrigModel.
col.palette	a list of color palettes. Color palettes can be defined by using the function CreatePalette.
...	plot.default other parameters.

Value

None

Examples

library(tess3r)

# Retrieve a dataset
data(data.at)

# Run of TESS3
obj <- tess3(data.at$X, coord = data.at$coord, K = 5,
                 ploidy = 1, method = "projected.ls", openMP.core.num = 4)
#> == Computing spectral decomposition of graph laplacian matrix: done
#> ==Main loop with 4 threads: done

# Get the ancestry matrix
Q.matrix <- qmatrix(obj, K = 5)

# Plot the spatial interpolation of the ancestry matrix
plot(Q.matrix, data.at$coord, method = "map.max",
     resolution = c(400,400),
     interpolation.model = FieldsKrigModel(10), cex = .4,
     xlab = "Longitude", ylab= "Latitude", main = "Ancestry coefficients")
#> Loading required package: maps
#> This function required to attach maps namespace.