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
|
\name{s.multinom}
\alias{s.multinom}
\title{Graph of frequency profiles (useful for instance in genetic)}
\description{
The main purpose of this function is to draw categories using scores and profiles by their gravity center.
Confidence intervals of the average position (issued from a multinomial distribution) can be superimposed.
}
\usage{
s.multinom(dfxy, dfrowprof, translate = FALSE, xax = 1, yax = 2,
labelcat = row.names(dfxy), clabelcat = 1, cpointcat = if (clabelcat == 0) 2 else 0,
labelrowprof = row.names(dfrowprof), clabelrowprof = 0.75,
cpointrowprof = if (clabelrowprof == 0) 2 else 0, pchrowprof = 20,
coulrowprof = grey(0.8), proba = 0.95, n.sample = apply(dfrowprof, 1, sum),
axesell = TRUE, ...)
}
\arguments{
\item{dfxy}{\code{dfxy} is a data frame containing at least two numerical variables.
The rows of \code{dfxy} are categories such as 1,2 and 3 in the triangular plot.}
\item{dfrowprof}{\code{dfrowprof} is a data frame whose the columns are the rows of \code{dfxy}.
The rows of \code{dfxy} are profiles or frequency distributions on the categories.
The column number of \code{dfrowprof} must be equal to the row number of \code{dfxy}.
\code{row.names(dfxy)} and \code{names(dfrowprof)} must be identical. }
\item{translate}{a logical value indicating whether the plot should be translated(TRUE) or not.
The origin becomes the gravity center weighted by profiles. }
\item{xax}{the column number of \code{dfxy} for the x-axis }
\item{yax}{the column number of \code{dfxy} for the y-axis }
\item{labelcat}{a vector of strings of characters for the labels of categories }
\item{clabelcat}{an integer specifying the character size for the labels of categories,
used with \code{par("cex")*clabelcat} }
\item{cpointcat}{an integer specifying the character size for the points showing the categories,
used with \code{par("cex")*cpointcat} }
\item{labelrowprof}{a vector of strings of characters for the labels of profiles (rows of \code{dfrowprof}) }
\item{clabelrowprof}{an integer specifying the character size for the labels of profiles
used with par("cex")*clabelrowprof}
\item{cpointrowprof}{an integer specifying the character size for the points representative
of the profiles used with par("cex")*cpointrowprof }
\item{pchrowprof}{either an integer specifying a symbol or a single character to be used for the profile labels }
\item{coulrowprof}{a vector of colors used for ellipses, possibly recycled}
\item{proba}{a value lying between 0.500 and 0.999 to draw a confidence interval }
\item{n.sample}{a vector containing the sample size, possibly recycled. Used \code{n.sample = 0} if the profiles
are not issued from a multinomial distribution and that confidence intervals have no sense. }
\item{axesell}{a logical value indicating whether the ellipse axes should be drawn}
\item{\dots}{further arguments passed from the \code{s.label} for the initial scatter plot. }
}
\value{
Returns in a hidden way a list of three components :
\item{tra}{a vector with two values giving the done original translation. }
\item{ell}{a matrix, with 5 columns and for rows the number of profiles, giving the means,
the variances and the covariance of the profile for the used
numerical codes (column of \code{dfxy})}
\item{call}{the matched call}
}
\author{Daniel Chessel
}
\examples{
par(mfrow = c(2,2))
par(mar = c(0.1,0.1,0.1,0.1))
proba <- matrix(c(0.49,0.47,0.04,0.4,0.3,0.3,0.05,0.05,0.9,0.05,0.7,0.25), ncol = 3, byrow = TRUE)
proba.df <- as.data.frame (proba)
names(proba.df) <- c("A","B","C") ; row.names(proba.df) <- c("P1","P2","P3","P4")
w.proba <- triangle.plot(proba.df, clab = 2, show = FALSE)
box()
w.tri = data.frame(x = c(-sqrt(1/2),sqrt(1/2),0), y = c(-1/sqrt(6),-1/sqrt(6),2/sqrt(6)))
L3 <- c("A","B","C")
row.names(w.tri) <- L3
s.multinom(w.tri, proba.df, n.sample = 0, coulrowprof = "black", clabelrowprof = 1.5)
s.multinom(w.tri, proba.df, n.sample = 30, coul = palette()[5])
s.multinom(w.tri, proba.df, n.sample = 60, coul = palette()[6], add.p = TRUE)
s.multinom(w.tri, proba.df, n.sample = 120, coul = grey(0.8), add.p = TRUE)
print(s.multinom(w.tri, proba.df[-3,], n.sample = 0, translate = TRUE)$tra)
}
\keyword{multivariate}
\keyword{hplot}
|
