1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
|
% Generated by roxygen2: do not edit by hand
% Please edit documentation in R/testPseudotime.R
\name{testPseudotime}
\alias{testPseudotime}
\alias{testPseudotime,ANY-method}
\alias{testPseudotime,SummarizedExperiment-method}
\title{Test for differences along pseudotime}
\usage{
testPseudotime(x, ...)
\S4method{testPseudotime}{ANY}(
x,
pseudotime,
df = 5,
get.lfc = TRUE,
get.spline.coef = FALSE,
trend.only = TRUE
)
\S4method{testPseudotime}{SummarizedExperiment}(x, ..., assay.type = "logcounts")
}
\arguments{
\item{x}{A numeric matrix-like object containing log-expression values for cells (columns) and genes (rows).
Alternatively, a \linkS4class{SummarizedExperiment} containing such a matrix.}
\item{...}{For the generic, further arguments to pass to specific method.
For the ANY method, further arguments to pass to \code{\link{fitLinearModel}}.
For the SummarizedExperiment method, further arguments to pass to the ANY method.}
\item{pseudotime}{A numeric matrix with one row per cell in \code{x} and one column per path (i.e., lineage).
A vector is treated the same as a 1-column matrix.}
\item{df}{Integer scalar specifying the degrees of freedom for the splines.}
\item{get.lfc}{Logical scalar indicating whether to return an overall log-fold change along each path.}
\item{get.spline.coef}{Logical scalar indicating whether to return the estimates of the spline coefficients.}
\item{trend.only}{Logical scalar indicating whether only differences in the trend should be considered
when testing for differences between paths.}
\item{assay.type}{String or integer scalar specifying the assay containing the log-expression matrix.}
}
\value{
A \linkS4class{DataFrame} is returned containing the statistics for each gene (row),
including the p-value and its BH-adjusted equivalent.
If \code{get.lfc=TRUE}, an overall log-fold change is returned for each path.
If \code{get.spline.coef=TRUE}, the estimated spline coefficients are also returned (single path)
or the differences in the spline fits to the first path are returned (multiple paths).
}
\description{
Implements a simple method of testing for significant differences with respect to pseudotime,
based on fitting linear models with a spline basis matrix.
This function is now deprecated as it has been moved to the \pkg{TSCAN} package itself.
}
\details{
For a single path in \code{pseudotime},
this function fits a natural spline to the expression of each gene with respect to pseudotime.
It then does an ANOVA to test whether any of the spline coefficients are non-zero.
In this manner, genes exhibiting a significant (and potentially non-linear) trend
with respect to the pseudotime can be detected as those with low p-values.
For multiple paths in \code{pseudotime},
the null hypothesis is that all paths have the same trend (if \code{trend.only=TRUE})
or the same trend and intercept (if \code{FALSE}).
This is done by effectively fitting a separate trend to each path
and performing an ANOVA to detect differences in the trend alone or in the trend and intercept.
In this manner, genes exhibiting differences in behavior between paths can be detected.
The expected format of \code{pseudotime} is the same as that returned by \code{\link{orderClusterMST}}.
Each cell is assigned to a path if it has a non-\code{NA} value in the corresponding column.
For single path testing, cells with \code{NA} values in \code{pseudotime} are ignored;
for multiple path testing, cells assigned to multiple paths are ignored.
By default, estimates of the spline coefficients are not returned as they are difficult to interpret.
Rather, a log-fold change of expression along each path is estimated
to provide some indication of the overall magnitude and direction of any change.
}
\examples{
y <- matrix(rnorm(10000), ncol=100)
# Testing for a difference in a single path:
u <- runif(100)
testPseudotime(y, u)
# Testing for differences in multiple paths
# by mocking up a pseudotime matrix.
p <- cbind(path1=u, path2=u)
path1 <- rbinom(length(u), 1, 0.5)==0
p[!path1,1] <- NA
p[path1,2] <- NA
testPseudotime(y, p)
}
\seealso{
\code{\link{orderClusterMST}}, to generate the pseudotime matrix.
\code{\link{testLinearModel}}, which performs the tests under the hood.
}
\author{
Aaron Lun
}
|
