hw1sol#
Here’s how I did hw1:
qui {
if 0 { //background: ph340.600.01, 04/21/2023-04/22/2023
1. hw1 .dofile: lookout for a blank -> noi di ""
2. output not exactly as required, instructed
3. spot the subtle differences for practice
4. these are more parsimonious & i'll certainly
5. ask next years' students to do it this way :)
}
if 1 { //methods: macros, logfile, settings
global workdir `pwd'
capture log close
log using "${workdir}hw1.lastname.firstname.log", replace
cls
set more off
}
if 2 { //results: data, shape, etc.
if 00 { //stata dropdown menu
1. with a mouse i selected (after trying other options)
2. file > import > text data (delimited...)
3. located hw1.txt in my downloads folder
4. imported it and copied & pasted the syntax, edited
5. from Stata's output window into this do file :)
}
import delimited "hw1.txt", clear
}
if 3 { //conclusions: questions, code, output
//q1
noi di "Question 1: There are `c(N)' records in the dataset."
noi di ""
//q2
forvalues i=1/2 {
sum init_age if female==`i'-1, detail
#delimit ;
local m_iqr`i': di %2.0f r(p50)
" [" %2.0f r(p25)
"-" %2.0f r(p75)
"]"
;
#delimit cr
}
noi di "Question 2: The median [IQR] age is `m_iqr1' among males and `m_iqr2' among females."
noi di ""
//q3
forvalues i=1/2 {
sum prev if female==`i'-1, detail
local perc`i': di %2.1f r(mean)*100
}
noi di "Question 3: `perc1'% among males and `perc2'% among females have history of previous transplant."
noi di ""
//q4
g htn=dx=="4=Hypertensive"
label define htn_lab 0 "No" 1 "Yes"
label values htn htn_lab
noi di "Question 4:"
noi tab htn
noi di ""
//q5
capture program drop question5
program define question5
#delimit ; //use throughout program
qui { ; //line 1 quietly doesn't apply inside this program
//label variables with output in mind
lab var init_age "Age, median [IQR]";
local age_lab: var lab init_age;
lab var prev "Previous transplant, %" ;
local prev_lab: var lab prev;
forvalues i=1/2 { ; //columns 1 & 2
count if !missing(dx) & female==`i'-1;
//row1
count if female==`i'-1;
local female`i'_N=r(N);
local row1: di "Question 5"
_col(30) "Males (N=`female1_N')"
_col(50) "Females (N=`female2_N')"
;
//row2
sum init_age if female==`i'-1,
detail; //copy&paste from q2, edit
local m_iqr`i': di %2.0f r(p50)
" [" %2.0f r(p25)
"-" %2.0f r(p75)
"]"
;
local row2: di "`age_lab'"
_col(30) "`m_iqr1'"
_col(50) "`m_iqr2'"
;
//row3
sum prev if female==`i'-1 ;
local per_prev`i': di %2.1f r(mean)*100 ;
local row3: di "`prev_lab'"
_col(30) "`per_prev1'"
_col(50) "`per_prev2'"
;
} ;
//rows4_13
split dx, p("=") ; //from chapter: delimit
destring dx1, replace ;
lab var dx1 "Cause of ESRD, %" ;
local varlab: var lab dx1 ;
label def varlab
1 "Glomerular"
2 "Diabetes"
3 "PKD"
4 "Hypertensive"
5 "Renovascular"
6 "Congenital"
7 "Tubulo"
8 "Neoplasm"
9 "Other"
;
lab values dx1 varlab;
local row4: di "`varlab'" ;
local vallab: value label dx1
; //debug: chatGPT moved this from line 137 to 152!!!
forvalues i=1/2 { ; //columns 1 & 2
levelsof dx1 if female==`i'-1,
local(diagnosis) ; //variable-level
global N_`i'=r(N) ;
local row=5 ; //based on Q5. template
foreach l of numlist `diagnosis' { ;
local dxcat: lab `vallab' `l' ; //alliterative
sum dx1 if dx1==`l' & female==(`i'-1) ;
local col_`i'_`row': di %2.1f r(N)*100/${N_`i'} ;
//indent the lab `dxcat'
local row`row': di " `dxcat'"
_col(30) "`col_1_`row''"
_col(50) "`col_2_`row''"
;
local row = `row' + 1 ; //tracks rows 5-13
} ;
} ;
forvalues i=1/13 { ; //rows1-13
noi di "`row`i''";
} ;
} ;
#delimit cr
end
noi question5
//q6
logistic received_kt init_age female
matrix define m=r(table)
//row1
noi di ""
noi di "Question 6"
//row2
local row2: di "Variable" _col(30) "OR" _col(35) "(95% CI)"
noi di "`row2'"
//rows3-4
lab var init_age "Age"
local age_lab: var lab init_age
lab var female "Female"
local female_lab: var lab female
local row=3
local col=1
foreach v of varlist init_age female {
local `v'_lab: var lab `v'
#delimit ;
local row`col': di "``v'_lab'" %3.2f _col(30) m[1,`col']
%3.2f _col(35)
"(" %3.2f m[5,`col']
"-" %3.2f m[6,`col']
")"
;
#delimit cr
//noi di "`row`num''"
local row=`row' + 1
local col=`col' + 1
}
noi di "`row1'"
noi di "`row2'"
noi di ""
//q7
noi di "Question 7: This regression included `e(N)' observations whereas the study dataset has `c(N)' observations in total."
noi di ""
//q8
noi di "Question 8: I estimate that it took me 24 hours to complete this assignment."
noi di ""
}
log close
}
Notes:
qui {: Click on the[-]in the margin, adjacent toqui:[-]qui {
qui {
* Does your entire .do file collapse as shown above? Good!
* Some code appears below `qui {` after you collapse using `[-]`?
* Then you've departed from the guidelines suggested in the first-half of this class
* All the code should be within the quietly braces
if X {: Click on the[-]in the margin, adjacent to eachif:[-] if {
qui {
if 0 {
if 1 {
if 2 {
$\vdots$
log close
}
* Does your .do file collapse as shown? Excellent!
* Any other lines of code visible after collapsing using `[-]`?
* Then you've departed from the .do file convenetions of this class
code-blocks * indentation within code-blocks * annotation and documentation! * how does this enhance readability of code?
loops *
forvalues i=1/2 {*foreachof numlistdiagnosis' {*forvalues i=1/12 {*foreach v of varlist init_age female {
5. unorthodox
* extensive use of `#delimit ;`
* how does this enhance vertical vs. horizontal growth of code?
* within `if 3 {` code-block
* specifically within `program define question5`
Cell In[1], line 1
5. unorthodox
^
SyntaxError: invalid syntax
macros * c() class macros in q1, q7 * r() class macros in q2, q3, q5, q6 * e() class macros in q7 global macro in
if 1 {local macros likei,l,vin loops local macros likem_iqr,perci’’,age_lab,prev_lab,femalei’_N’, etc.
miscellaneous *
split dx, p("=")*matrix define m=r(table)*local row=row’ + 1*noi di “”*levelsof`