InitialEst: InitialEst 'InitialEst' is a function to build an initial...

Description Usage Arguments Value Examples

View source: R/FastLORS_Functions.R

Description

InitialEst InitialEst is a function to build an initial estimate for B

Usage

1
InitialEst(Y, X, lambda = NULL)

Arguments

Y

gene expression matrix

X

matrix of SNPs

lambda

tuning parameter

Value

B An initial estimate of the coefficient matrix

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
## 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

Init_est <- InitialEst(Y,X)

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