Hello all:
I was trying to do a forest plot of HR and 95%CI and trying to sort by HR in the plot. But I am struggling. Was wondering if it had to do with me assigning those spaced column names. Tried xsvmat and mkmat etc but I am not getting ahead. Would appreciate any input. have attached a sample of my code with the stan3 dataset.
I was trying to do a forest plot of HR and 95%CI and trying to sort by HR in the plot. But I am struggling. Was wondering if it had to do with me assigning those spaced column names. Tried xsvmat and mkmat etc but I am not getting ahead. Would appreciate any input. have attached a sample of my code with the stan3 dataset.
Code:
webuse stan3, clear
gen byte age40 = age>40
stset stime, failure(died)
*Direct forest plot of univariable Cox models without resorting to R for the final work.
local list age40 transplant wait surgery
local n: word count `list'
local names ""Age>40 yrs" "Transplant" "Wait Time" "Surgery""
matrix drop _all
matrix C = J(3,`n',.)
matrix rownames C = HR ll95 ul95
matrix colnames C = `names'
local i 0
foreach v of var `list' {
local ++ i
stcox `v'
matrix C[1,`i'] = r(table)[1,1] \ r(table)[5,1] \ r(table)[6,1]
}
matrix list C
coefplot ///
(matrix(C), ci((2 3)) msymbol(S) ciopts(recast(rcap))), ///
xline(1, lwidth(small) lcolor(gray) lpattern(dash)) ///
xlabel(,notick) ///
ylabel(,notick) ///
title("HR of uni-variable models as a forest plot", ///
span size(*.9) linegap(1.5) margin(medium)) ///
graphregion(margin(0 10 0 0)) ///
caption(" " "testing it out to see if I can do HR, 95% CI here", justification(right) position(3))

Comment