| byclade {LOST} | R Documentation |
Simulate missing morphometric data with taxonomic bias
Description
This function simulates higher frequency of missing data points in groups that are less numerically well represented in the whole sample, relative to other group. These groups may represent taxa (as used in Brown et al., 2012), but may also represent any other group of interest (e.g. populations, trials, subsamples, etc.). From a morphometric dataset, this function first selects a number of specimens to have data points removed from at random. A vector containing the number of measurements to remove from each specimen is sorted into descending order. Specimens are then sampled without replacement with a probability relative to the sum of the entire sample sizes divided by the number of specimens its respective group. The order the specimens are sampled determines the number of data points to be removed (i.e. the first to be sampled has the most removed). A complete mathematical description may be found in Brown et al. (2012).
Usage
byclade(x, remperc , groups)
Arguments
x |
A n X m matrix of morphometric data with n specimens and m variables. Or an l X 2 or 3 X n array of geometric morphometric coordinates (2D or 3D), where l is the number of landmarks. |
remperc |
The percentage of data to be removed from the matrix, expressed as a decimal (ex: 30 percent would be entered as 0.3) |
groups |
A vector of length n specifying taxonomic group membership as integers (ex: c(1,1,2,2,3,3,...) ) |
Value
returns a matrix or array (depending on input) of morphometric data with missing variables input as 'NA'
Author(s)
J. Arbour and C. Brown
References
Brown, C., Arbour, J. and Jackson, D. 2012. Testing of the Effect of Missing Data Estimation and Distribution in Morphometric Multivariate Data Analyses. Systematic Biology 61(6):941-954.