breakpointDist {Trendy} | R Documentation |
calculates number of breakpoints at each time.
breakpointDist(topTrendyData, NDigits = 0)
topTrendyData |
results from topTrendy() function |
NDigits |
how many digits to be used when rounding (default is 0 (return integers)) |
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.
Ning Leng
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)