Hi everyone I want help in graduating the attached rates and checking for goodness of fit and smoothness using R please help. Many thnk TOo every one around the world This message and attachments are subject to a disclaimer. Please refer to www.it.up.ac.za/documentation/governance/disclaimer/ for full details. / Hierdie boodskap en aanhangsels is aan 'n vrywaringsklousule onderhewig. Volledige besonderhede is by www.it.up.ac.za/documentation/governance/disclaimer/ beskikbaar. -------------- next part -------------- Age Males females 55 0.00000000 0.00000000 56 0.00000000 0.00000000 57 0.00888883 0.00000000 58 0.00000000 0.00803209 59 0.00000000 0.00000000 60 0.00901373 0.00296296 61 0.00684929 0.00000000 62 0.00770214 0.00182648 63 0.00386847 0.00165426 64 0.00476568 0.00141543 65 0.00298507 0.00127796 66 0.00229885 0.00000000 67 0.00203046 0.00210084 68 0.00448229 0.00000000 69 0.00602667 0.00103252 70 0.00092208 0.00000000 71 0.00078895 0.00105988 72 0.00537152 0.00097040 73 0.00405268 0.00219539 74 0.00434216 0.00212992 75 0.00296004 0.00117855 76 0.00707182 0.00125235 77 0.00316957 0.00126984 78 0.00665924 0.00284091 79 0.01002777 0.00404040 80 0.01138508 0.00460475 81 0.01070940 0.00162470 82 0.00995017 0.00372439 83 0.00715560 0.00198610 84 0.00701751 0.00703397 85 0.02061782 0.00232829 86 0.01508267 0.00282087 87 0.01733058 0.00000000 88 0.02020133 0.00389863 89 0.02227078 0.00000000
On Sep 23, 2009, at 4:54 AM, MKHABELA,SN wrote:> Hi everyone > > I want help in graduating the attached rates and checking for > goodness of fit and smoothness using R please help. > > <Mortality rates for males and > females.txt>______________________________________________You have provided the rates but not the death counts or person-years of exposure that would be needed to create these rates. Such numbers would be useful in doing a graduation that had a statistical basis. At the moment I think you should be looking at the data and deciding whether you can model as a Gompertz function which is the canonical "mortality law": mrates <- read.table(file="<path>/Mortality rates for males and females.txt", header=T) with(mrates, plot(Age, log(Males+.001))) # added 1/2 the minimum rate to the zero entries with(mrates, plot(Age, log(females+.0005))) But I do not think you have given enough information or background to take further steps. The plots are not entirely supportive of a Gompertz model and is "telling" you that you need to think more deeply about what the underlying situation really represents. Of course, if this is homework, you have not been completely forthright with us and should review your educational institution's regulations. -- David Winsemius, MD Heritage Laboratories West Hartford, CT