genio (version 1.0.9)

write_bed: Write a genotype matrix into plink BED format

Description

This function accepts a standard R matrix containing genotypes (values in c(0,1,2,NA)) and writes it into a plink-formatted BED (binary) file. Each genotype per locus (m loci) and individual (n total) counts the number of alternative alleles or NA for missing data. No *.fam or *.bim files are created by this basic function.

Usage

write_bed(file, X, verbose = TRUE)

Arguments

file

Output file path. .bed extension may be omitted (will be added automatically if it is missing).

X

The \(m \times n\) genotype matrix. Row and column names, if present, are ignored.

verbose

If TRUE (default) function reports the path of the file being written (after autocompleting the extension).

Value

Nothing

Details

Genotypes with values outside of \([0,2]\) cause an error, in which case the partial output is deleted. However, beware that decimals get truncated internally, so values that truncate to 0, 1, or 2 will not raise errors. The BED format does not accept fractional dosages, so such data will not be written as expected.

See Also

write_plink for writing a set of BED/BIM/FAM files.

Plink BED format reference: https://www.cog-genomics.org/plink/1.9/formats#bed

Examples

# NOT RUN {
file_out <- tempfile('delete-me-example', fileext = '.bed') # will also work without extension
# create 10 random genotypes
X <- rbinom(10, 2, 0.5)
# replace 3 random genotypes with missing values
X[sample(10, 3)] <- NA
# turn into 5x2 matrix
X <- matrix(X, nrow = 5, ncol = 2)
# write this data to file in BED format
# (only *.bed gets created, no *.fam or *.bim in this call)
write_bed(file_out, X)
# delete output when done
file.remove(file_out)

# }