pleio_test function

Sequential Wald test for pleiotropy

Sequential Wald test for pleiotropy

Performs the sequential test of pleiotropic effects using results of pleioR().

pleio_test( pleio_results, loop_breaker = 1, save_at = NULL, contrast_matrices_list = NULL )

Arguments

  • pleio_results: pleio_class object returned by pleioR().
  • loop_breaker: numeric value for a maximum p-value used to stop the sequence if a higher p-value is obtained. This saves computation time if there are many tests to perform.
  • save_at: character with directory and/or file name (.rdata) to save the results. This is useful when handling multiple results such as in parallel jobs.
  • contrast_matrices_list: user-specified contrast matrices within a list of lists, or a single contrast matrix (see example). Each matrix must have the same number of columns, and must be equal to the number of traits.

Returns

list of p-values, indices, and trait numeric identifier.

Examples

# Example of user-specified contrast matrices with 3 traits cm1 <- matrix(c(-1, 0, 1), ncol = 3) cm2 <- matrix(c(0, -1, 1), ncol = 3) contrast_matrices <- list('1vs3' = list(cm1), '2vs3' = list(cm2)) # or a single contrast matrix as: contrast_matrices <- cm1