NewTimeEst: Estimate New Time Points

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

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

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
tmOut <- rmRNAseq:::TimeMin(v, Subject, Time, ncores)
TimeMinOut <-tmOut$MinMaxTime
NewTimeOut <- rmRNAseq:::NewTimeEst(v, Subject, Time, TimeMinOut, ncores)
NewTimeOut

ntyet/rmRNAseq documentation built on July 7, 2023, 11:10 a.m.