MCPAN (version 1.1-21)

poly3estf: Only for internal use.

Description

Poly-3- adjusted point and variance estimates for long term carcinogenicity data if data are given as a numeric time vector, a logical status vector and a factor containing a grouping variable

Usage

poly3estf(time, status, tmax=NULL, f, method = "BP", k=NULL)

Arguments

time

a numeric vector of times of death of the individuals

status

a logical (or numeric, consisting of 0,1 only) vector giving the tumour status at time of death of each individual, where TRUE (1) = tumour present, FALSE (0) = no tumour present

tmax

a single numeric value, the time of sacrifice in the trial, or the last last time of death, defaults to the maximal value observed in time

f

a factor of the same length as time, status, giving the levels of a grouping variable in a one-way layout

method

a single charcter string, specifying the method for adjustment, with options: "BP" (Bailer Portier: assuming poly-3-adjusted rates are binomial variables), "BW" (Bieler, Williams: delta method as in Bieler-Williams (1993)) "ADD1" (as Bailer Portier, including an add1-adjustment on the raw tumour rates) "ADD2" (as Bailer Portier, including an add2-adjustment on the raw tumour rates following Agresti Caffo (2000) for binomials)

k

a single numeric value, the exponent to calculate survival adjusted proportions according to Bailer and Portier (1988), defaults to 3

Value

A list containing:

Y

a numeric vector, groupwise number of tumours

n

a numeric vector, groupwise number of individuals

estimate

a numeric vector, groupwise poly-3-adjusted rates according to Bailer, Portier (1988)

weight

a numeric vector of poly-3-adjusted weights

estp

a numeric vector, groupwise poly-3-adjusted rate (according to method)

nadj

adjusted n (sum of weights)

varp

a numeric vector, groupwise variance estimate (according to method)

varcor

a numeric vector, groupwise variance estimate, if necessary corrected such that estimates of 0 can not occur

names

a character vector, the levels of the grouping variable f

k

a single numeric value, as input

Details

For internal use.

References

See poly3est

Examples

# NOT RUN {
data(bronch)

poly3estf(status=bronch$Y, time=bronch$time, f=bronch$group, k=3)

poly3estf(status=bronch$Y, time=bronch$time, f=bronch$group, k=5)

# }