permute.BACI {emon}

Does non-parametric randomisation test for the interaction term in a BACI design.

Description

We have control and treatment data from time 1 in a BACI design, plus control and treatment data from time 2. The interaction the amount that the difference in the control and treatment meansis different between times 1 and 2.

Usage

permute.BACI(t1, c1, t2, c2, nreps=999)

Arguments

t1

Data vector for the treatment at time 1

c1

Data vector for the control at time 1

t2

Data vector for the treatment at time 2

c2

Data vector for the control at time 2

nreps

Number of replications used in the randomisation and generation of the p-value. Default is nreps=999

Details

There are several permutation that can be used to generate the null distribution for the interaction (see Manly, 2006 and Anderson and Terr Braak, 2003). The method used here is to do a complete randomisation of the raw data.

The p-value is calculated as suggested by Manly (2006).

Value

The p-value is returned as $p.value

Author(s)

Jon Barry: [e-mail]

References

Manly BFJ (2006) Randomization, Bootstrap And Monte Carlo Methods in Biology: 3rd edition. Chapman and Hall.

Anderson, M.J. and Ter Braak, C.J.F. (2003). Permutation tests for multi-factorial analysis of variance. Journal of Computation and Simulation, 73, 85-113.

See Also

power.BACI, permute.groups