pleiotropyh2_nocor_computing_single: Compute pleioh2g before bias correction for target disease

In pleioh2g: Estimation of Pleiotropic Heritability from Genome-Wide Association Studies (GWAS) Summary Statistics

Compute pleioh2g before bias correction for target disease

Description

This function is used to compute pleioh2g after bias correction for target disease

Usage

pleiotropyh2_nocor_computing_single(
  G,
  phenotype,
  h2_vector,
  h2_vector_mat,
  Results_full_rg,
  Results_full_rg_array
)

Arguments

G	index of target disease.
phenotype	Vector of the phenotype name
h2_vector	h2g vector for all traits - aligned as the order in phenotype file
h2_vector_mat	h2g array from jackknife-block estimates for all traits - aligned as the order in phenotype file
Results_full_rg	genetic correlation matrix.- aligned as the order in phenotype file
Results_full_rg_array	genetic correlation jackknife-block array.- aligned as the order in phenotype file

Value

A 'list' containing the following elements: - 'target_disease' (character): The value "401.1". - 'target_disease_h2_est' (numeric): target disease h2g. - 'target_disease_h2_se' (numeric): target disease h2g_se. - 'selected_auxD' (character): auxiliary diseases. - 'h2pleio_uncorr' (numeric): pre-correction pleiotropic heritability estimate. - 'h2pleio_uncorr_se' (numeric): pre-correction pleiotropic heritability jackknife s.e. estimate. - 'percentage_h2pleio_uncorr' (numeric): pre-correction percentage of pleiotropic heritability estimate. - 'percentage_h2pleio_uncorr_se' (numeric): pre-correction percentage of pleiotropic heritability jackknife s.e. estimate. - 'percentage_h2pleio_jackknife_uncorr' (numeric): vector of all pre-correction percentage of pleiotropic heritability jackknife estimates.

Examples

G <- 1
data(Results_full_rg)
data(Results_full_rg_array)
data(h2_vector)
data(h2_vector_mat)
Results_full_rg<-Results_full_rg[1:15,1:15]
Results_full_rg_array<-Results_full_rg_array[1:15,1:15,]
h2_vector<-t(as.matrix(h2_vector[1,1:15]))
h2_vector_mat<-h2_vector_mat[,1:15]
phenotype<-c("401.1","244.5","318","735.3","411.4",
"427.2","454.1","278.1","250.2","550.1","530.11",
"296.22","519.8","562.1","763")
h2pleiobeforecorr<-pleiotropyh2_nocor_computing_single(G,phenotype,h2_vector,
h2_vector_mat,Results_full_rg,Results_full_rg_array)

pleioh2g documentation built on March 9, 2026, 5:07 p.m.