intern.chromo {globalSeq} | R Documentation |
Communicates between cursus
and omnibus
by coordinating a chromosome-wide analysis.
intern.chromo(Y, Ystart, Yend, X, Xloc, window, offset, group, perm, nodes, phi, kind)
Y |
RNA-Seq data:
numeric matrix with |
Ystart |
start location of genes:
numeric vector of length |
Yend |
end location of genes:
NULL or numeric vector of length |
X |
genomic profile:
numeric matrix with |
Xloc |
location covariates:
numeric vector of length |
window |
maximum distance: non-negative real number |
offset |
numeric vector of length |
group |
confounding variable:
factor of length |
perm |
number of iterations: positive integer |
nodes |
number of cluster nodes for parallel computation |
phi |
dispersion parameters: vector of length |
kind |
computation : number between 0 and 1 |
The function returns a dataframe, with the p-value in the first column, and the test statistic in the second column.
# simulate high-dimensional data n <- 30 q <- 10 p <- 100 set.seed(1) Y <- matrix(rnbinom(q*n,mu=10, size=1/0.25),nrow=q,ncol=n) X <- matrix(rnorm(p*n),nrow=p,ncol=n) Yloc <- seq(0,1,length.out=q) Xloc <- seq(0,1,length.out=p) window <- 1 # hypothesis testing cursus(Y,Yloc,X,Xloc,window)