GetMaxRho: GetMaxRho

In jdrhyne2/FastLORS: Joint Modeling for eQTL Mapping in R

Description

GetMaxRho is a function used to determine the maximum value as a candidate for rho. See the parmeter tuning section of Yang et al. (2013) Note: This function is adapted from the LORS MATLAB implementation

Usage

GetMaxRho(X, Y, L, Omega0)

Arguments

X	matrix of SNPs
Y	gene expression matrix
L	matrix of hidden factors
Omega0	Boolean matrix of observed entries

Value

MaxRho The maximum value of rho to be used in parameter tuning

Examples

##Example

## Generate some data
n <- 50
p <- 200
q <- 100
k <- 10
set.seed(123)
X <- matrix(rbinom(n*p,1,0.5),n,p)
L <- matrix(rnorm(n*k),n,k) %*% t(matrix(rnorm(q*k),q,k))
B <- matrix(0, ncol(X), ncol(L))
activeSNPs <- sort(sample(c(1:nrow(B)), 20))
for(i in 1:length(activeSNPs)){
genes_influenced <- sort(sample(c(1:ncol(B)),5))
B[activeSNPs[i], genes_influenced] <- 2
}
E <- matrix(rnorm(n*q),n,q)
Y <- X %*% B + L + E
Omega0 <- !(is.na(Y))

## Usage
GetMaxRho(X, Y, L, Omega0)

jdrhyne2/FastLORS documentation built on March 5, 2020, 6:50 a.m.