This function writes a tibble with the right columns into a standard *.phen file.
readr::write_tsv() to do it efficiently.
GCTA and EMMAX use this format.
Output file (whatever is accepted by
The tibble or data.frame to write.
It must contain these columns:
tib invisibly (what
GCTA PHEN format reference: https://cnsgenomics.com/software/gcta/#GREMLanalysis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21
# to write an existing table `phen` into file "data.phen", run like this: # write_phen("data", phen) # this also works # write_phen("data.phen", phen) # The following example is more detailed but also more awkward # because (only for these examples) the package must create the file in a *temporary* location # create a dummy tibble with the right columns library(tibble) tib <- tibble( fam = 1:3, id = 1:3, pheno = 1 ) # a dummy file file_out <- tempfile('delete-me-example', fileext = '.phen') # will also work without extension # write the table out in *.phen format (no header, columns in right order) write_phen(file_out, tib) # delete output when done file.remove(file_out)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.