nordpred.getpred: Gets the observed and predicted incidence rates on matrix...

In haraldwf/nordpred: Prediction of cancer incidence (as used in the Nordpred project)

Description

nordpred.getpred uses a nordpred object (see nordpred.object) to extract the observed and predicted incidence rates

Usage

nordpred.getpred(nordpred.object, incidence = TRUE, standpop = NULL,
  excludeobs = FALSE, byage, agegroups = "all")

Arguments

nordpred.object	An object based on nordpred or nordpred.prediction.
incidence	Indicates whether to give incidence or number of cases
standpop	A vector of weights for age standardisation. Default is no standardisation (crude rates), but using a standardisation (for the suitable no of age groups) is recommended
excludeobs	Exclude number for observed periods and only give numbers for predicted periods
byage	Report numbers by age groups. If false, crude or age standardised rates are given
agegroups	Which agegroups to include. E.g. c(5:18) includes age groups five to eighteen

Value

an object of class nordpred (see nordpred.object).

Note for S-plus

Powerlink is made via a special modification in S-PLUS. This works fine for the point estimates, but the variance estimates found via the glm-objects are wrong. For variance estimates, we would rather recommend using R.

Author(s)

Harald Fekjaer and Bjorn Moller (Cancer Registry of Norway)

References

• A website for nordpred is available at: http://www.kreftregisteret.no/software/nordpred/

• Background for the methods can be found in: Moller B., Fekjaer H., Hakulinen T., Sigvaldason H, Storm H. H., Talback M. and Haldorsen T 'Prediction of cancer incidence in the Nordic countries: Empirical comparison of different approaches' Statistics in Medicine 2003; 22:2751-2766

• An application of the function, using all the default settings, can be found in: Moller B, Fekjaer H, Hakulinen T, Tryggvadottir L, Storm HH, Talback M, Haldorsen T. Prediction of cancer incidence in the Nordic countries up to the year 2020. Eur J Cancer Prev Suppl 2002; 11: S1-S96

See Also

Other nordpred: is.nordpred.estimate, nordpred.estimate; is.nordpred, nordpred; nordpred.estimate.object; nordpred.object; nordpred.prediction; plot.nordpred; print.nordpred.estimate; print.nordpred; summary.nordpred

Examples

# data (Colon cancer for Norwegian males)
indata

# Create dataset with observed and predicted population
inpop <- cbind(inpop1, inpop2)

# Fit model & predict new incidence:
res <- nordpred(indata, inpop, startestage = 5, startuseage = 6, cuttrend = c(0,.25,.5,.75,.75))
res2 <- nordpred(indata, inpop, startestage = 5, startuseage = 6,
         cuttrend = c(0, .25, .5, .75, .75), linkfunc = 'poisson')

# Print / get results:
print(res)
nordpred.getpred(res)
summary(res, printpred = FALSE)

# Get results with standardisation:'
wstand <- c(0.12, 0.1, 0.09, 0.09, 0.08, 0.08, 0.06, 0.06, 0.06, 0.06,0.05,
            0.04, 0.04, 0.03, 0.02, 0.01, 0.005, 0.005)
round(nordpred.getpred(res, incidence = TRUE, standpop = NULL), 2)
round(nordpred.getpred(res, incidence = TRUE, standpop = wstand), 2)

haraldwf/nordpred documentation built on May 17, 2019, 2:41 p.m.