-
Notifications
You must be signed in to change notification settings - Fork 0
/
hfp.Rd
executable file
·62 lines (54 loc) · 1.55 KB
/
hfp.Rd
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
\name{hfp}
\alias{hfp}
\title{
Function to construct a heatmap of the hidden variation in the gene
expression data.
}
\description{
The function \code{hfp} produces a plot of the PLS imputed estimate of the
hidden variability in the data, derived from the optimal model, corresponding
to an user-specified set of genes and subjects/samples.
}
\usage{
hfp(obj, gen, ind, Y)
}
\arguments{
\item{obj}{
An \code{svpls} object.
}
\item{gen}{
An user-specified set of genes.
}
\item{ind}{
An user-specified set of subjects.
}
\item{Y}{
A log transformed gene expression matrix with genes along the rows
and subjects/samples along the columns.
}
}
\value{
A heatmap of the hidden variability corresponding to the specified set of
genes and subjects, attributable to the unknown subject-specific
factors in the gene expression data.
}
\references{
Sutirtha Chakraborty, Somnath Datta and Susmita Datta. (2012)
Surrogate Variable Analysis Using Partial Least Squares (SVA-PLS) in Gene Expression Studies. Bioinformatics, 28(6): 799-806.
}
\author{
Sutirtha Chakraborty, Somnath Datta and Susmita Datta.
}
\seealso{
\code{\link{heatmap}}, \code{\link{fitModel}}, \code{\link{svpls}}
}
\examples{
## Fitting the optimal ANCOVA model to the data gives:
data(hidden_fac.dat)
fit <- svpls(10,10,hidden_fac.dat,pmax = 5)
## Specifying the set of genes and subjects
genes <- c(1,20,55,70,100,150,250,450)
subjects <- c(1,4,7,10,11,15,17,20)
hfp(fit,genes,subjects,hidden_fac.dat)
}
\keyword{print}