[19121] in s-news-athena
[S] Where are the actions in coxph?
daemon@ATHENA.MIT.EDU (Haicheng Tang)
Mon Sep 6 22:21:31 1999
Date: Mon, 6 Sep 1999 22:17:33 -0400 (EDT)
From: Haicheng Tang <haicheng@umich.edu>
To: s-news@wubios.wustl.edu
Message-Id: <Pine.GSO.4.10.9909062210250.1960-200000@csrvr11.sph.umich.edu>
Mime-Version: 1.0
Content-Type: MULTIPART/MIXED; BOUNDARY="-559023410-959030623-936670653=:1960"
This message is in MIME format. The first part should be readable text,
while the remaining parts are likely unreadable without MIME-aware tools.
Send mail to mime@docserver.cac.washington.edu for more info.
---559023410-959030623-936670653=:1960
Content-Type: TEXT/PLAIN; charset=US-ASCII
Hi all,
I checked how the function coxph was written
(attached). I'm very surprised to see it's so
simple and can't find any key action (maximizing
partial likelihood). All I can see are some
administrative oprations, not even dynamically
loaded programs from other language. Can someone
point out the key commands in this function?
Thanks in advance!
------------------------------
Haicheng Tang |
Dept. of Biostatistics |
1420 Washington Heights |
Ann Arbor, MI 48109 |
(734)480-3229 (H) |
936-4023(O), 936-4018(Voice) |
------------------------------
---559023410-959030623-936670653=:1960
Content-Type: TEXT/PLAIN; charset=US-ASCII; name="cox.s"
Content-Transfer-Encoding: BASE64
Content-ID: <Pine.GSO.4.10.9909062217330.1960@csrvr11.sph.umich.edu>
Content-Description:
Content-Disposition: attachment; filename="cox.s"
ZnVuY3Rpb24oZm9ybXVsYSA9IGZvcm11bGEoZGF0YSksIGRhdGEgPSBzeXMu
cGFyZW50KCksIHdlaWdodHMsIHN1YnNldCwgDQogICAgICAgIG5hLmFjdGlv
biwgZXBzID0gMC4wMDAxLCBpbml0LCBpdGVyLm1heCA9IDEwLCBtZXRob2Qg
PSBjKCJlZnJvbiIsIA0KICAgICAgICAiYnJlc2xvdyIsICJleGFjdCIpLCBz
aW5ndWxhci5vayA9IFQsIHJvYnVzdCA9IEYsIG1vZGVsID0gRiwgeCA9IEYs
IHkNCiAgICAgICAgID0gVCkNCnsNCiAgICAgICAgbWV0aG9kIDwtIG1hdGNo
LmFyZyhtZXRob2QpDQogICAgICAgIGNhbGwgPC0gbWF0Y2guY2FsbCgpDQog
ICAgICAgIG0gPC0gbWF0Y2guY2FsbChleHBhbmQgPSBGKQ0KICAgICAgICBt
JG1ldGhvZCA8LSBtJG1vZGVsIDwtIG0keCA8LSBtJHkgPC0gbSQuLi4gPC0g
TlVMTA0KICAgICAgICBtJGVwcyA8LSBtJGluaXQgPC0gbSRpdGVyLm1heCA8
LSBtJHJvYnVzdCA8LSBOVUxMDQogICAgICAgIFRlcm1zIDwtIGlmKG1pc3Np
bmcoZGF0YSkpIHRlcm1zKGZvcm11bGEsIGMoInN0cmF0YSIsICJjbHVzdGVy
IikpIGVsc2UgDQogICAgICAgICAgICAgICAgICAgICAgICB0ZXJtcyhmb3Jt
dWxhLCBjKCJzdHJhdGEiLCAiY2x1c3RlciIpLCBkYXRhID0gZGF0YSkNCiAg
ICAgICAgbSRmb3JtdWxhIDwtIFRlcm1zDQogICAgICAgIG1bWzFdXSA8LSBh
cy5uYW1lKCJtb2RlbC5mcmFtZSIpDQogICAgICAgIG0gPC0gZXZhbChtLCBz
eXMucGFyZW50KCkpDQogICAgICAgIFkgPC0gbW9kZWwuZXh0cmFjdChtLCAi
cmVzcG9uc2UiKQ0KICAgICAgICBpZighaW5oZXJpdHMoWSwgIlN1cnYiKSkN
CiAgICAgICAgICAgICAgICBzdG9wKCJSZXNwb25zZSBtdXN0IGJlIGEgc3Vy
dml2YWwgb2JqZWN0IikNCiAgICAgICAgd2VpZ2h0cyA8LSBtb2RlbC5leHRy
YWN0KG0sICJ3ZWlnaHRzIikNCiAgICAgICAgb2Zmc2V0IDwtIGF0dHIoVGVy
bXMsICJvZmZzZXQiKQ0KICAgICAgICB0dCA8LSBsZW5ndGgob2Zmc2V0KQ0K
ICAgICAgICBvZmZzZXQgPC0gaWYodHQgPT0gMCkgcmVwKDAsIG5yb3coWSkp
IGVsc2UgaWYodHQgPT0gMSkNCiAgICAgICAgICAgICAgICBtW1tvZmZzZXRd
XQ0KICAgICAgICBlbHNlIHsNCiAgICAgICAgICAgICAgICBmZiA8LSBtW1tv
ZmZzZXRbMV1dXQ0KICAgICAgICAgICAgICAgIGZvcihpIGluIDI6dHQpDQog
ICAgICAgICAgICAgICAgICAgICAgICBmZiA8LSBmZiArIG1bW29mZnNldFtp
XV1dDQogICAgICAgICAgICAgICAgZmYNCiAgICAgICAgfQ0KICAgICAgICBh
dHRyKFRlcm1zLCAiaW50ZXJjZXB0IikgPC0gMSAgICNDb3ggbW9kZWwgYWx3
YXlzIGhhcyBcTGFtYmRhXzANCiAgICAgICAgc3RyYXRzIDwtIGF0dHIoVGVy
bXMsICJzcGVjaWFscyIpJHN0cmF0YQ0KICAgICAgICBjbHVzdGVyIDwtIGF0
dHIoVGVybXMsICJzcGVjaWFscyIpJGNsdXN0ZXINCiAgICAgICAgZHJvcHgg
PC0gTlVMTA0KICAgICAgICBpZihsZW5ndGgoY2x1c3RlcikpIHsNCiAgICAg
ICAgICAgICAgICBpZihtaXNzaW5nKHJvYnVzdCkpDQogICAgICAgICAgICAg
ICAgICAgICAgICByb2J1c3QgPC0gVA0KICAgICAgICAgICAgICAgIHRlbXBj
IDwtIHVudGFuZ2xlLnNwZWNpYWxzKFRlcm1zLCAiY2x1c3RlciIsIDE6MTAp
DQogICAgICAgICAgICAgICAgb3JkIDwtIGF0dHIoVGVybXMsICJvcmRlciIp
W3RlbXBjJHRlcm1zXQ0KICAgICAgICAgICAgICAgIGlmKGFueShvcmQgPiAx
KSkNCg0KDQoJICAgICAgICAgICAgICAgICAgICAgIHN0b3AoIkNsdXN0ZXIg
Y2FuIG5vdCBiZSB1c2VkIGluIGFuIGludGVyYWN0aW9uIikNCiAgICAgICAg
ICAgICAgICBjbHVzdGVyIDwtIHN0cmF0YShtWywgdGVtcGMkdmFyc10sIHNo
b3J0bGFiZWwgPSBUKQ0KICAgICAgICAjYWxsb3cgbXVsdGlwbGVzDQogICAg
ICAgICAgICAgICAgZHJvcHggPC0gdGVtcGMkdGVybXMNCiAgICAgICAgfQ0K
ICAgICAgICBpZihsZW5ndGgoc3RyYXRzKSkgew0KICAgICAgICAgICAgICAg
IHRlbXAgPC0gdW50YW5nbGUuc3BlY2lhbHMoVGVybXMsICJzdHJhdGEiLCAx
KQ0KICAgICAgICAgICAgICAgIGRyb3B4IDwtIGMoZHJvcHgsIHRlbXAkdGVy
bXMpDQogICAgICAgICAgICAgICAgaWYobGVuZ3RoKHRlbXAkdmFycykgPT0g
MSkNCiAgICAgICAgICAgICAgICAgICAgICAgIHN0cmF0YS5rZWVwIDwtIG1b
W3RlbXAkdmFyc11dDQogICAgICAgICAgICAgICAgZWxzZSBzdHJhdGEua2Vl
cCA8LSBzdHJhdGEobVssIHRlbXAkdmFyc10sIHNob3J0bGFiZWwgPSBUKQ0K
ICAgICAgICAgICAgICAgIHN0cmF0cyA8LSBhcy5udW1lcmljKHN0cmF0YS5r
ZWVwKQ0KICAgICAgICB9DQogICAgICAgIGlmKGxlbmd0aChkcm9weCkpDQog
ICAgICAgICAgICAgICAgWCA8LSBtb2RlbC5tYXRyaXgoVGVybXNbIC0gZHJv
cHhdLCBtKVssIC0xLCBkcm9wID0gRl0NCiAgICAgICAgZWxzZSBYIDwtIG1v
ZGVsLm1hdHJpeChUZXJtcywgbSlbLCAtMSwgZHJvcCA9IEZdDQogICAgICAg
IHR5cGUgPC0gYXR0cihZLCAidHlwZSIpDQogICAgICAgIGlmKG1ldGhvZCA9
PSAiYnJlc2xvdyIgfHwgbWV0aG9kID09ICJlZnJvbiIpIHsNCiAgICAgICAg
ICAgICAgICBpZih0eXBlID09ICJyaWdodCIpDQogICAgICAgICAgICAgICAg
ICAgICAgICBmaXR0ZXIgPC0gZ2V0KCJjb3hwaC5maXQiKQ0KICAgICAgICAg
ICAgICAgIGVsc2UgaWYodHlwZSA9PSAiY291bnRpbmciKQ0KICAgICAgICAg
ICAgICAgICAgICAgICAgZml0dGVyIDwtIGdldCgiYWdyZWcuZml0IikNCiAg
ICAgICAgICAgICAgICBlbHNlIHN0b3AocGFzdGUoIkNveCBtb2RlbCBkb2Vz
bid0IHN1cHBvcnQgXCIiLCB0eXBlLCANCiAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgIlwiIHN1cnZpdmFsIGRhdGEiLCBzZXAgPSAiIikpDQog
ICAgICAgIH0NCiAgICAgICAgZWxzZSBpZihtZXRob2QgPT0gImV4YWN0IikN
CiAgICAgICAgICAgICAgICBmaXR0ZXIgPC0gZ2V0KCJhZ2V4YWN0LmZpdCIp
DQogICAgICAgIGVsc2Ugc3RvcChwYXN0ZSgiVW5rbm93biBtZXRob2QiLCBt
ZXRob2QpKQ0KICAgICAgICBpZihtaXNzaW5nKGluaXQpKQ0KICAgICAgICAg
ICAgICAgIGluaXQgPC0gTlVMTA0KICAgICAgICBmaXQgPC0gZml0dGVyKFgs
IFksIHN0cmF0cywgb2Zmc2V0LCBpbml0ID0gaW5pdCwgaXRlci5tYXggPSBp
dGVyLm1heCwgDQogICAgICAgICAgICAgICAgZXBzID0gZXBzLCB3ZWlnaHRz
ID0gd2VpZ2h0cywgbWV0aG9kID0gbWV0aG9kLCByb3cubmFtZXMobSkpDQog
ICAgICAgIGlmKGlzLmNoYXJhY3RlcihmaXQpKSB7DQogICAgICAgICAgICAg
ICAgZml0IDwtIGxpc3QoZmFpbCA9IGZpdCkNCiAgICAgICAgICAgICAgICBh
dHRyKGZpdCwgImNsYXNzIikgPC0gImNveHBoIg0KICAgICAgICB9DQogICAg
ICAgIGVsc2Ugew0KICAgICAgICAgICAgICAgIGlmKGFueShpcy5uYShmaXQk
Y29lZikpKSB7DQogICAgICAgICAgICAgICAgICAgICAgICB2YXJzIDwtICgx
Omxlbmd0aChmaXQkY29lZikpW2lzLm5hKGZpdCRjb2VmKV0NCg0KDQoJICAg
ICAgICAgICAgICAgICAgICAgIG1zZyA8LSBwYXN0ZSgiWCBtYXRyaXggZGVl
bWVkIHRvIGJlIHNpbmd1bGFyOyB2YXJpYWJsZSIsDQogICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgIHBhc3RlKHZhcnMsIGNvbGxhcHNlID0gIiAi
KSkNCiAgICAgICAgICAgICAgICAgICAgICAgIGlmKHNpbmd1bGFyLm9rKQ0K
ICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICB3YXJuaW5nKG1zZykN
CiAgICAgICAgICAgICAgICAgICAgICAgIGVsc2Ugc3RvcChtc2cpDQogICAg
ICAgICAgICAgICAgfQ0KICAgICAgICAgICAgICAgIGZpdCRuIDwtIG5yb3co
WSkNCiAgICAgICAgICAgICAgICBhdHRyKGZpdCwgImNsYXNzIikgPC0gZml0
JG1ldGhvZA0KICAgICAgICAgICAgICAgIGZpdCR0ZXJtcyA8LSBUZXJtcw0K
ICAgICAgICAgICAgICAgIGZpdCRhc3NpZ24gPC0gYXR0cihYLCAiYXNzaWdu
IikNCiAgICAgICAgICAgICAgICBpZihyb2J1c3QpIHsNCiAgICAgICAgICAg
ICAgICAgICAgICAgIGZpdCRuYWl2ZS52YXIgPC0gZml0JHZhcg0KICAgICAg
ICAgICAgICAgICAgICAgICAgZml0JG1ldGhvZCA8LSBtZXRob2QNCiAgICAg
ICAgIyBhIGxpdHRsZSBzbmVha3kgaGVyZTogYnkgY2FsbGluZyByZXNpZCBi
ZWZvcmUgYWRkaW5nIHRoZQ0KIyAgIG5hLmFjdGlvbiBtZXRob2QsIEkgYXZv
aWQgaGF2aW5nIG1pc3NpbmdzIHJlLWluc2VydGVkDQojIEkgYWxzbyBtYWtl
IHN1cmUgdGhhdCBpdCBkb2Vzbid0IGhhdmUgdG8gcmVjb25zdHJ1Y3QgWCBh
bmQgWQ0KICAgICAgICAgICAgICAgICAgICAgICAgZml0MiA8LSBjKGZpdCwg
bGlzdCh4ID0gWCwgeSA9IFksIHdlaWdodHMgPSB3ZWlnaHRzKSkNCiAgICAg
ICAgICAgICAgICAgICAgICAgIGlmKGxlbmd0aChzdHJhdHMpKQ0KICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgICBmaXQyJHN0cmF0YSA8LSBzdHJh
dGEua2VlcA0KICAgICAgICAgICAgICAgICAgICAgICAgaWYobGVuZ3RoKGNs
dXN0ZXIpKSB7DQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIHRl
bXAgPC0gcmVzaWR1YWxzLmNveHBoKGZpdDIsIHR5cGUgPSAiZGZiZXRhIiwg
DQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgY29sbGFwc2Ug
PSBjbHVzdGVyLCB3ZWlnaHRlZCA9IFQpDQogICAgICAgICMgICAgICAgICAg
ICAgICBpbmR4IDwtIG1hdGNoKGNsdXN0ZXIsIHVuaXF1ZShjbHVzdGVyKSkN
CiMgICAgICAgICAgICAgICBrICAgIDwtIGFzLnZlY3Rvcih0YWJsZShpbmR4
KSkNCiMgICAgICAgICAgICAgICBpZiAoYW55KGs+MSkpIHsNCiMgICAgICAg
ICAgICAgICAgICAgI2NvbXB1dGUgdGhlIElDQyBmb3IgdGhlIG0tcmVzaWR1
YWxzDQojICAgICAgICAgICAgICAgICAgIE4gICAgPC0gc3VtKGsgKiAoay0x
KSkNCiMgICAgICAgICAgICAgICAgICAgbXUgICA8LSBzdW0oZml0JHJlc2lk
ICogKGstMSlbaW5keF0pL04gICAjZ3JhbmQgbWVhbg0KIyAgICAgICAgICAg
ICAgICAgICBtdTIgIDwtIHRhcHBseShmaXQkcmVzaWQsIGluZHgsIG1lYW4p
ICAgICMgZ3JvdXAgbWVhbnMNCiMgICAgICAgICAgICAgICAgICAgc2lnICA8
LSB0YXBwbHkoKGZpdCRyZXNpZCAtIG11KV4yLCBpbmR4LCBzdW0pICAjU1MN
CiMgICAgICAgICAgICAgICAgICAgaWNjMSA8LSBzdW0oIChrKihtdTItbXUp
KV4yIC0gc2lnKSAvIHN1bSgoay0xKSpzaWcpDQojICAgICAgICAgICAgICAg
ICAgICNyYW5rIHJlc2lkdWFscw0KIyAgICAgICAgICAgICAgICAgICByciA8
LSByYW5rKGZpdCRyZXNpZCkNCiMgICAgICAgICAgICAgICAgICAgbXUgICA8
LSBzdW0ocnIgKiAoay0xKVtpbmR4XSkvTiAgICNncmFuZCBtZWFuDQojICAg
ICAgICAgICAgICAgICAgIG11MiAgPC0gdGFwcGx5KHJyLCBpbmR4LCBtZWFu
KSAgICAjIGdyb3VwIG1lYW5zDQojICAgICAgICAgICAgICAgICAgIHNpZyAg
PC0gdGFwcGx5KChyciAtIG11KV4yLCBpbmR4LCBzdW0pICAjU1MNCiMgICAg
ICAgICAgICAgICAgICAgaWNjMiA8LSBzdW0oIChrKihtdTItbXUpKV4yIC0g
c2lnKSAvIHN1bSgoay0xKSpzaWcpDQojDQojICAgICAgICAgICAgICAgICAg
IGZpdCRpY2MgPC0gYyhsZW5ndGgoayksIGljYzEsIGljYzIpDQojICAgICAg
ICAgICAgICAgICAgIG5hbWVzKGZpdCRpY2MpIDwtIGMoIm5jbHVzdCIsICJp
Y2MocmVzaWQpIiwgDQojICAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgImljYyhyYW5rKHJlc2lkKSkiKQ0KIyAg
ICAgICAgICAgICAgICAgICB9DQogICAgICAgICAgICAgICAgICAgICAgICB9
DQogICAgICAgICAgICAgICAgICAgICAgICBlbHNlIHRlbXAgPC0gcmVzaWR1
YWxzLmNveHBoKGZpdDIsIHR5cGUgPSAiZGZiZXRhIiwgDQogICAgICAgICAg
ICAgICAgICAgICAgICAgICAgICAgICAgd2VpZ2h0ZWQgPSBUKQ0KICAgICAg
ICAgICAgICAgICAgICAgICAgZml0JHZhciA8LSB0KHRlbXApICUqJSB0ZW1w
DQogICAgICAgICAgICAgICAgfQ0KICAgICAgICAgICAgICAgIG5hLmFjdGlv
biA8LSBhdHRyKG0sICJuYS5hY3Rpb24iKQ0KICAgICAgICAgICAgICAgIGlm
KGxlbmd0aChuYS5hY3Rpb24pKQ0KICAgICAgICAgICAgICAgICAgICAgICAg
Zml0JG5hLmFjdGlvbiA8LSBuYS5hY3Rpb24NCiAgICAgICAgICAgICAgICBp
Zihtb2RlbCkNCiAgICAgICAgICAgICAgICAgICAgICAgIGZpdCRtb2RlbCA8
LSBtDQogICAgICAgICAgICAgICAgZWxzZSB7DQogICAgICAgICAgICAgICAg
ICAgICAgICBpZih4KSB7DQogICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgIGZpdCR4IDwtIFgNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAg
ICAgaWYobGVuZ3RoKHN0cmF0cykpDQogICAgICAgICAgICAgICAgICAgICAg
ICAgICAgICAgICAgZml0JHN0cmF0YSA8LSBzdHJhdGEua2VlcA0KICAgICAg
ICAgICAgICAgICAgICAgICAgfQ0KICAgICAgICAgICAgICAgICAgICAgICAg
aWYoeSkNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgZml0JHkg
PC0gWQ0KICAgICAgICAgICAgICAgIH0NCiAgICAgICAgfQ0KICAgICAgICBp
ZighaXMubnVsbCh3ZWlnaHRzKSAmJiBhbnkod2VpZ2h0cyAhPSAxKSkNCiAg
ICAgICAgICAgICAgICBmaXQkd2VpZ2h0cyA8LSB3ZWlnaHRzDQogICAgICAg
IGZpdCRmb3JtdWxhIDwtIGFzLnZlY3RvcihhdHRyKFRlcm1zLCAiZm9ybXVs
YSIpKQ0KICAgICAgICBmaXQkY2FsbCA8LSBjYWxsDQogICAgICAgIGZpdCRt
ZXRob2QgPC0gbWV0aG9kDQogICAgICAgIGZpdA0KfQ0K
---559023410-959030623-936670653=:1960--
-----------------------------------------------------------------------
This message was distributed by s-news@wubios.wustl.edu. To unsubscribe
send e-mail to s-news-request@wubios.wustl.edu with the BODY of the
message: unsubscribe s-news
