经管之家App
让优质教育人人可得
立即打开
全部版块
我的主页
›
论坛
›
数据科学与人工智能
›
数据分析与数据科学
›
SAS专版
求多因素Logistic回归结果做森林图程序。
楼主
959712305
10266
2
收藏
2016-01-18
悬赏
200
个论坛币
未解决
各位大神,我想利用多因素Logistic回归后得到的各因素的OR值、OR值的95%CI和P值做森林图,如这张网上找到的图。我想要的结果与其类似,但是想要每个因素的第一个水平作参照,起OR值为1或无,参照水平无置信区间。感谢帮忙的人。
%$X2T7$NH2THQ~7ZJO2FLH1.jpg
原图尺寸 89.13 KB
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
全部回复
沙发
wang1839
2016-6-30 23:03:27
这个简单,我做一个类似的吧
data or;
input param $10. or min max p;
cards;
sex 0.81 1.23 0.53 0.33
age 0.70 1.53 0.52 0.06
bmi 0.85 1.39 0.64 0.29
eduacation 0.89 1.10 0.70 0.76
;
run;
data tot;
set or;
CI=put(min,7.2)||' to '||put(max,7.2. -l);
ord=_n_;
retain constant 1;
run;
ods path work.templ(update) sashelp.tmplmst(read) sashelp.styles(read);
proc template;
define style mystyle;
parent = Styles.default;
style parskip / fontsize = 0pt;
class text / 'continued' = " ";
class Continued /
frame=void
pretext = _undef_
borderspacing = 0
padding = 0 ;
class GraphBackground/
backgroundcolor = white
color = white;
class LayoutContainer /
padding = 0
borderspacing = 0
borderwidth = 0
frame = void
rules = none
backgroundcolor = _undef_;
style GraphFonts from GraphFonts/
'GraphDataFont'=('Times New Roman' ,12pt)
'GraphUnicodeFont'=('Times New Roman' ,12pt)
'GraphValueFont'=('Times New Roman' ,12pt)
'GraphLabelFont'=('Times New Roman' ,12pt)
'GraphFootnoteFont'=('Times New Roman' ,12pt) /*Times New Roman*/
'GraphTitleFont'=('Times New Roman' ,12pt)
'GraphAnnoFont'=('Times New Roman' ,12pt);
class graph/
fontweight=bold;
end;
run;
proc template;
define statgraph ForestPlot;
begingraph / designwidth=9in designheight=16in;
layout lattice / columns=3 columngutter=0 columnweights=(0.10 .60 .30 );
layout overlay / walldisplay=none border=false
yaxisopts=(reverse=true type=discrete display=none)
xaxisopts=(display=none offsetmin=0.25 offsetmax=0.25);
entry halign=left " Parameter" / location=outside valign=top ;
scatterplot y=ord x=eval(constant*1) / markercharacter=param;
endlayout;
layout overlay /walldisplay=(fill)
yaxisopts=(display=none reverse=true offsetmin=0.00 offsetmax=0.10)
xaxisopts=(type=log offsetmin=0.00 offsetmax=0.00
label="<-----Protect Factor Harm Factor----->"
logopts=(base=10 minorticks=true viewmin=0.01 viewmax=100 ) );
entry halign=center "OR and 95% CIs" / location=outside valign=top ;
scatterplot y=ord x=or / xerrorupper=max xerrorlower=min
markerattrs=graphdata1(symbol=diamondfilled size=1pct);
referenceline x=1 / lineattrs=(pattern=solid);
referenceline x=0.1 / lineattrs=(pattern=shortdash);
referenceline x=10 / lineattrs=(pattern=shortdash);
endlayout;
layout overlay / walldisplay=none border=false
yaxisopts=(reverse=true type=discrete display=none)
xaxisopts=(display=none offsetmin=0.25 offsetmax=0.25);
entry halign=left " OR 95% CI P" / location=outside valign=top ;
scatterplot y=ord x=eval(constant*1) / markercharacter=or;
scatterplot y=ord x=eval(constant*2) / markercharacter=ci;
scatterplot y=ord x=eval(constant*3) / markercharacter=p;
endlayout;
endlayout;
endgraph;
end;
run;
/*ods graphics on / height=5.5 in width=10 in border=off;*/
ods rtf file="d:\test.rtf" style=mystyle;
proc sgrender data=tot template=ForestPlot;
run;
ods rtf close;
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
藤椅
一吴玉
2017-12-16 20:17:15
请教大神,用sgplot的程序是什么样的,谢谢!
扫码加我 拉你入群
请注明:姓名-公司-职位
以便审核进群资格,未注明则拒绝
相关推荐
logistic回归的问题(或因子分析问题)
[分享]多元统计分析之logistic回归
关于Logistic回归中的结果解释问题,急切等待回复!多谢!
恳请各位解惑,一个关于logistic回归方程检验的问题
sas logistic回归最大似然比不存在是怎么原因啊
如何用R软件做多类别Logistic回归
Logistic回归
求助,sas做配对logistic回归,不知道问题出在哪
加入资产对数后变得不显著
有做信用评分logistic回归,需要数据分箱(Binning)求帮助的吗?
栏目导航
SAS专版
真实世界经济学(含财经时事)
悬赏大厅
数据交流中心
金融实务版
金融学(理论版)
热门文章
2026“课题申报”抢跑号角的已吹响!国社科 ...
Nature点赞!哈佛MIT最新作:AI科学家时代来 ...
CDA 认证考试大纲 2025 重磅更新:一二级考 ...
GTAP11运行扩展数据库出错,希望高手指点。
2025秋季大摩宏观团队闭门会议纪要
英文书籍
建筑的想象之整理补充笔记
英文书籍
中国移动:智能体互联网技术白皮书2025
超越普里瓦洛夫无穷乘积与它对解析函数的应 ...
推荐文章
10月重磅来袭|《打造Coze/Dify专属学术智能 ...
高校老师和学生都在偷偷上的智能体课,到底 ...
最快1年拿证,学费不足5W!热门美国人工智能 ...
关于如何利用文献的若干建议
关于学术研究和论文发表的一些建议
关于科研中如何学习基础知识的一些建议 (一 ...
一个自编的经济学建模小案例 --写给授课本科 ...
AI智能体赋能教学改革: 全国AI教育教学应用 ...
2025中国AIoT产业全景图谱报告-406页
关于文献求助的一些建议
说点什么
分享
微信
QQ空间
QQ
微博
扫码加好友,拉您进群
各岗位、行业、专业交流群