NewTimeEst: Estimate New Time Points

In rmRNAseq: RNA-Seq Analysis for Repeated-Measures Data

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

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

rmRNAseq documentation built on Nov. 8, 2019, 5:06 p.m.