breakpointDist: Distribution of breakpoints

In rhondabacher/Trendy: Breakpoint analysis of time-course expression data

Distribution of breakpoints

Description

calculates number of breakpoints at each time.

Usage

breakpointDist(topTrendyData, NDigits = 0)

Arguments

topTrendyData	results from topTrendy() function
NDigits	how many digits to be used when rounding (default is 0 (return integers))

Value

The function takes significant genes called from the topTrendyData() function. For any time point, this function calculates how many genes have a breakpoint at this time point. The output is the numbers of genes sorted by time point.

Author(s)

Ning Leng

Examples

 m1 <- matrix(c(c(rnorm(50,5,1),sort(rnorm(50, 15, 5))), rnorm(100, 50,10)), 2, 100, TRUE)
 rownames(m1) <- c("g1","g2")
 colnames(m1) <- paste0("time", seq_len(100))
 myTrends <- results(trendy(m1))
 topGenes <- topTrendy(myTrends)
 bpDist <- breakpointDist(topGenes)

rhondabacher/Trendy documentation built on Oct. 26, 2023, 3:46 a.m.