Scribd
Upload
23 views5 pages

SMG Makeham's Law

ok

Uploaded by

Nandini Rachcha
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
23 views5 pages

SMG Makeham's Law

ok

Uploaded by

Nandini Rachcha
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as DOCX, PDF, TXT or read online on Scribd
You are on page 1/ 5

Name : Nandini Rachcha Practical No.

: 2

Title of Practical : Construction of life table using Makeham’s Mortality Law Date: 06/08/2024

> #Q1)

> lb=seq(40,70,10);lb

[1] 40 50 60 70

> ub=seq(49,79,10);ub

[1] 49 59 69 79

> age=(lb+ub)/2

> births = c(100,50,20,10)

> deaths = c(10,20,30,40)

> # Constants for Makeham's law (These values would need to be determined or estimated)

> A = 0.0001

> B = 0.01

> C = 0.1

> # Calculated mortality rates using Makeham's Law

> mortality_rate = A + B * exp(C * age)

> # These calculations are similar to those in traditional life table construction

> # Display the life table

> lifeotable = data.frame(Age = age, Birth = births, Deaths = deaths, Mortality_Rate = mortality_rate)

> print(lifeotable)

Age Birth Deaths Mortality_Rate

1 44.5 100 10 0.8563694

2 54.5 50 20 2.3276817

3 64.5 20 30 6.3271229

4 74.5 10 40 17.1987315

> plot(age,mortality_rate,"b",main="")
> #Q2)

> age = c(0, 1, 2, 3, 4, 5, 6)

> deaths = c(10, 8, 7, 6, 5, 4, 3)

> person_lived = c(10000, 9992, 9985, 9978, 9971, 9964, 9957)

> # Constants foe Makeham's Law

> A = 0.0001

> B = 0.01

> C = 0.1

> # Calculated mortality rates using Makeham's Law

> mortality_rate = A + B * exp(C * age)

> # Display the life table


> lifeotable = data.frame(Age = age, Person_Lived = person_lived, Deaths = deaths, Mortality_Rate =
mortality_rate)

> print(lifeotable)

Age Person_Lived Deaths Mortality_Rate

1 0 10000 10 0.01010000

2 1 9992 8 0.01115171

3 2 9985 7 0.01231403

4 3 9978 6 0.01359859

5 4 9971 5 0.01501825

6 5 9964 4 0.01658721

7 6 9957 3 0.01832119

> plot(age,mortality_rate,"b",main="")

>
> #Q3)

> age3 = 0:95

> deaths = c(10, 12, 15, 18, 21, 24, 28, 33, 39, 45,

+ 52, 60, 69, 79, 90, 102, 115, 129, 144, 160,

+ 177, 195, 214, 234, 255, 277, 300, 324, 349, 375,

+ 402, 430, 459, 489, 520, 552, 585, 619, 654, 690,

+ 727, 765, 804, 844, 885, 927, 970, 1014, 1059, 1105,

+ 1152, 1200, 1249, 1299, 1350, 1402, 1455, 1509, 1564,

+ 1620, 1677, 1735, 1794, 1854, 1915, 1977, 2040, 2104,

+ 2169, 2235, 2302, 2370, 2439, 2509, 2580, 2652, 2725,

+ 2799, 2874, 2950, 3027, 3105, 3184, 3264, 3345, 3427,

+ 3510, 3594, 3679, 3765, 3852, 3940, 4029, 4119, 4210)

> length(deaths)

[1] 95

> person_years_lived

[1] 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000

[13] 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000

[25] 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000

[37] 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000

[49] 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000

[61] 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000

[73] 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000

[85] 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000 10000

> #Constants for Makeham's Law

> A=0.0001

> B=0.01

> C=0.1

> #Calculate mortality rates using Makeham's Law

> mortality3rates

[1] 0.01010000 0.01115171 0.01231403 0.01359859 0.01501825

[6] 0.01658721 0.01832119 0.02023753 0.02235541 0.02469603


[11] 0.02728282 0.03014166 0.03330117 0.03679297 0.04065200

[16] 0.04491689 0.04963032 0.05483947 0.06059647 0.06695894

[21] 0.07399056 0.08176170 0.09035013 0.09984182 0.11033176

[26] 0.12192494 0.13473738 0.14889732 0.16454647 0.18184145

[31] 0.20095537 0.22207951 0.24542530 0.27122639 0.29974100

[36] 0.33125452 0.36608234 0.40457304 0.44711184 0.49412449

[41] 0.54608150 0.60350288 0.66696331 0.73709794 0.81460869

[46] 0.90027131 0.99494316 1.09957172 1.21520418 1.34299780

[51] 1.48423159 1.64031907 1.81282242 2.00346810 2.21416416

[56] 2.44701932 2.70436407 2.98877401 3.30309560 3.65047468

[61] 4.03438793 4.45867770 4.92759041 5.44581910 6.01855038

[66] 6.65151633 7.35105189 8.12415825 8.97857292 9.92284716

[71] 10.96643158 12.11977074 13.39440764 14.80309928 16.35994430

[76] 18.08052414 19.98205895 22.08357992 24.40611978 26.97292328

[81] 29.80967987 32.94478075 36.40960307 40.23882394 44.47076748

[86] 49.14778840 54.31669591 60.02922217 66.34254006 73.31983539

[91] 81.03093928 89.55302703 98.97139059 109.38029208 120.88390730

[96] 133.59736830

>

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy