Nothing
Pplot <-
function( rates,
age = as.numeric( dimnames( rates )[[1]] ),
per = as.numeric( dimnames( rates )[[2]] ),
grid = FALSE,
p.grid = grid,
ygrid = grid,
col.grid = gray( 0.9 ),
p.lim = range( per, na.rm=TRUE ) + c(0,diff(range(per))/30),
ylim = range( rates[rates>0], na.rm=TRUE ),
p.lab = names( dimnames( rates ) )[2],
ylab = deparse( substitute( rates ) ),
at = NULL,
labels = paste( at ),
type = "l",
lwd = 2,
lty = 1,
col = par( "fg" ),
log.ax = "y",
las = 1,
ann = FALSE,
cex.ann = 0.8,
xannx = 1/20,
a.thin = seq( 1, length( age ), 2 ),
...
)
{
# Plot the frame
if( ann ) p.lim <- p.lim + c(0,diff( range( age ) ) * xannx )
matplot( per, t(rates), type="n",
xlim=p.lim, ylim=ylim, xlab=p.lab, ylab=ylab,
log=log.ax, las=las, yaxt=if( !is.null( at ) ) "n" else "s" )
if( !is.null( at ) ) axis( side=2, at=at, labels=labels, yaxt="s", las=las )
# and the grid if required
if( !missing( p.grid ) | !missing( grid ) )
{
if( is.logical( p.grid ) & p.grid[1] ) p.grid <- nice( per, log=par("xlog") )
abline( v=p.grid, col=col.grid )
}
if( !missing( ygrid ) | !missing( ygrid ) )
{
if( is.logical( ygrid ) & ygrid[1] ) ygrid <- nice( rates[!is.na(rates)], log=par("ylog") )
abline( h=ygrid, col=col.grid )
}
box()
# then the curves
matlines( per, t(rates), lwd=lwd, lty=lty, col=col, type=type, ... )
# annotate them if required (every second by default )
if( ann )
{
nr <- nrow( rates )
nc <- ncol( rates )
text( rep( per[nc], nr )[a.thin], rates[,nc][a.thin],
paste( "", age[a.thin] ), adj=c(0,0.5),
cex=cex.ann, col=if( length(col)==1 ) col else col[a.thin] )
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.