I don't know if this is useful of not.
But here's some dis-aggregated US data.
It's about 10 days old, and I need to improve the dis-aggregation algorithm.
------------------------------------------------------
age <- 0:80
rel.mort.rate <- c (0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5.97119198173175e-14, 1.49279799543293e-12,
1.80329997848298e-11, 1.40323011570696e-10, 7.91302361419091e-10,
3.45021443896442e-09, 0.0000000121231319447505, 0.0000000353442018829079,
0.0000000874198628343878, 0.000000186810950489511, 0.000000350496608960251,
0.000000586292710567257, 0.00000088809586113832, 0.00000123828572996647,
0.00000161636339411418, 0.00000200770373248443, 0.00000240656741448015,
0.00000281336507897391, 0.00000323131179486445, 0.00000366643530132155,
0.00000412942874923419, 0.00000463520698193433, 0.00000519914074326175,
0.00000583299029582447, 0.00000654330938109396, 0.00000733211405261902,
0.00000819881232642656, 0.00000914397689953956, 0.0000101752369848173,
0.0000113121562971092, 0.0000125849353938131, 0.0000140250299768815,
0.0000156519020717643, 0.0000174632527695486, 0.000019434044086,
0.0000215255604245307, 0.0000237027683240887, 0.0000259554459399679,
0.0000283147468915187, 0.0000308551986409676, 0.0000336771894650306,
0.0000368753737436115, 0.0000405063605263262, 0.0000445688466005247,
0.0000490032796563514, 0.0000537118848883967, 0.0000585949512128272,
0.0000635935967055111, 0.0000687239859205501, 0.0000740884186923496,
0.0000798570919656766, 0.0000862259973389807, 0.0000933650267485426,
0.000101373734773254, 0.000110260898647345, 0.0001199575254273,
0.000130361177821875, 0.000141396704797265, 0.000153071840957868,
0.00016550903061107, 0.000178943792317475, 0.000193689554660746,
0.000210077052171998, 0.00022838282771009, 0.000248763951569319,
0.000271212123374848, 0.000295531695409173, 0.000321339382531619,
0.000348084036019435, 0.000375091815336575, 0.000401647747211891)
names (rel.mort.rate) <- age
options (scipen=4)
mar <- par ("mar")
mar [2] <- mar [2] + 1
p0 <- par (mar=mar)
barplot (rel.mort.rate [41:81], ylim = c (0, 0.0004), xlab="age",
ylab="relative mortality rate\n(ndeaths / npeople | age)")
par (p0)
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
