well, this is just a reminder and a place to look for the exact syntax
require(tidyverse)
require(ggplot2)
peds=read_table2("renadd05.ped",
col_names=c("nid","ns","nd","varMS","yob",
"numberKnownParents","numberOfRecords",
"progenyAsSire","progenyAsDam","origId"),
col_types="iiidiiiiic"
)
# this file has BLUP_noUPG and BLUP_MF already joined
rels=read_table2("acc_bf90_both",
col_names=TRUE,
col_types="iiiddiiidd"
)
rels$nid=rels$levelNoUPG
relpeds= rels %>% inner_join(peds,by="nid")
write_delim(relpeds,"acc_bf90_both_ped")
No comments:
Post a Comment