Visually presents and compares different multivariate location estimates and their confidence ellipsoids.
Usage
## S3 method for class 'mvloc'
plot(x, est2 = NULL, est3 = NULL, X = NULL, ...)
Arguments
x
an object of class mvl1lm.
est2
An optional additional location estimate. A list with the components location, vcov and est.name, for example an object of class 'mvloc'.
est3
An optional additional location estimate. A list with the components location, vcov and est.name, for example an object of class 'mvloc'.
X
a numeric data frame or matrix. Optional data points on which the estimates could have been based.
...
optional plotting arguments. For details see plotMvloc .
Details
The figure can be used to compare different multivariate location estimates. The location of the legend is currently problematic and it is
recommended that the user should provide the coordinates for the legend. The function calls plotMvloc.
R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
Copyright (C) 2016 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)
R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.
R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.
Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.
> library(MNM)
Loading required package: ICSNP
Loading required package: mvtnorm
Loading required package: ICS
Loading required package: SpatialNP
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/MNM/plot.mvloc.Rd_%03d_medium.png", width=480, height=480)
> ### Name: plot.mvloc
> ### Title: Plotting Method for mvloc Objects
> ### Aliases: plot.mvloc
> ### Keywords: methods
>
> ### ** Examples
>
> X <- rmvt(50, diag(c(1, 2)), 3)
>
> est1 <- mv.1sample.est(X)
> est2 <- mv.1sample.est(X, "sign")
> est3 <- mv.1sample.est(X, "rank", "inner")
>
> plot(est1)
> plot(est1, est2, est3, X, alim="b", lty.ell=1:3, pch.ell=14:16)
> plot(est1, est2, est3, X, alim="e")
>
>
>
>
>
> dev.off()
null device
1
>