random.pm: Random square definite positive matrix

In gaston: Genetic Data Handling (QC, GRM, LD, PCA) & Linear Mixed Models

Random square definite positive matrix

Description

Generate a random definite positive matrix with specified dimension

Usage

 random.pm(n, values) 

Arguments

n	Dimension of matrix
values	(Optional) A numeric vector of dimension n : the eigenvalues of the matrix

Details

If values isn't given, it is chosen (deterministically) so that the eigenvalues of the resulting matrix are similar to eigenvalues observed on Genetic Relationship Matrices.

The random matrix is generated as U %*% diag(values) %*% t(U) with U a random orthogonal matrix.

Value

A named list with members:

K	A n x n symmetric positive matrix
eigen	The eigen decomposition of K as eigen(K) would output it

See Also

lmm.simu, eigen

Examples

# generate a random positive matrix 
set.seed(1)
R <- random.pm(500)
str(R)

gaston documentation built on Dec. 28, 2022, 1:30 a.m.