About desireatem

desireatem · ‎02-11-2013

I wish to have the survival probability of both censored and events as Stata does

desireatem · ‎02-11-2013

The following code gives me the estimated survival function for just the events, How can I get the baseline survival for both events and censored observation. I need the estimated baseline survival' proc phreg data=ph; model zz*status(0)=x1/RL TIES=EFRON; baseline out=survs survival=s; run;

desireatem · ‎02-10-2013

Thanks a lot Tom

desireatem · ‎02-10-2013

How do I delete the rows and column that are in yellow bolow. I want table 1 to look like table 2 below. proc means data=para noprint; class Variable; var estimate; output out=two mean(estimate)=col1 std(estimate)=col2 mean(StdErr)=col3; run; Table 1 Obs Variable _TYPE_ _FREQ_ col1 col2 col3 1 0 40 5.91033 3.42236 0.92869 2 Intercept 1 20 2.61495 0.08503 0.13475 3 TIMES 1 20 9.20571 1.08266 1.72264 Table 2 Obs Variable _ FREQ_ col1 col2 col3 2 Intercept 20 2.61495 0.08503 0.13475 3 TIMES 20 9.20571 1.08266 1.72264

desireatem · ‎02-07-2013

Please, How do I fix this; proc iml; use Complete; read all var _all_ into DM; close; frequency =DM[,5]; y =DM[,4]; status=DM[,3]; Censored= DM[,2]; time = DM[,1]; use GIB; read all var _all_ into Col; timeS=time; locC=loc(censored=1); matrix=J(1,6,-99); do im=1 to 25; do j= 1 to frequency[1]; A=RANUNI(-1); i=max(loc(Col[,j]-A>0)); timeS[locC ] = time+ (time[i+ 1]-time)*((Col[i,j]-A)/(Col[i,j]+col[i+1,j])); end; sim=repeat(im,nrow(time),1); temp=status||censored||time||Y||timeS||sim; Matrix=matrix//temp; end; matrix=matrix[2:nrow(matrix),]; create impute_t var{status censored time Y timeS sim}; append from matrix; quit; 843 proc iml; NOTE: IML Ready 844 use Complete; 845 read all var _all_ into DM; 846 close; NOTE: Closing WORK.COMPLETE 847 frequency =DM[,5]; 847! y =DM[,4]; 847! status=DM[,3]; 847! Censored= DM[,2]; 847! time = DM[,1]; 848 use GIB; 849 read all var _all_ into Col; 850 timeS=time; 851 locC=loc(censored=1); 851! matrix=J(1,6,-99); 852 do im=1 to 25; 853 854 do j= 1 to frequency[1]; 855 A=RANUNI(-1); 856 i=max(loc(Col[,j]-A>0)); 857 timeS[locC ] = time+ (time[i+1]-time)*((Col[i,j]-A)/(Col[i,j]+col[i+1,j])); 858 end; 858! sim=repeat(im,nrow(time),1); 859 temp=status||censored||time||Y||timeS||sim; 860 Matrix=matrix//temp; 861 end; ERROR: (execution) Invalid subscript or subscript out of range. operation : [ at line 857 column 35

desireatem · ‎02-07-2013

I fixed it; proc means data=para noprint; class Variable; var estimate; output out=two mean(estimate)=col1 std(estimate)=col2 mean(StdErr)=col3; run; The SAS System 0 50 6.7 4.45480 1.11542 2 Intercept 1 25 2.44651 0.09489 0.30363 3 TIMES 1 25 10.99305 1.57000 1.92721 1 0 50 6.71978 4.45480 1.11542 2 Intercept 1 25 2.44651 0.09489 0.30363 3 TIMES 1 25 10.99305 1.57000 1.92721

desireatem · ‎02-07-2013

This does not work; I need the standarddeviation of estimate; proc reg data=impute_t; by sim; model y=times; ods output ParameterEstimates=Para; run; proc means data=para noprint; class Variable; var estimate; output out=two mean(estimate)=col1 standarddeviation(estimate)=col2 mean(StdErr)=col3; run; proc print data=two; run; /*proc reg da

desireatem · ‎02-07-2013

Please help, Hi I need both the mean and standard deviation of the estimate mean(estimate) proc reg data=impute_t; by sim; model y=times; ods output ParameterEstimates=Para; run; proc means data=para noprint; class Variable; var estimate; output out=two mean(estimate)=col1 mean(StdErr)=col2; run;

desireatem · ‎02-05-2013

THanks Hutch, M = int(frequency[1]) works well. Programis good to go. Your help has been invaluable.

desireatem · ‎02-05-2013

Thanks , a lot, please do you know how I could set M to the (integer) number of columns that you need to read in. I do not know how to do that, googled but no helpful answer

desireatem · ‎02-05-2013

Thanks foryour help Hutch, I just started programming in imland having issues. I did what you said but still got problem; prociml; use Complete; read all var _all_ into DM; frequency = DM[,4]; M=frequency[1]; fixed_names = {"frequency""Censored" "time""y"}; cols = "COL1" "COL"+strip(char(M))); read all var(fixed_names || cols) into testA; close; print testA; 2382 proc iml; NOTE: IML Ready 2383 use Complete; 2384 read all var _all_ into DM; 2385 frequency = DM[,4]; 2386 M=frequency[1]; 2387 fixed_names = {"frequency" "Censored" "time" "y"}; 2388 cols = "COL1": ("COL"+strip(char(M))); ERROR: (execution) Invalid operand to operation. operation : : at line 2388 column 14 operands : *LIT1004, _TEM1003 *LIT1004 1 row 1 col (character, size 4) COL1 _TEM1003 1 row 1 col (character, size 15)

desireatem · ‎02-04-2013

Thanks foryou help. just started using proc iml on last Thursday. This is mymodification base on your statement; proc iml; use Complete; read all var _all_ into DM; M=frequency[1]; fixed_names = {"frequency""Censored" "time""y"}; cols = "COL1": ("COL"+strip(char(M))); read all var(fixed_names || cols) into testA; close; print testA; LOG FILE; 261 proc iml; NOTE: IML Ready 262 use Complete; 263 read all var _all_ into DM; 264 M=frequency[1]; ERROR: (execution) Matrix has not been set to a value. operation : [ at line 264 column 18 operands : frequency, *LIT1001 Frequency 0 row 0 col (numeric) 　 *LIT1001 1 row 1 col (numeric) 1 statement : ASSIGN at line 264 column 7 265 fixed_names = {"frequency" "Censored" "time" "y"}; 266 cols = "COL1" : ("COL"+strip(char(M))); ERROR: (execution) Matrix has not been set to a value. operation : CHAR at line 266 column 40 operands : M M 0 row 0 col (type ?, size 0) 　 statement : ASSIGN at line 266 column 7 267 read all var(fixed_names || cols) into testA; 268 close; NOTE: Closing WORK.COMPLETE 269 270 print testA;

desireatem · ‎02-04-2013

I am tried what you said this is what I did; proc iml; use Complete; read all var _all_ into DM; close; fixed_names = {"frequency" "Censored" "time" "y"}; M=frequency[1]; cols = "COL1" "COL"+strip(char(M))); read all var(fixed_names || cols); print ; This the log file; 1931 proc iml; NOTE: IML Ready 1932 use Complete; 1933 read all var _all_ into DM; 1934 close; NOTE: Closing WORK.COMPLETE 1935 fixed_names = {"frequency" "Censored" "time" "y"}; 1936 M=frequency[1]; ERROR: (execution) Matrix has not been set to a value. operation : [ at line 1936 column 16 operands : frequency, *LIT1002 Frequency 0 row 0 col (numeric) 　 *LIT1002 1 row 1 col (numeric) 1 statement : ASSIGN at line 1936 column 5 1937 cols = "COL1":("COL"+strip(char(M))); ERROR: (execution) Matrix has not been set to a value. operation : CHAR at line 1937 column 36 operands : M M 0 row 0 col (type ?, size 0) 　 statement : ASSIGN at line 1937 column 5 1938 read all var(fixed_names || cols); ERROR: No data set is currently open for input. statement : READ at line 1938 column 5 1939 print ; - 22 76 ERROR 22-322: Syntax error, expecting one of the following: a name, a quoted string, a numeric constant, a datetime constant, a missing value, ',', /. ERROR 76-322: Syntax error, statement will be ignored.

desireatem · ‎02-03-2013

The value for M=Frequency[1]. That is M is always equal to M=Frequency[1], the first element in frequencybut this m will vary from in each simulation. The column are name COL1-COLM, where M=M=Frequency[1]. I know there names, I just want a generic read statement that will accomodata COL1-COL4 in one matrix and COL1-COL6 in anothere proc iml; use Complete; read all var _all_ into DM; close; frequency =DM[,4]; y =DM[,3]; Censored= DM[,2]; time = DM[,1]; Cols=DM(n,m); /*Survivalival is second column*/ /*time is first column*/ n = nrow(DM); M=frequency[1]; Completed=J(nrow(survival1),M+4,1); do j=1 to M; proc iml; use Complete; read all var _all_ into DM; close; frequency =DM[,4]; y =DM[,3]; Censored= DM[,2]; time = DM[,1]; Cols=DM(n,m); /*Survivalival is second column*/ /*time is first column*/ n = nrow(DM); M=frequency[1]; Completed=J(nrow(survival1),M+4,1); do j=1 to M; read all var{COLj}

desireatem · ‎02-03-2013

I wish create a matrix C from matrix A and B by putting them side by side. For example using the following will place A to the right of B; Data C; Set A; Set B; Run; Can I have create GIB from BIG[colname={COL1}] and kaplan var{time Survival survival1}? Just one to place GIB to the of Kaplan. The number of rows wont change since kaplan GIB have the same row but column=Column (GIB)+Column(Kaplan)

Online Status	Offline
Date Last Visited	‎07-27-2024 12:54 AM

Re: Compute length of stay from 17 character date

Re: Compute length of stay from 17 character date

Re: Compute length of stay from 17 character date

Compute length of stay from 17 character date

Re: Checking if two points agree or disagree in data

Re: Checking if two points agree or disagree in data

Re: Checking if two points agree or disagree in data

Re: Checking if two points agree or disagree in data

Re: Checking if two points agree or disagree in data

Checking if two points agree or disagree in data

Re: Checking if two points agree or disagree in data

Re: Coefficient plots

Re: Plot of sample size versus year

Re: Manipulation with data and time

Re: Manipulation with data and time

Re: SAS Macro

Re: Transforming data long to wide, issue with code

Re: Using Substr function

Re: ERROR: An exception has been encountered. What type of error is t...

Re: How do I automatically clear log and output

Re: baseline Survival function for both events and censored observatio...

baseline Survival function for both events and censored observation

Re: how to delete row and column from the output

how to delete row and column from the output

ERROR: (execution) Invalid subscript or subscript out of range.

Re: I need the mean and standard deviation

Re: I need the mean and standard deviation

I need the mean and standard deviation

Re: How to name an unknown dimention column

Re: How to name an unknown dimention column

Re: How to name an unknown dimention column

Re: How to name an unknown dimention column

Re: How to name an unknown dimention column

Re: How to name an unknown dimention column

Create a matrix from two matrices in proc iml