Last data update: 2014.03.03

R: Permutation t-test
perm.t.testR Documentation

Permutation t-test

Description

Two Sample t-test via monte-carlo permutation

Usage

perm.t.test(x,y,statistic=c("t","mean"),
			alternative=c("two.sided", "less", "greater"), midp=TRUE, B=10000)

Arguments

x

a numeric vector containing the first sample

y

a numeric vector containing the second sample

statistic

The statistic to be permuted. See details

alternative

The alternative hypothesis

midp

should the mid p-value be used

B

The number of monte-carlo samples to be generated

Details

This function performs a two sample permutation test. If the mean is permuted, then the test assumes exchangability between the two samples. if the t-statistic is used, the test assumes either exchangability or a sufficiently large sample size. Because there is little lost in the way of power, and the assumptions are weaker, the t-statistic is used by default.

Value

A list with class "htest" containing the following components:

statistic

The observed value of the statistic.

p.value

the p-value for the test.

method

a character string indicating the type of test performed.

data.name

a character string giving the name(s) of the data.

B

The number of samples generated

alternative

the direction of the test

See Also

t.test

Examples

perm.t.test(rnorm(100),runif(100,-.5,.5))

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(Deducer)
Loading required package: ggplot2
Loading required package: JGR
Loading required package: rJava
Loading required package: JavaGD
Loading required package: iplots

Please type JGR() to launch console. Platform specific launchers (.exe and .app) can also be obtained at http://www.rforge.net/JGR/files/.


Loading required package: car
Loading required package: MASS


Note Non-JGR console detected:
	Deducer is best used from within JGR (http://jgr.markushelbig.org/).
	To Bring up GUI dialogs, type deducer().

> png(filename="/home/ddbj/snapshot/RGM3/R_CC/result/Deducer/perm.t.test.Rd_%03d_medium.png", width=480, height=480)
> ### Name: perm.t.test
> ### Title: Permutation t-test
> ### Aliases: perm.t.test
> 
> ### ** Examples
> 
> perm.t.test(rnorm(100),runif(100,-.5,.5))

	Two-Sample permutation test using Welsh's t (B= 10000 )

data:  rnorm(100) and runif(100, -0.5, 0.5)
Welsh t-statistic = -0.32383, p-value = 0.7397
alternative hypothesis: two.sided

> 
> 
> 
> 
> 
> dev.off()
null device 
          1 
>