ggtess3Q: Display geographic maps of ancestry coefficients using the...

View source: R/plotQ.R

ggtess3QR Documentation

Display geographic maps of ancestry coefficients using the ggplot grammar

Description

This function displays geographic maps of ancestry coefficients using the ggplot syntax.

Usage

ggtess3Q(Q, coord, resolution = c(300, 300), window = NULL,
  background = TRUE, map.polygon = NULL, raster.filename = NULL,
  interpolation.model = FieldsKrigModel(10), col.palette = CreatePalette())

Arguments

Q

an object of class tess3Q corresponding to an ancestry coefficient matrix obtained 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.

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 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.

Value

None

Author(s)

Kevin Caye, Flora Jay, Olivier François

Examples

library(tess3r)

# Load Arabidopsis data
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)

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

# Show a spatial interpolation of the ancestry matrix
ggtess3Q(Q.matrix, data.at$coord)

bcm-uga/TESS3_encho_sen documentation built on June 30, 2023, 3:08 a.m.