rmRNAseq (version 0.1.0)

NewTimeEst: Estimate New Time Points

Description

This function estimate new time points to fit the glsCAR1 for the voom transformed data. Note that this function is very specific for this dataset, with only 4 time points. If there are more than 4 time points, the method needs to be updated.

Usage

NewTimeEst(v, Subject, Time, TimeMinOut, ncores)

Arguments

v

output of voom function.

Subject

a vector of subjects or experimental units.

Time

a vector of time points.

TimeMinOut

output from the TimeMin function

ncores

number of cores for embarrassingly parallel procedure. Default value of ncores is 1.

Value

New time points.

Examples

# NOT RUN {
data(res)
data(covset)
v <- res$ori.res$v
v$E <- v$E[1:2,]
v$weights <- v$weights[1:2,]
Subject <- covset$ear
Time <- covset$time
ncores <- 1
NewTimeOut <- rmRNAseq:::NewTimeEst(v, Subject, Time, ncores)
NewTimeOut
# }