Description Usage Arguments Details Note Author(s) References See Also Examples
View source: R/apc_plot_data.R
Produces plots showing age, period and cohort sums. As a default this is done both for responses and dose, giving a total of six plots.
1 2 3 | apc.plot.data.sums(apc.data.list,data.type="a",
average=FALSE,keep.incomplete=TRUE,apc.index=NULL,
type="o",log="",main.outer=NULL,main.sub=NULL)
|
apc.data.list |
List. See |
data.type |
Optional. Character. "r","d","m","a" if sums are computed for responses, dose, (mortality rates), all. Rates are computed as responses/doses. Default is "a". |
average |
Optional. Logical. Sums are reported if FALSE, Averages are reported if TRUE. Default is FALSE. |
keep.incomplete |
Optional. Logical. If true perform calculation for incomplete sequences by removing NA.
If false incomplete sequences are NA. See example in |
apc.index |
Optional. List. See |
type |
Optional |
log |
Optional |
main.outer |
Optional. Character. Main title for plot, to be shown in outer margin. Default is NULL, in which case a title is generated internally. |
main.sub |
Optional. Titles for sub plots. Use with data.type "r","d","m". For data.type "a" use default. Default is NULL, in which case a title is generated internally. |
The data sums are computed using apc.data.sums
. Then plotted as requested.
Use apc.data.sums
if numerical values needed.
Bent Nielsen <bent.nielsen@nuffield.ox.ac.uk> 15 Aug 2018 (15 Dec 2013)
Martinez Miranda, M.D., Nielsen, B. and Nielsen, J.P. (2015) Inference and forecasting in the age-period-cohort model with unknown exposure with an application to mesothelioma mortality. Journal of the Royal Statistical Society A 178, 29-55. Download: Article, Nuffield DP.
The example below uses Japanese breast cancer data, see data.Japanese.breast.cancer
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 | #####################
# EXAMPLE with artificial data
# Generate a 3x4 matrix in "AP" data.format with the numbers 1..12
# Then make a data list
# Then plot data sums.
# Note only 3 plots are made as there are no doses
m.data <- matrix(data=seq(length.out=12),nrow=3,ncol=4)
m.data
data.list <- apc.data.list(m.data,"AP")
apc.plot.data.sums(data.list)
apc.plot.data.sums(data.list,average=TRUE)
apc.plot.data.sums(data.list,keep.incomplete=FALSE)
#####################
# EXAMPLE with Japanese breast cancer data
# Make a data list
# Then plot data sums for both responses and doses.
data.list <- data.Japanese.breast.cancer()
apc.plot.data.sums(data.list)
# Or plot data sums for responses only
apc.plot.data.sums(data.list,data.type="r")
#####################
# EXAMPLE with asbestos data
# Miranda Martinex, Nielsen and Nielsen (2013).
# This is Figure 1,a-c
data.list <- data.asbestos()
apc.plot.data.sums(data.list,type="l")
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.