Constructing covariance matrix from local family graph for multi trait analysis — graph_based_covariance_construction

Function that constructs the genetic covariance matrix given a graph around a proband and extracts the threshold information from the graph.

graph_based_covariance_construction_multi(
  fam_id,
  pid,
  cur_proband_id,
  cur_family_graph,
  h2_vec,
  genetic_corrmat,
  phen_names,
  add_ind = TRUE
)

Arguments

fam_id: Name of column with the family ID
pid: Name of column of personal ID
cur_proband_id: id of proband
cur_family_graph: local graph of current proband
h2_vec: vector of liability scale heritabilities
genetic_corrmat: matrix with genetic correlations between considered phenotypes. Must have same order as h2_vec.
phen_names: Names of the phenotypes, as given in cur_family_graph.
add_ind: whether to add genetic liability of the proband or not. Defaults to true.

Value

list with three elements. The first element is temp_tbl, which contains the id of the current proband, the family ID and the lower and upper thresholds for all phenotypes. The second element, cov, is the covariance matrix of the local graph centred on the current proband. The third element is newOrder, which is the order of ids from pid and phen_names pasted together, such that order can be enforced elsewhere too.