git-svn-id: file:///home/git/hedgehog.fhcrc.org/bioconductor/trunk/madman/Rpacks/MineICA@73179 bc3139a8-67e5-0310-9ffc-ced21a209358
Marc Carlson authored on 05/02/2013 22:15:53
| 1 | 1 |
new file mode 100644 |
| ... | ... |
@@ -0,0 +1,54 @@ |
| 1 |
+\name{plotMix}
|
|
| 2 |
+\alias{plotMix}
|
|
| 3 |
+\title{Plots an histogram and Gaussian fitted by \code{\link[mclust]{Mclust}}}
|
|
| 4 |
+\usage{
|
|
| 5 |
+ plotMix(mc, data, nbBreaks, traceDensity = TRUE, |
|
| 6 |
+ title = "", xlim, ylim, ...) |
|
| 7 |
+} |
|
| 8 |
+\arguments{
|
|
| 9 |
+ \item{mc}{The result of Mclust function applied to
|
|
| 10 |
+ argument \code{data}}
|
|
| 11 |
+ |
|
| 12 |
+ \item{data}{A vector of numeric values}
|
|
| 13 |
+ |
|
| 14 |
+ \item{nbBreaks}{The number of breaks for the histogram}
|
|
| 15 |
+ |
|
| 16 |
+ \item{traceDensity}{If TRUE (default) density are
|
|
| 17 |
+ displayed on the y-axis, else if FALSE counts are |
|
| 18 |
+ displayed on the y-acis} |
|
| 19 |
+ |
|
| 20 |
+ \item{title}{A title for the plot}
|
|
| 21 |
+ |
|
| 22 |
+ \item{xlim}{x-axis limits to be used in the plot}
|
|
| 23 |
+ |
|
| 24 |
+ \item{ylim}{y-axis limits to be used in the plot}
|
|
| 25 |
+ |
|
| 26 |
+ \item{...}{additional arguments for hist}
|
|
| 27 |
+} |
|
| 28 |
+\value{
|
|
| 29 |
+ NULL |
|
| 30 |
+} |
|
| 31 |
+\description{
|
|
| 32 |
+ Given a result of function \code{Mclust} applied to a
|
|
| 33 |
+ numeric vector, this function draws the fitted Gaussian |
|
| 34 |
+ on the histogram of the data values. |
|
| 35 |
+} |
|
| 36 |
+\details{
|
|
| 37 |
+ A shapiro test p-value is added to the plot title. This |
|
| 38 |
+ function can only deal with at the most three Gaussian. |
|
| 39 |
+} |
|
| 40 |
+\examples{
|
|
| 41 |
+## create a mix of two Gaussian |
|
| 42 |
+v <-c(rnorm(80,mean=-0.5,sd=1),rnorm(80,mean=1,sd=0.2)) |
|
| 43 |
+## apply Mclust |
|
| 44 |
+mc <- Mclust(v) |
|
| 45 |
+## plot fitted Gaussian on histogram of v |
|
| 46 |
+plotMix(mc=mc,data=v,nbBreaks=30) |
|
| 47 |
+} |
|
| 48 |
+\author{
|
|
| 49 |
+ Anne Biton |
|
| 50 |
+} |
|
| 51 |
+\seealso{
|
|
| 52 |
+ \code{\link{hist}}, \code{\link[mclust]{Mclust}}
|
|
| 53 |
+} |
|
| 54 |
+ |