proc logistic过程步整理

一、代码示例(部分参数)

proc logistic data=Adpasi outest=betas covout;

class SEX DSACTST CMMTXFL CMAPTFL;

model AVALC(event='Y')= AGE SEX BMIBL PROGCLS DSACTST BSABL CMMTXFL CMAPTFL

/selection=stepwise slentry=0.05 slstay=0.1 details lackfit;

output out=pred p=phat lower=lcl upper=ucl predprob=(individual crossvalidate);

ods output Association=Association;

run;

二、proc logistic中的参数解析

1、data= :指示输入数据集名称

2、outest=:给出系数矩阵

三、class语句解析

        指示后面的model语句中模型中变量,有哪些是字符型的分类变量。该语句必须在model之前生效。class语句后面可以跟参数,参数有两种表达方式。第一种是跟在每一个变量后面,如此不同的变量可以跟不同的参数,以参数missing,变量var为例,形如:class var1 missing var2 var3 missing;第二种是跟在所有变量后面,作为通用的参数。形如:class var1 var2 var3/missing;

        class语句可以在各个过程步通用。

1、param=:该参数有若干可选项

GLM: 该取值作用目前不得而知,但是它是后面使用 lsmeans 语句的前提,只有该选项设定后 lsmeans语句才会执行。

四、model语句及其参数解析

model语句中的主体是一个等式,它是多过程步公用的语句。等式的左边(AVALC)是反应变量,右边是用来预测的控制变量(AGE SEX BMIBL PROGCLS DSACTST BSABL CMMTXFL CMAPTFL)。

1、event=:在AVALC后面的(event='Y')作用是指示把AVALC中的哪个值作为分析的目标类别。event=除了直接设置值外,还有两种选项:event=first和event=last,分别表示选取AVALC中排序在第一个的那个值和最后一个值。event选项只能在反应变量是二分类的情况下使用。

2、selection=:用来设定把变量加入模型的方法。当等式右边的控制变量不是一个的时候,简单的说,就会涉及到一个问题:哪个变量先进模型?因为一般来说自变量越少,对结果的解释力度就越弱,所以简单的想法是尽量多的在等号右边加入变量。然而,变量虽多,但每个变量对最终解释力度的贡献不同,有些甚至可能完全没有贡献,所以需要对变量要有个筛选。

选项有5个,对应的是五种筛选方法,括号里是选项简写:

        (1)BACKWARD(B):向后剔除法。将所有变量(AGE SEX BMIBL PROGCLS DSACTST BSABL CMMTXFL CMAPTFL)全部放到模型中检验,检查解释力度最低的那个变量,然后和事先设定的剔除的检验水准进行比较,进行剔除。然后重复上述过程,知道没有可以剔除的变量为止。

        (2)FORWARD(F):向前剔除法。与向后剔除法相反,先选择其中一个变量进入模型中,检查该变量的解释力度,如果符合进入的检验水准则加入模型,然后重复上述过程,知道新变量对总体解释力度没有贡献为止。

        (3)NONE(N):SAS的默认方法。

        (4)STEPWISE (S):逐步回归法。是对向后剔除法和向前剔除法的综合应用。先用向前剔除法加入变量 ,然后再用向后剔除法排除不合适的变量。这是一个反复的过程,直到既没有显著的解释变量选入回归方程,也没有不显著的解释变量从回归方程中剔除为止。以保证最后所得到的解释变量集是最优的。

        (5)SCORE:最优子集法。顾名思义,遍历了所有可能的变量组合,因此筛选出的特征必定是最优的。但缺点是计算量非常大。

3、slentry=:只有当selection=选项为STEPWISE和FORWARD的时候才能用,用来设定进入的检验水准。可简写成 sle=。

4、slstay=:只有当selection=选项为STEPWISE和BACKWARD的时候才能用,用来设定剔除的检验水准。可简写成 sls=。

5、details:显示逐步选择过程,但是当SELECTION=NONE的时候无效

6、lackfit:执行Hosmer和Lemeshow拟合优度检验。

五、lsmeans 语句及其参数解析

lsmeans 语句的语法  lsmeans 变量1 变量2 ....... /参数1 参数2;他的使用需要class中使用 param=GLM 才能使用。该语句可以在一个过程步中重复使用,这样有一个好处,因为斜杠后面参数是所有变量共享的,如果想给每个变量定制个性化参数就不能只用一个lsmeans语句。lsmeans语句参数有:

1、ilink:用来设置生成均值,均值标准误。

2、cl:用来设置上下限,包括均值上下限,均值标准误的上下限。

3、pdiff:用来生成差分表,设置该参数会增加一个差分表。

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 214,172评论 6 493
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 91,346评论 3 389
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 159,788评论 0 349
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 57,299评论 1 288
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 66,409评论 6 386
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 50,467评论 1 292
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 39,476评论 3 412
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 38,262评论 0 269
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 44,699评论 1 307
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 36,994评论 2 328
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 39,167评论 1 343
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 34,827评论 4 337
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 40,499评论 3 322
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 31,149评论 0 21
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 32,387评论 1 267
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 47,028评论 2 365
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 44,055评论 2 352

推荐阅读更多精彩内容