由买买提看人间百态

boards

本页内容为未名空间相应帖子的节选和存档,一周内的贴子最多显示50字,超过一周显示500字 访问原贴
Statistics版 - 请教HAZARD MODEL的sas code
相关主题
censored observation 在hazard model里面是怎么处理的?请教一个fit Cox PHREG model的问题
survival 中关于censoring的问题need some help about survival analysis in SAS code
问个censored的survival data概念请教R有没有package做Linear Mixed Model Survival Analysis (continuous,right censored DV)
About the model adequacy test on Parametric survival analysisconverge problem in proc phreg
time-to-event survival model vs time sliced classification modelsSAS question (紧急求助,在线等)
为什么验证survival predictive model 都用c-index?survival analysis
问一个SAS Life-table的问题 !求survival analysis课本中的数据
[急求助] survival analysis (SAS)【请教】关于survival的问题
相关话题的讨论汇总
话题: censor话题: time话题: sas话题: model话题: hazard
进入Statistics版参与讨论
1 (共1页)
D*********Y
发帖数: 3382
1
proc phreg data=asa.hmohiv;
model time*censor(0)=drug age;
run;
不是很懂。我猜data是不是这样的。
data 在不同time 收集。初始有N个observation。
比如,time 1的时候有n1个observation 死了,那么就设censor=0,没死的censor=1;
死掉的从smaple里面剔除。
time 2的时候就只有N-N1个observation。 N2个死了的就射censor=0,没死的censor=1;
。。。。。。。。。。。etc
这个time是不是也只是个event的indicator。具体数据是多少是不影响regression结果
的。
我的理解对不对呢?
A*****a
发帖数: 1091
2
不是您想的这个意思,
survival里面时间很重要,censor为0的那些observation提供的是完整的数据,censor
为1的那些提供的是不完整的数据。

=1;

【在 D*********Y 的大作中提到】
: proc phreg data=asa.hmohiv;
: model time*censor(0)=drug age;
: run;
: 不是很懂。我猜data是不是这样的。
: data 在不同time 收集。初始有N个observation。
: 比如,time 1的时候有n1个observation 死了,那么就设censor=0,没死的censor=1;
: 死掉的从smaple里面剔除。
: time 2的时候就只有N-N1个observation。 N2个死了的就射censor=0,没死的censor=1;
: 。。。。。。。。。。。etc
: 这个time是不是也只是个event的indicator。具体数据是多少是不影响regression结果

D*********Y
发帖数: 3382
3
能详细说说嘛?能给个reference吗?
我看那个proportional hazard里面的survival time就是用来计算likelihood,只是比
大小,没有具体time 计入公式。
不知道只是怎么在SAS里实现的。
censor是什么意思?如果数据不完整,怎么处理的呢?

censor

【在 A*****a 的大作中提到】
: 不是您想的这个意思,
: survival里面时间很重要,censor为0的那些observation提供的是完整的数据,censor
: 为1的那些提供的是不完整的数据。
:
: =1;

h******s
发帖数: 3420
4
你这个完全不对。censor 是 lost following up, 跟死不死没干系。censor 的要从
denominator 去掉。
建议找个浅显一点的survival analysis, 多看几遍,直到看懂为止。

=1;

【在 D*********Y 的大作中提到】
: proc phreg data=asa.hmohiv;
: model time*censor(0)=drug age;
: run;
: 不是很懂。我猜data是不是这样的。
: data 在不同time 收集。初始有N个observation。
: 比如,time 1的时候有n1个observation 死了,那么就设censor=0,没死的censor=1;
: 死掉的从smaple里面剔除。
: time 2的时候就只有N-N1个observation。 N2个死了的就射censor=0,没死的censor=1;
: 。。。。。。。。。。。etc
: 这个time是不是也只是个event的indicator。具体数据是多少是不影响regression结果

D*********Y
发帖数: 3382
5
你和楼上的说的不一样啊。楼上说censor的是数据缺失?难道他说的数据缺失就是要被
去掉的意思?

【在 h******s 的大作中提到】
: 你这个完全不对。censor 是 lost following up, 跟死不死没干系。censor 的要从
: denominator 去掉。
: 建议找个浅显一点的survival analysis, 多看几遍,直到看懂为止。
:
: =1;

D*********Y
发帖数: 3382
6
http://en.wikipedia.org/wiki/Proportional_hazards_models
我看的是这个,我觉得我看懂了。
不知道SAS是怎么处理的。数据如何安排。还有censor这个variable怎么建立的?

【在 h******s 的大作中提到】
: 你这个完全不对。censor 是 lost following up, 跟死不死没干系。censor 的要从
: denominator 去掉。
: 建议找个浅显一点的survival analysis, 多看几遍,直到看懂为止。
:
: =1;

g**u
发帖数: 205
7
不过censor(0), censor(1) indicate 什么 都是可以自己规定的吧?
一般censor(1)是代指那些drop off,lost follow up和survive to the end的, 但
是如果自己规定censor(0)才是的话也没办法。time一般是指survival time,是一个
时间段,比如7天,30天。
D*********Y
发帖数: 3382
8
搞不清楚,抓心饶肝。能给个data sample (for sas programming)看看什么摸样吗
?给包子。

【在 g**u 的大作中提到】
: 不过censor(0), censor(1) indicate 什么 都是可以自己规定的吧?
: 一般censor(1)是代指那些drop off,lost follow up和survive to the end的, 但
: 是如果自己规定censor(0)才是的话也没办法。time一般是指survival time,是一个
: 时间段,比如7天,30天。

A*****a
发帖数: 1091
9
举个列子,灯泡的寿命,灯泡迟早会坏,质量好的能用久些,质量差的很快就坏,某些
质量好的,可能观察期结束了还没坏,我们只知道它们用得超过了这个时间段,但是真
实的寿命并不知道,这些就是被censor了的数据(这里专指right censored)。
survival中真正想知道的是某种灯泡到底能用多久(质量有多好),所以我说没censor
的数据是完整的,censor了的是不完整的。
具体sas怎么处理,要不您去看下sas官网?http://support.sas.com/documentation/cdl/en/statug/63347/HTML/default/viewer.htm#statug_phreg_sect001.htm

【在 D*********Y 的大作中提到】
: http://en.wikipedia.org/wiki/Proportional_hazards_models
: 我看的是这个,我觉得我看懂了。
: 不知道SAS是怎么处理的。数据如何安排。还有censor这个variable怎么建立的?

D*********Y
发帖数: 3382
10
您咋这么客气啊!

censor

【在 A*****a 的大作中提到】
: 举个列子,灯泡的寿命,灯泡迟早会坏,质量好的能用久些,质量差的很快就坏,某些
: 质量好的,可能观察期结束了还没坏,我们只知道它们用得超过了这个时间段,但是真
: 实的寿命并不知道,这些就是被censor了的数据(这里专指right censored)。
: survival中真正想知道的是某种灯泡到底能用多久(质量有多好),所以我说没censor
: 的数据是完整的,censor了的是不完整的。
: 具体sas怎么处理,要不您去看下sas官网?http://support.sas.com/documentation/cdl/en/statug/63347/HTML/default/viewer.htm#statug_phreg_sect001.htm

相关主题
为什么验证survival predictive model 都用c-index?请教一个fit Cox PHREG model的问题
问一个SAS Life-table的问题 !need some help about survival analysis in SAS code
[急求助] survival analysis (SAS)请教R有没有package做Linear Mixed Model Survival Analysis (continuous,right censored DV)
进入Statistics版参与讨论
A*****a
发帖数: 1091
11
哈哈,我这不是半桶水么,心里发虚着呢:P

【在 D*********Y 的大作中提到】
: 您咋这么客气啊!
:
: censor

D*********Y
发帖数: 3382
12
在sas dataset里面是不是要在每个time (比如说第7天,第八天)都要把survivals
列出来。然后time大于7天的就是censor的?设个0.是这样吗?

【在 g**u 的大作中提到】
: 不过censor(0), censor(1) indicate 什么 都是可以自己规定的吧?
: 一般censor(1)是代指那些drop off,lost follow up和survive to the end的, 但
: 是如果自己规定censor(0)才是的话也没办法。time一般是指survival time,是一个
: 时间段,比如7天,30天。

g**u
发帖数: 205
13
比如说你design了一个study,是10天的,如果能活过10天的,就是censor的,一般是
用1来表示,censor(1)。在SAS中,你可以自己输入一个column来专门indicate
censor,0 或者 1。
data a;
input ID time censor;
datalines;
1 5 0
2 6 0
3 10 1
4 8 0
5 4 0
;
run;
大概意思就是这样的.

【在 D*********Y 的大作中提到】
: 在sas dataset里面是不是要在每个time (比如说第7天,第八天)都要把survivals
: 列出来。然后time大于7天的就是censor的?设个0.是这样吗?

y**3
发帖数: 267
14
Lz need understand how Time is calculated. Here Time is a duration, not a
time point.
For example, for a study which needs recuit 100 lung cancer patients in a
medical center.The whole study duration is set 5 years. Suppose the starting
time point fr computing Time is the date when the patient is diagnosed lung
cancer( they are different for different patients), the event time point is
the patient dies of lung cancer(censor=1 as in the model statement LZ gave)
. At the date of the study closes, some lung cancer patients are still
alive, so they are censored(censor=0).
how to caluculte the variable Time in the model statement:
for the patient died,Time=date of death-date of diagnosis;
for patient cnesored, Time=date of last avaiable follow_up -date of
diagnosis
Hope this helps
D*********Y
发帖数: 3382
15
原来是这样,我以为1,2,3,4.。。。。10,都得check一下呢。

【在 g**u 的大作中提到】
: 比如说你design了一个study,是10天的,如果能活过10天的,就是censor的,一般是
: 用1来表示,censor(1)。在SAS中,你可以自己输入一个column来专门indicate
: censor,0 或者 1。
: data a;
: input ID time censor;
: datalines;
: 1 5 0
: 2 6 0
: 3 10 1
: 4 8 0

1 (共1页)
进入Statistics版参与讨论
相关主题
【请教】关于survival的问题time-to-event survival model vs time sliced classification models
请教个用SAS lifereg interval-censored survival analysis后作图的问题为什么验证survival predictive model 都用c-index?
请教个多元回归中的问题问一个SAS Life-table的问题 !
包子求助Survival问题,很菜,请大家包涵[急求助] survival analysis (SAS)
censored observation 在hazard model里面是怎么处理的?请教一个fit Cox PHREG model的问题
survival 中关于censoring的问题need some help about survival analysis in SAS code
问个censored的survival data概念请教R有没有package做Linear Mixed Model Survival Analysis (continuous,right censored DV)
About the model adequacy test on Parametric survival analysisconverge problem in proc phreg
相关话题的讨论汇总
话题: censor话题: time话题: sas话题: model话题: hazard