-
Notifications
You must be signed in to change notification settings - Fork 0
/
MultiVariAnalysis.Rd
84 lines (80 loc) · 2.15 KB
/
MultiVariAnalysis.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
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
\name{MultiVariAnalysis}
\alias{MultiVariAnalysis}
%- Also NEED an '\alias' for EACH other topic documented here.
\title{
Function for Multivariant Analysis
}
\description{
%% ~~ A concise (1-5 lines) description of what the function does. ~~
}
\usage{
MultiVariAnalysis(object, method = "clustering", ...)
}
%- maybe also 'usage' for other objects documented here.
\arguments{
\item{object}{
preprocessing data from gcms, lcms and rnaseq
}
\item{method}{
clustering methods
}
\item{\dots}{
%% ~~Describe \code{\dots} here~~
}
}
\details{
%% ~~ If necessary, more details than the description above ~~
}
\value{
%% ~Describe the value returned
%% If it is a LIST, use
%% \item{comp1 }{Description of 'comp1'}
%% \item{comp2 }{Description of 'comp2'}
%% ...
}
\references{
%% ~put references to the literature/web site here ~
}
\author{
%% ~~who you are~~
}
\note{
%% ~~further notes~~
}
%% ~Make other sections like Warning with \section{Warning }{....} ~
\seealso{
%% ~~objects to See Also as \code{\link{help}}, ~~~
}
\examples{
##---- Should be DIRECTLY executable !! ----
##-- ==> Define data, use random,
##-- or do help(data=index) for the standard data sets.
## The function is currently defined as
function (object, method = "clustering", ...)
{
if (method == "clustering") {
clust.genes <- amap::hcluster(x = object@qData, method = "pearson",
link = "average")
clust.arrays <- amap::hcluster(x = t(object@qData), method = "pearson",
link = "average")
heatcol <- colorRampPalette(c("Green", "Red"))(32)
heatmap(x = object@qData, Rowv = as.dendrogram(clust.genes),
Colv = as.dendrogram(clust.arrays), col = heatcol)
}
if (method == "pca") {
p <- amap::acp(object@qData)
plot(p)
return(p)
}
if (method == "kmeans") {
cl <- amap::Kmeans(object@qData, 2)
plot(object@qData, col = cl$cluster)
points(cl$centers, col = 1:2, pch = 8, cex = 2)
return(cl)
}
}
}
% Add one or more standard keywords, see file 'KEYWORDS' in the
% R documentation directory.
\keyword{ ~kwd1 }
\keyword{ ~kwd2 }% __ONLY ONE__ keyword per line