-
Notifications
You must be signed in to change notification settings - Fork 4
/
AKI project data reshaping in STATA
45 lines (39 loc) · 1.41 KB
/
AKI project data reshaping in STATA
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
ssh username/domain
module load apps/binapps/stata/14-iph
stata
insheet using "SIR_Rinput.csv"
set more {off} [, permanently]
destring studyid dos stcreat age sex mdrd , replace
#DEFINE 2012 AS INDEX YEAR WHERE DATE IS PROVIDED IN STATA TIME
gen yr12=1 if year(dos)==2012
label var yr12 "sample from 2012"
gen post2009=1 if dos>=17899
duplicates tag studyid dos, gen(dup)
gen dupyn=dup>=1
tab dupyn
ren stcreat oldstcreat
gsort studyid dos -oldstcreat
by studyid dos: gen n=_n
gen newstcreat=oldstcreat if n==1
ren newstcreat stcreat
drop if stcreat==.
count if yr12==1
tabstat studyid if age>39& age<70& protagonistind==1 & yr12==1, stat(n)
label var post2009 "sample after 1st jan 2009"
bysort studyid post2009(dos): gen protagonist=_n if post2009==1
bysort studyid yr12(dos): gen protagonist2=_n if yr12==1
tabstat protagonist, stat (n, mean, median, min, max)
*the maximum value of "protagonist is the number of times you need to do your loop to ensure each sample is compared with all others.
tabstat mdrd if protagonist==1 & mdrd>=30&mdrd<=45, stat(count)
#CALL NUMBER OF BLOOD TESTS IN INDEX YEAR
tabstat protagonist2, stat (n, mean, median, min, max)
duplicates tag studyid dos, gen(dup)
gen dupyn=dup>=1
tab dupyn
ren stcreat oldstcreat
gsort studyid dos -oldstcreat
by studyid dos: gen n=_n
gen newstcreat=oldstcreat if n==1
ren newstcreat stcreat
drop if stcreat==.
tabstat protagonist2, stat (n, mean, median, min, max)