Compare activity level estimates

Share:

Description

Wald test for the statistical difference between two or more activitiy level estimates.

Usage

1
compareAct(fits)

Arguments

fits

A list of fitted actmod objects

Details

Uses a Wald test to ask whether the difference between estimates a1 and a2 is significantly different from 0: statistic W = (a1-a2)^2 / (SE1^2+SE2^2) tested on chi-sq distribution with 1 degree of freedom.

Value

A matrix with 4 columns: 1. differences between estimates; 2. SEs of the differences; 3. Wald statistics; 4. p-values (H0 is no difference between estimates). Matrix rows give all possible pairwise comparisons, numbered in the order in which they entered in the list fits.

Examples

1
2
3
4
5
6
7
8
#Test whether paca have a sigificantly different activity level from rat.
#Bootstrap reps limited to speed up example.
data(BCItime)
tPaca <- 2*pi*BCItime$time[BCItime$species=="paca"]
tRat <- 2*pi*BCItime$time[BCItime$species=="rat"]
(fPaca <- fitact(tPaca, reps=10))
(fRat <- fitact(tRat, reps=10))
compareAct(list(fPaca,fRat))