plotsibs: A function to plot sibships (full and half-sibs) in an x-y...

In rcolony: Tools for implementation and visualisation of Colony2 and its outputs.

Description

This function plots the estimated sibships from Colony data. Both the x and y axis represent individuals in the population and points indicate a sibship assignment between individuals x and y.

Usage

plotsibs(colony.object, maintitle = "sibships", prob = 0.8, pairwise = FALSE)

Arguments

colony.object	A list object derived from the colony data (get.colony.data).
maintitle	The plot's title.
prob	Assignment probability threshold. The probability of sibship must be greater than this value to be included in the plot.
pairwise	TRUE/FALSE: should pairwise sibship estimates be used. (TRUE= pairwise, FALSE=full likelihood)

Details

Full sibships are shown in the upper left part of the plot, and half-sibships are shown in the lower right part. The data can be derived from pairwise methods, or from the full likelihood method.

Value

A plot showing full- or half-sibships.

Author(s)

Owen R. Jones

References

Wang, J. (2004) Sibship reconstruction from genetic data with typing errors. Genetics 166: 1963-1979.

See Also

get.colony.data,run.colony

Examples

data(testdata)

plotsibs(testdata,maintitle="Sibships (pairwise)",pairwise=TRUE)
plotsibs(testdata,maintitle="Sibships (full)",pairwise=FALSE)

rcolony documentation built on May 2, 2019, 6:42 p.m.