C***l
发帖数: 2625 | 1 at proc lifetest statement, specify outsurv=,
then print the data set. |
|
s**********l
发帖数: 395 | 2 刚才上传的文档是word2007, 现在改成03版了
哪位高手能帮忙看看这个数据库对RX(drug)能不能用lifetest, 如果可以的话， SAS
Code 如何写？因为这个数据库的时间分别用年、月、日 表示的，我不知道如何写CODE
.有人跟我说可以先用transpose 处理一下，但是具体怎么用？
谢谢了。 |
|
A*********u
发帖数: 8976 | 3 rxbegday=mdy(rxbegmm, rxbegdd, rxbegyrx);
刚才上传的文档是word2007, 现在改成03版了
哪位高手能帮忙看看这个数据库对RX(drug)能不能用lifetest, 如果可以的话， SAS
Code 如何写？因为这个数据库的时间分别用年、月、日 表示的，我不知道如何写CODE
.有人跟我说可以先用transpose 处理一下，但是具体怎么用？
谢谢了。 |
|
m******3
发帖数: 734 | 4 proc phreg貌似有个proportional hazard的assumption
lifetest可以有别的option, 不一定ph |
|
c**d
发帖数: 104 | 5 Those ODS graphics in the statistical procedures are predefined template
using Graphics template language(GTL). Using ods trace on to find the name
of that graph template. Then you have two ways to modify them.
1) output SAS codes of the survival graph template.
proc template;
source Stat.Lifetest.Graphics.ProductLimitSurvival;
run;
You can find and change them.
(a) "layout overlay / yaxisopts=(label="Survival Probability")".
(b) DiscreteLegend "Survival" / title=GROUPNAME across=1 lo |
|
g******d
发帖数: 231 | 6 我用以下的SAS code作Survival Curve:
proc lifetest data=models plots=(s) cs=none graphics;
time duration*event(0);
strata age1;
run;
比较讨厌的情况是output plot中的legend没法编辑，结果图表中的horizontal axis
title和variable legend都是按照输入的数据的原本模样。还有其他一些讨厌的地方。
请问有什么办法可以改善这种状况? 我现在用的是SAS9.2.
我能想得出来的是将survival function output拷到Excel中，然后完全让Excel来做
graph。但问题是数据的样本很大，case非常多，没办法简单的copy and paste。
所以请高手予以指教. 多谢多谢! |
|
a********a
发帖数: 346 | 7 In proc lifetest, how to overlay the graph of KM from real observed survival
time and KM from predicated survival time?
Thanks |
|
f********m
发帖数: 197 | 8 相同的数据和program，在9.1.3和9.2上得出的结果不一样，差别还挺大的。有大牛遇
到过类似的问题吗？是因为改版后需要设置些option,还是9.1.3本身的lifetest有问题
啊？ |
|
A*********u
发帖数: 8976 | 9 用SAS
proc lifetest缺省计算 25% 50% %75
如果想要知道 33% 和 66% 的 point estimate 以及 95% CI
有什么option吗
还是不得不手算？
多谢多谢 |
|
h****9
发帖数: 26 | 10 Hi,当我用LIFETEST Procedure 时，我只知道Survival Standard error, 如果我想知
道 failure Standard error，然后算95%CI for failure，SAS 有办法算吗？
非常感谢
The LIFETEST Procedure
Product-Limit Survival Estimates
Survival
The LIFETEST Procedure
Product-Limit Survival Estimates
Survival
|
|
c*********n
发帖数: 87 | 11 我觉得我犯了一个大错，你在1楼提到使用lifereg对吧，但是我把lifereg和lifetest
混用了，因为lifetest是nonparametric的，常用的方法包括KM和life table （2个
group），log rank test对多组，但是lifereg是parametric的，所以可能要改成：
proc lifereg ;
class trt;
model time*status(0) = trt / dist=???;
run;
你得specify一个survial time的distribution，常用的有exponential，gamma等等，
你看看这个网页：http://support.sas.com/documentation/cdl/en/statug/63033/HTML/default/statug_lifereg_sect011.htm
run完后的output里有Analysis of Parameter Estimates，你可以看到对trt effect，
看他们有没有作用，但是貌似没有pairwise的test。如果是我的话我就用l |
|
s*****a
发帖数: 2735 | 12 试试把lifetest放到macro里？
%macro test1;
proc lifetest XXX;
strata &trt / order=internal;
....;
run;
%mend test1;
%test1;
看看这样行不行？
--------------------
不知道你的&trt是怎么来的，能把&trt产生那一步也给出吗？ |
|
|
|
g**r
发帖数: 425 | 15 DATA STEP是必须要搞熟的，除了读入文件DO LOOP什么的，
retain, first. last.； date format,这些都不在话下
MACRO很重要，很重要很重要；但MACRO学起来是很容易的，比DATA STEP
容易好多
PROC TABULATE (or PROC REPORT);要搞搞熟的，这个很重要
Haworth的 PROC TEBULATE BY Example是值得花时间读一读的
我痛恨胡乱用PROC TABULATE的PROGRAMMER，那么强的工具，经常被人
糟塌--明明三五行就解决的事情，楞写出百十行来，让人看的那个难受。
直观易懂的PROC MEANS, PROC FREQ当然不用说了
PROC GPLOT和PROC GCHART要翻一翻的，AXIS,PATTERN, LENGEND,SYMBLO什么的
临时学倒是来的及。但人家问你的时候你要知道SAS GRAPH这回事情。
如果是搞TRIAL的，PROC LIFETEST, PROC PHREG要看一看的；
否则的话，PROC REG, PROG LOGISTIC, PROC GLM翻翻 |
|
c**********e
发帖数: 2007 | 16 Use proc lifetime. It is the parametric approach.
proc reliability in SAS/QC is the same as proc lifetime
in SAS/STAT. Be careful for left censored case, one
of them requires left bound as 0 while the other requires
missing (use .).
For interval censored data, SAS only has parametric
approach. For other data, SAS has nonparametric, i.e.,
Kaplan-Meier (proc lifetest), and even more popular
semi-parametric approach. i.e., Cox regression (proc
phreg). |
|
G********t
发帖数: 334 | 17 用的KM方法，想要Product-Limit Survival Estimates的置信区间！为什么输出里面没
有？用了CONFTYPE= ALPHA= 的OPTION都没用？请高人指教！在线等！！谢谢！！ |
|
|
G********t
发帖数: 334 | 19 btw,为啥它不直接放在OUTPUT里面？难道是我SAS版本低了？ |
|
|
a********a
发帖数: 346 | 21 Do you know how to test or see the proportinality assumpiton in PROC
Lifetest?
Thanks |
|
p*****0
发帖数: 3104 | 22 看到了一些的区别 但是不知道有没有看出本质来
多谢！ |
|
|
m*******s
发帖数: 469 | 24 proc univariate
proc freq
proc glm
proc lifetest
proc phreg
proc mixed
真正弄懂玩透这几个SAS命令，药厂95%的project可以应付了。 |
|
l****g
发帖数: 304 | 25 1: How can I change the Y-axis name from default setting 'Survival
Districution Frunction' to something else I needed such as freedom from CAD
death?
2. How I can change the position of the legend? (From the top of the right
side to some other place).
Thank you so much! |
|
|
w****a
发帖数: 114 | 27 这个问题不好回答，不过你应该看看SAS online help里面的gplot部分。
其实很简单，第一个，你需要定义axis，然后在你的proc gplot statement里面用你定
义的axis；
第二个，有一个legend statement，其中有定义位置的option，你可以在这里定义位置
。 |
|
l******3
发帖数: 23 | 28 生物转到统计的，背景不好，是个作业题目，实在是做不出来了，好心人帮忙看看，先
谢谢了。
In a clinical trial of leukemia patients in remission, 12 patients were
randomized into treatment and placebp arms with 6 on each arm. The observed
waiting times to relapse (positive numbers) or censored times (negative
numbers) in the two arms, recorded in weeks are arranged in ascending order
as follows:
T: 2, 2, 3, -3, -5, 10
P: 1, 1, 2, -3, 5, -7
Using longrank test and the generalized wilconxin test by SAS proc IML.
不会用proc iml求，只会用proc lifetest. |
|
j*****e
发帖数: 182 | 29 Your instructor is asking you to apply the formula you learned in class and
write a program in SAS to perform the test.
Your can first compute the test statistic by hand. Compare it with the
output from lifetest. If it is correct, convert your hand calculation into
SAS IML.
The formula for this type of test is not as easy as it seems. You wouldn't
appreciate it until you do your own calculation. |
|
|
n*****n
发帖数: 3123 | 31 you need to re-arrange your data first (can be easily done by sas)
use your data as an example
day censorflag count
1 0 5
1 1 10
2 0 5
2 1 8
3 0 1
3 1 3
4 1 8
then use the sas statements
proc lifetest data=**** method=lt;
time day*censorflag(0);
freq count;
run; |
|
c*********n
发帖数: 87 | 32 我觉得你是要看group的差异的话，应该还需要一个group indicator，然后用starta，
这样才能test group 差异。
day censorflag count group
1 0 5 1
1 1 10 0
2 0 5 0
2 1 8 1
proc lifetest data=**** method=lt;
time day*censorflag(0);
freq count;
strata group；
test group;
run;
我要看的是 group 之间的差异，例如两个survival function是不是有显著性差异。 |
|
j*****e
发帖数: 182 | 33 If you use lifetest to check the group effect, do not use method=it. It
corresponds to lifetable, which is used for summarizing large data set and
no tests could be performed.
Use method=pl. This is also the default. |
|
c*********n
发帖数: 87 | 34 test两个group的survival function是不是一样，也就是test两个group有没有差异，
我们把两个group看成两个treatment，就是看treatment的差异。
你可以用plots=（s）画出2组的KM estimates，就比较直观了。
proc lifetest plots=(s);
....... |
|
s*******e
发帖数: 1385 | 35 已经很感谢了，parametric 和 nonparametric我都用，就是为了看看哪个更好一些。
lifetest |
|
s*******e
发帖数: 1385 | 36 我有一组survival data
（days） （days）
startpoint endpoint count status
1 3 3 1
4 6 5 0
8 10 3 0
10 12 4 1
我不知道animal具体是什么时候死或者censored的，只知道在startpoint和endpoint之
间。count是在这个区间failed或者censored的个体数，status是censor indicator。
在SAS里面用lifetest，how to specify interval for each event?
intervals statement不识别variable,而且好像都是固定的interval。还有就是在time
statement里面，应该用endpoint还是startpoint的时间？
望达人指教，包子答谢。 |
|
B******y
发帖数: 9065 | 37 PROC REPORT, MEANS, FREQ, GLM and LOGISTICS，包括换汤不换药的有UNIVARIATE，
MIXED, GENMOD等，基本上涵括了90％以上的工作。剩下的10％就为PROC LIFETEST，
PHREG，POWER，PLAN等等了。 |
|
a****a
发帖数: 3411 | 38 刚刚开始学sas
没用过ods
想用ods生成survival function,但是没有显示出来，谁能帮看看问题出在哪？
ODS GRAPHICS ON;
PROC LIFETEST DATA=mylib.survival3 PLOTS=(C S);
TIME years*status(0);
RUN;
ODS GRAPHICS OFF;
log的显示 |
|
m*********x
发帖数: 127 | 39 紧急请教一个问题：
一个医学study, 我用PROC LIFETEST生成了Survival curve，但是因为观测的时间比较
短，大部分subjects（90%左右）都没有fail,都在最后be censored, 得到的survival
curve都挤在上边一窄条地方，根本看不出阶梯的样子。老板要求我把Y坐标改一改，比
如把最小值改为60%。
我看了半天sas use's guide, 找不到可以改Y坐标的option。
请问有谁知道吗，请给点线索，多谢了！ |
|
J********J
发帖数: 571 | 40 请问哪些proc 可以算出confidence interval?
proc lifetest好像可以,还有别的吗 |
|
|
y***n
发帖数: 51 | 42 try this:
proc lifetest data='c:\alda\rearrest' plot=(lls); % some version is plot=
loglogsurv
model months*censor(1)=treat;
run; |
|
m******u
发帖数: 277 | 43 First, you need to know commonly used parametric, non-parametric and semi-
parametric methods. If you cannot remember, please go back to check your
notes.
Second, know some SAS procedures will be helpful, such as PROC LIFETEST,
PROC PHREG.
If you did a project before, know how to make a K-M curve and provide Log-
rank test is also very basic.
Hope that helps. |
|
c**d
发帖数: 104 | 44 in SAS 9.2, the default transformation type is LOGLOG .In SAS 9.1.3 the
transformation type is ASINSQRT |
|
a******p
发帖数: 414 | 45 我有个问题想请教一下这里的高手:
我在用proc lifetest 画survival plot时,
我想把strata &trt; 后面加个option: strata &trt / order=internal;
但&trt是character,所以我需要用format 和 informat.
问题是&trt是macro variable, 我不知道如何用format 和 informat.
多谢! |
|
s*****a
发帖数: 2735 | 46 有两个问题：
1. 如果不是要放在loop里，或者要生成很多同样类型的图，为什么要生成&trt。生成&
trt的目的是什么？
2. 如果只是单纯的想在最下面legend里头改变treatment label，可以把survival
curve从proc lifetest output到一个dataset，然后用proc gplot 画出来。在gplot里
头用legend statement：
legend1 across=1 frame
label=(h=2 position=bottom j=c 'Treatment label:')
value=( h=2 'One' 'Three' 'Two');
不好意思，理解不太清楚。如果你能够把你最终想要的描述一下，或者能帮你更多。 |
|
a******p
发帖数: 414 | 47 shincha, 谢谢你这么热心的给我解答。
1. 我在写一个survival macro, 类似下面的：
Macro SYNTAX
%SurvTab(data=,
stratum_list=,
time=,
censor=,
estimate=,
graphfile=,
rtffile=,
title= );
&trt是macro parameter, 类似于上面的stratum_list.
2.我没有用gplot来画图，我用的是proc lifetest + plot(s)
所以我想通过在strata statement 后加 order=internal option 来控制我的
treatment label.
不知我说清楚没有，谢谢！
成& |
|
s******8
发帖数: 102 | 48 不清楚究竟是啥问题。“问题是&trt是macro variable, 我不知道如何用format 和
informat.”
是否可以照猫画虎？
%macro test1;
data xxx;
format &trt yourmatnamehere.;
set xxx;
run;
proc lifetest XXX;
strata &trt / order=internal;
....;
run;
%mend test1;
%test1; |
|
a******p
发帖数: 414 | 49 谢谢你的回答。我也是赶鸭子上架，表述不清楚。
下面的code是我从网上找到的。我的&trt 相当于下面code的里的group,
我的问题是，&trt是macro variable, 我不知怎么向下面一样写format
proc format;
invalue $bmtifmt 'ALL' = 1 'AML-Low Risk' = 2 'AML-High Risk' = 3;
value bmtfmt 1 = 'ALL' 2 = 'AML-Low Risk' 3 = 'AML-High Risk';
run;
data Bmt2;
set sashelp.BMT(rename=(Group=G));
Group = input(input(G, $bmtifmt.), 1.);
label Group = 'Disease Group';
format Group bmtfmt.;
run;
proc LIFETEST data=Bmt2 plots=s(atrisk(outside maxlen=13)=0 to 2500 by 500);... 阅读全帖 |
|
|
