make_bim {genio} | R Documentation |
Create a Plink BIM tibble
Description
This function simplifies the creation of Plink BIM-formatted tibbles, which autocompletes missing information if a partial tibble is provided, or generates a completely made up tibble if the number of individuals is provided. The default values are most useful for simulated genotypes, where IDs can be made up but must be unique, and there are no chromosomes, positions, or particular reference or alternative alleles.
Usage
make_bim(tib, n = NA)
Arguments
tib |
The input tibble (optional).
Missing columns will be autocompleted with reasonable values that are accepted by Plink and other external software.
If missing, all will be autocompleted, but |
n |
The desired number of loci (rows).
Required if |
Details
Autocompleted column values:
-
chr
:1
(all data is on a single chromosome) -
id
:1:n
-
posg
:0
(missing) -
pos
:1:n
-
ref
:1
-
alt
:2
Note that n
is either given directly or obtained from the input tibble.
Value
The input tibble with autocompleted columns and columns in default order, or the made up tibble if only the number of individuals was provided. The output begins with the standard columns in standard order: chr, id, posg, pos, ref, alt. Additional columns in the input tibble are preserved but placed after the standard columns.
See Also
Plink BIM format reference: https://www.cog-genomics.org/plink/1.9/formats#bim
Examples
# create a synthetic tibble for 10 loci
# (most common use case)
bim <- make_bim(n = 10)
# manually create a partial tibble with only chromosomes defined
library(tibble)
bim <- tibble(chr = 0:2)
# autocomplete the rest of the columns
bim <- make_bim(bim)