Subsetting from a bed.matrix
Returns subset of individuals satisfying a condition.
select.inds(x, condition)
x |
|
condition |
Condition used to select individuals |
The conditions can involve global variables and all variables defined
in the data frame x@ped
, in particular
famid
, id
, father
, mother
, sex
, pheno
If basic stats have been computed (see set.stats
),
N0
, N1
, N2
, NAs
, callrate
, etc.
If some condition evaluate to NA
(e.g. sex == 1
when sex
is undefined for some individuals),
a warning is issued and the corresponding individuals are removed.
A bed.matrix
similar to x
, containing the selected individuals only
Hervé Perdry and Claire Dandine-Roulland
# Load data data(LCT) x <- as.bed.matrix(LCT.gen, LCT.fam, LCT.bim) # Select individuals with a call rate > 95% # and more than 5% of heterozygous genotypes y <- select.inds(x, callrate > 0.95 & N1/(N0+N1+N2) > 0.05) y
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.