Epi (version 2.50)

plot.apc: Plot the estimates from a fitted Age-Period-Cohort model

Description

This function plots the estimates created by apc.fit in a single graph. It just calls apc.frame after computing some sensible values of the parameters, and subsequently plots the estimates using apc.lines.

Usage

# S3 method for apc
plot( x, r.txt="Rate", ...)
          apc.plot( x, r.txt="Rate", ...)

Value

A numerical vector of length two, with names

c("cp.offset","RR.fac"). The first is the offset for the cohort period-axis, the second the multiplication factor for the rate-ratio scale. Therefore, if you want to plot at (x,y) in the right panel, use (x-res["cp.offset"],y/res["RR.fac"])

=(x-res[1],y/res[2]). This vector should be supplied for the parameter frame.par to

apc.lines if more sets of estimates is plotted in the same graph, however see cp.points.

Arguments

x

An object of class apc.

r.txt

The text to put on the vertical rate axis.

...

Additional arguments passed on to apc.lines.

Author

Bendix Carstensen, Steno Diabetes Center, http://bendixcarstensen.com

Details

plot.apc is just a wrapper for apc.plot.

See Also

apc.lines, lines.apc, apc.frame, apc.fit

Examples

data( lungDK )
apc1 <- apc.fit( transform( lungDK,
                            A = Ax, P = Px, Y = Y/10^5 ),
                 ref.c = 1920 )
fp <- apc.plot( apc1 )
apc.lines( apc1, frame.par=fp, drift=1.01, col="red" )
for( i in 1:11 )
  apc.lines( apc1, frame.par=fp, drift=1+(i-6)/100, col=rainbow(12)[i] )