missingPersonPlot {forrel} | R Documentation |
Missing person plot
Description
Visualises the competing hypotheses of a family reunion case. A plot with two panels is generated. The left panel shows a pedigree in which the person of interest (POI) is identical to the missing person (MP). The right panel shows the situation where these two are unrelated. See Details for further explanations.
Usage
missingPersonPlot(
reference,
missing,
labs = labels(reference),
marker = NULL,
hatched = typedMembers(reference),
MP.label = "MP",
POI.label = "POI",
POI.sex = getSex(reference, missing),
POI.col = "red",
POI.hatched = FALSE,
POI.height = 8,
titles = c(expression(H[1] * ": POI = MP"), expression(H[2] * ": POI unrelated")),
width = NULL,
cex = 1.2,
newdev = interactive(),
...
)
Arguments
reference |
A |
missing |
The ID label of the missing pedigree member. |
labs |
A character vector with labels for the pedigree members. See
|
marker |
Optional vector of marker indices to be included in the plot. |
hatched |
A vector of ID labels indicating who should appear with hatched symbols in the plot. By default, all typed members. |
MP.label |
The label of the missing member. Default: "MP". |
POI.label |
The label of the person of interest. Default: "POI". |
POI.sex |
The sex of POI. This defaults to that of the missing person, but may be set explicitly. This is particularly useful when the missing person has unknown sex. |
POI.col |
The plot colour of POI. Default: red. |
POI.hatched |
A logical: If TRUE (default), the POI is plotted with a hatched symbol. |
POI.height |
A numeric controlling the vertical placement of the POI singleton (in the right panel). |
titles |
A character of length 2, with subtitles for the two frames. |
width |
A positive number controlling the width of the plot. More specifically this number is the relative width of the reference pedigree, compared to a singleton. |
cex |
Expansion factor for pedigree symbols and font size. |
newdev |
A logical: If TRUE the plot is created in a new plot window. |
... |
Extra parameters passed on to |
Details
A standard family reunification case involves the following ingredients:
A reference family in which a single member ("MP") is missing.
Some of the family members have been genotyped
A person of interest ("POI") is to be matched against the reference family
After genotyping of POI, the genetic evidence is typically assessed by computing the likelihood ratio of the following hypotheses:
H1: POI is MP
H2: POI is unrelated to the family
The goal of this function is to illustrate the above hypotheses, using labels, colours and shading to visualise the different aspects of the situation.
This function cannot handle cases with more complicated hypotheses (e.g.
multiple missing persons, or where H2 specifies a different relationship).
However, as it is basically a wrapper of pedtools::plotPedList()
, an
interested user should be able to extend the source code to such cases
without too much trouble.
Value
None
Examples
x = nuclearPed(father = "fa", mother = "mo", children = c("b1", "b2"))
# Default plot
missingPersonPlot(x, missing = "b2")
# Exploring various options
missingPersonPlot(x,
missing = "b2",
hatched = "b1",
POI.hatched = TRUE,
POI.sex = 0,
cex = 1.5, # larger symbols and label font (see ?par())
cex.main = 1.3, # larger frame titles (see ?par())
dev.width = 7, # device width (see ?plotPedList())
dev.height = 3 # device height (see ?plotPedList())
)