Impact Calculation Tool for R

From Opasnet
Revision as of 10:52, 8 July 2013 by Virpi Kollanus (talk | contribs) (Created page with "=Life table calculations= ==Input data== ===Population structure in the beginning of the assessment follow-up period (pop_data)=== <t2b name="population" index="Age" desc="...")
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

Life table calculations

Input data

Population structure in the beginning of the assessment follow-up period (pop_data)

population(-)
ObsAgeResultKuvaus
1056683
2156683
3256683
435668
5456683
6560615
7660615
8760615
9860615
10960615
111066167
121166167
131266167
14136616
151466167
161563786
171663786
181763786
191863786
201963786
212066423
222166423
232266423
242366424
252466423
262565882
272665882
282765882
292865882
302965882
313061495
323161495
333261495
343361495
353461495
363572474
373672474
383772474
393872474
403972474
414075917
424175917
434275917
444375917
454475917
464576977
474676977
484776977
494876977
504976977
515080206
525180206
535280206
545380206
555480206
565580291
575680291
585780291
595880291
605980291
616054300
626154300
636254300
646354300
656454300
666548077
676648077
686748077
696848077
706948077
717041475
727141475
737241475
747341475
757441475
767534987
777634987
787734987
797834987
807934987
818023300
828123300
838223300
848323300
858423300
868511292
878611292
888711292
898811292
908911292
91904394
92914394
93924394
94934394
95944394
9695886
9796886
9897886
9998886
10099886

Annual birth rate (birth_rate)

You have error(s) in your data:

Number of indices and result cells does not match

birth.rate(-)
ObsFollow-up periodUnitResultKuvaus
12010 57000
22011 57000
32012 57000
42013 57000
52014 57000
62015 57000
72016 57000
82017 57000
92018 57000
102019 57000
112020 57000
122021 57000
132022 57000
142023 57000
152024 57000
162025 57000
172026 57000
182027 57000
192028 57000
202029 57000

Mortality rate

Mortality risk (mort_risk)

Mortality rate/Population

Start year (start-year)

2010

Follow-up time in years (followup_time)

20


Analytica codes

Follow-up year (fu_year)

sequence(Start_year,Start_year+(Followup_time-1),1)


Year in life table (year_lt)

sequence(Start_year,Start_year+Followup_time+99,1)


Follow-up period in 5-year time steps (fu_period)

sequence(Start_year,Start_year+(Followup_time-1),5)


===5-year period in life table (period_lt)

sequence(Start_year,Start_year+Followup_time+99,5)


Age of child (age_child)

sequence(0,4,1)



Population in time, child (pop_in_time_child)

var k: Birth_rate[Fu_year=Year_lt]; k:= if k = null then 0 else k; var a:= if @Year_lt = 1 then Pop_data else (if @Age=1 then k else 0); a:= a[Age=age_child];

var j:= Mort_risk[Age=age_child];

j:=j[Fu_period=Period_lt]; j:= Si_pi(j, 5, Period_lt, Year_lt, Year_help)*5;

j:= if j < 0 then 0 else j; j:= if j > 1 then 1 else j;

j:= 1-j;

var x:= 1; while x<= min([size(age_child),size(Year_lt)]) do ( var b:= a*j; b:= b[@age_child=@age_child-1, @Year_lt=@Year_lt-1]; a:= if b=null then a else b; x:= x+1);

sum(if Year_lt = period_vs_year then a else 0,Year_lt)


Population in time, beginning of time step (pop_in_time_beg)

var a:= sum(if floor(Age/5)+1 = @Age_cat then Pop_data else 0 , Age); a:= if @Age_cat=1 then sum(Pop_in_time_child, Age_child) else (if @period_lt = 1 then a else 0);

var j:= sum(Mort_risk_scenarios,Mort_endpoint); j:=if Mort_effect = 'Acute' then j[Exposure='BAU'] else j; j:=j[Fu_period=Period_lt]; j:= if j = null then j[Exposure='Reference',Period_lt=max(Fu_period)] else j;

j:= if j < 0 then 0 else j; j:= if j > 1 then 1 else j;

j:= 1-j;

j:= sum(if floor(Age/5)+1 = @Age_cat then j else 0 , Age)/5; var m:=j[@Age_cat=@Age_cat+1]; m:= if m=null then 0 else m; var n:=((j^5)+(j^4*m)+(j^3*m^2)+(j^2*m^3)+(j*m^4))/5;

var x:= 1; while x<= min([size(Age_cat),size(Period_lt)]) do ( var b:= a*n; b:= b[@Age_cat=@Age_cat-1, @Period_lt=@Period_lt-1]; a:= if b=null then a else b; x:= x+1); a