read_phen {genio} | R Documentation |
Read *.phen files
Description
This function reads a standard *.phen file into a tibble.
It uses readr::read_table()
to do it efficiently.
GCTA and EMMAX use this format.
Usage
read_phen(file, verbose = TRUE)
Arguments
file |
Input file (whatever is accepted by |
verbose |
If |
Value
A tibble with columns: fam
, id
, pheno
.
See Also
GCTA PHEN format reference: https://cnsgenomics.com/software/gcta/#GREMLanalysis
Examples
# to read "data.phen", run like this:
# phen <- read_phen("data")
# this also works
# phen <- read_phen("data.phen")
# The following example is more awkward
# because package sample data has to be specified in this weird way:
# read an existing plink *.phen file
file <- system.file("extdata", 'sample.phen', package = "genio", mustWork = TRUE)
phen <- read_phen(file)
phen
# can specify without extension
file <- sub('\\.phen$', '', file) # remove extension from this path on purpose
file # verify .phen is missing
phen <- read_phen(file) # load it anyway!
phen