R2ucare (version 1.0.2)

test2ct: Test2.CT

Description

This function performs Test2.CT

Usage

test2ct(X, freq, verbose = TRUE, rounding = 3)

Value

This function returns a list with first component the overall test and second component a data.frame with 5 columns for components i (2:K-2) (in rows) of test2.Cti: component, degree of freedom, statistic of the test, p-value, signed test, test performed.

Arguments

X

is a matrix of encounter histories with K occasions

freq

is a vector of the number of individuals with the corresponding encounter history

verbose

controls the level of the details in the outputs; default is TRUE for all details

rounding

is the level of rounding for outputs; default is 3

Author

Olivier Gimenez <[e-mail]>, Jean-Dominique Lebreton, Rémi Choquet, Roger Pradel

Examples

# read in the classical dipper dataset
dipper = system.file("extdata", "ed.inp", package = "R2ucare")
dipper = read_inp(dipper,group.df=data.frame(sex=c('Male','Female')))

# Get encounter histories, counts and groups:
dip.hist = dipper$encounter_histories
dip.freq = dipper$sample_size
dip.group = dipper$groups

# split the dataset in males/females
mask = (dip.group == 'Female')
dip.fem.hist = dip.hist[mask,]
dip.fem.freq = dip.freq[mask]
mask = (dip.group == 'Male')
dip.mal.hist = dip.hist[mask,]
dip.mal.freq = dip.freq[mask]

# for females
X = dip.fem.hist
freq = dip.fem.freq
res.females = test2ct(X,freq)
res.females