R: Plots the results from a permutation based version of...
plot.hotelling.test
R Documentation
Plots the results from a permutation based version of Hotelling's
T-squared test for the difference in two multivariate sample means
Description
Plots a histogram of the distribution of the permuted test statistics
for a permutation version of Hotelling's T-squared
Usage
## S3 method for class 'hotelling.test'
plot(x,...)
Arguments
x
an object of type hotelling.test
...
any additional arguments to be passed to the hist command
Details
This function only works if you have performed a permutation test. It
will return an error message if not. It could be programmed to draw
the relevant F distribution in the standard case, but this seems
rather pointless.
Author(s)
James M. Curran
Examples
data(bottle.df)
bottle.df = subset(bottle.df, Number == 1)
bottle.df$Number = rep(1:2,c(10,10))
fit = hotelling.test(.~Number, bottle.df, perm = TRUE)
plot(fit)
plot(fit, col = "lightblue")
Results
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(Hotelling)
Loading required package: corpcor
> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Hotelling/plot.hotelling.test.Rd_%03d_medium.png", width=480, height=480)
> ### Name: plot.hotelling.test
> ### Title: Plots the results from a permutation based version of
> ### Hotelling's T-squared test for the difference in two multivariate
> ### sample means
> ### Aliases: plot.hotelling.test
> ### Keywords: plot
>
> ### ** Examples
>
> data(bottle.df)
> bottle.df = subset(bottle.df, Number == 1)
> bottle.df$Number = rep(1:2,c(10,10))
> fit = hotelling.test(.~Number, bottle.df, perm = TRUE)
| | | 0% | |= | 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% | |============================================================ | 85% | |============================================================ | 86% | |============================================================= | 87% | |============================================================== | 88% | |============================================================== | 89% | |=============================================================== | 90% | |================================================================ | 91% | |================================================================ | 92% | |================================================================= | 93% | |================================================================== | 94% | |================================================================== | 95% | |=================================================================== | 96% | |==================================================================== | 97% | |===================================================================== | 98% | |===================================================================== | 99% | |======================================================================| 100%> plot(fit)
> plot(fit, col = "lightblue")
>
>
>
>
>
> dev.off()
null device
1
>