This code below creat a uniform (0,1), how can I modify to creat uniform (0.2, 1) data t; do i = 1 to 50; a=uniform(0); output; end; run; ... View more

Hey, How do I fixed this so that nothing is left our: PROC SGPLOT data=nn1; label Kappa_A="A, 1-R10=01, P=14/69, Sens_n=.65, Sens_o=.60, Spec_o=.60, Spec_n=.60"  Kappa_B="B, 1-R10=.90, P=14/69 Sens_n=.65, Sens_o=.60, Spec_o=.60, Spec_n=.60"  Kappa_C="C, 1-R10=.85, P=14/69 Sens_n=.65, Sens_o=.60, Spec_o=.60, Spec_n=.60" Kappa_D="D, 1-R10=.80, P=14/69 Sens_n=.65, Sens_o=.60, Spec_o=.60, Spec_n=.60" Kappa_E="E, 1-R10=.30, P=14/69 Sens_n=.65, Sens_o=.60, Spec_o=.60, Spec_n=.60" Kappa_F="F, 1-R10=.25, P=14/69 Sens_n=.65, Sens_o=.60, Spec_o=.60, Spec_n=.60" Kappa_G="G, 1-R10=.20, P=14/69 Sens_n=.65, Sens_o=.60, Spec_o=.60, Spec_n=.60" Kappa_H="H, 1-R10=.10, P=14/69 Sens_n=.65, Sens_o=.60, Spec_o=.60, Spec_n=.60"; series x=R y=Kappa_A / markers ; series x=R y=Kappa_B / markers; series x=R y=Kappa_C / markers; series x=R y=Kappa_D / markers; series x=R y=Kappa_E / markers ; series x=R y=Kappa_F / markers; series x=R y=Kappa_G / markers; series x=R y=Kappa_H / markers; xaxis label = '1-R01' grid; yaxis label = 'Kappa' grid; run; NOTE: There were 10 observations read from the data set WORK.NEW3. NOTE: The data set WORK.NN3 has 10 observations and 23 variables. NOTE: DATA statement used (Total process time):       real time           0.04 seconds       cpu time            0.01 seconds 4439  PROC SGPLOT data=nn2; 4440  label Kappa_A="A, 1-R10=01, P=14/69, Sens_n=.90, Sens_o=.60, Spec_o=.60, Spec_n=.60" 4441  Kappa_B="B, 1-R10=.90, P=14/69 Sens_n=.90, Sens_o=.60, Spec_o=.60, Spec_n=.60" 4442  Kappa_C="C, 1-R10=.85, P=14/69 Sens_n=.90, Sens_o=.60, Spec_o=.60, Spec_n=.60" 4443  Kappa_D="D, 1-R10=.80, P=14/69 Sens_n=.90, Sens_o=.60, Spec_o=.60, Spec_n=.60" 4444  Kappa_E="E, 1-R10=.30, P=14/69 Sens_n=.90, Sens_o=.60, Spec_o=.60, Spec_n=.60" 4445  Kappa_F="F, 1-R10=.25, P=14/69 Sens_n=.90, Sens_o=.60, Spec_o=.60, Spec_n=.60" 4446  Kappa_G="G, 1-R10=.20, P=14/69 Sens_n=.90, Sens_o=.60, Spec_o=.60, Spec_n=.60" 4447  Kappa_H="H, 1-R10=.10, P=14/69 Sens_n=.90, Sens_o=.60, Spec_o=.60, Spec_n=.60"; 4448  series x=R y=Kappa_A / markers ; 4449  series x=R y=Kappa_B / markers; 4450  series x=R y=Kappa_C / markers; 4451  series x=R y=Kappa_D / markers; 4452  series x=R y=Kappa_E / markers ; 4453  series x=R y=Kappa_F / markers; 4454  series x=R y=Kappa_G / markers; 4455  series x=R y=Kappa_H / markers; 4456  xaxis label = '1-R01' grid; 4457  yaxis label = 'Kappa' grid; 4458  run; NOTE: PROCEDURE SGPLOT used (Total process time):       real time           0.40 seconds       cpu time            0.11 seconds NOTE: Some graph legends have been dropped due to size constraints. Try adjusting the       MAXLEGENDAREA=, WIDTH= and HEIGHT= options in the ODS GRAPHICS statement. NOTE: There were 10 observations read from the data set WORK.NN2. ... View more

When I plot thes plots, the name on the plot suppose to be " Plot _SENSIT_A      _SENSIT_B" however I am getting " Plot sensitivity      _SENSIT_B. How do I correct it? GOPTIONS reset=all cback=white border; DATA new; INPUT _SENSIT_A _SENSIT_B SPEC ; _1MSPEC_ =1-SPEC; CARDS; 1 1 0 .98 .95 0.1 .97  .91  .2 .95 .88 .3 .87  .70 .5 .70  .4 .7 .55 .30 .8 .1 0.1 1 ; DATA anno; FUNCTION='move'; xsys='1'; ysys='1'; x=0; y=0; OUTPUT; FUNCTION='draw'; xsys='1'; ysys='1'; color='red'; x=100; y=100; OUTPUT; RUN; PROC GPLOT data=new; proc gplot data=new;    plot _SENSIT_A*_1MSPEC_  _SENSIT_B*_1MSPEC_ / overlay legend=legend1                              vref=0 to 1 by .1                              lvref=2                              haxis=axis1 hminor=1                              vaxis=axis2 vminor=2; symbol2 i=join v=dot c=red line=2; SYMBOL1 i=j v=dot c=black; AXIS1 length=5 in; AXIS2 length=5 in; LABEL _SENSIT_A ='Sensitivity'; LABEL _1MSPEC_='1 - Specificity'; run;ke on the plote ... View more

How do I compute the P-value using SAS if I I have the following; 1) Estimate of Beta 2)Standard error of estimate of Beta 3) Degree of freedom of t distribution. ... View more

I am interested in computing the simulated error for the parameter t in the regression.  SAS gives me the regression error ot t which is StdEr I wish to know if sas has an easy way to get it. THanks ... View more

Hello, I need the baseline covariates fro PHREG for each of the random sample with replacement. THis code does see to work well. The problem is with baseline covariates=Pattern_&pt out=survs_&u survival= survival;. NB: I need to run this within a macro because I have some matrices to compute. Or within a for loop. CODE. data neww; seed=-1;   alpha1 = 10;   beta1 = 10;   beta2=10;   do i = 1 to 141;   Age=40+(5)*rannor(1);    lambdaT = 0.025; *baseline hazard;    lambdaC= .03;  *heavy=0.15; *light=0.03;    er=0+(0.01)*rannor(1);     t = rand("WEIBULL", 0.75, lambdaT);  * time of event;     c = rand("WEIBULL", 1.25, lambdaC) ;* time of censoring;    monset = min(t, c);    * which came first?;       obs=(t lt c); * creating observation vaiable from censored when observeed obs=1;   pibflr= alpha1 + beta1*t + beta2*Age +er;     output;   end; run; %let pt=pt; %macro pat(patt=); %do &pt=1 %to &patt;  /*** I NEED THE BASELINE COVARIATE FOR AGE=0, So for each data below new_1, new_2, I want the first to set the covariates=Pattern_1, covariates=Pattern_2*/ data Pattern_&pt; age=0; output; run; %end; %mend pat; %MACRO DO_BRANCH(replicates=); %do u=1 %to &replicates; proc surveyselect data=neww  method = urs sampsize = 141 OUTHITS      reps=1 seed=1234100 out=new_&u; run; %pat(patt=20); ods trace on /listing; ods graphics on; ods output ParameterEstimates=CPHH_&u; proc phreg data=new_&u ; model monset*obs(0)=age; *baseline out=survs survival= survival; baseline covariates=Pattern_&pt out=survs_&u survival= survival;  /* The issue is that SAS doesnot seem to understand that I want covariates=pattern_1 when running the data set new_1 and covariates=pattern_2 for new_2*/ run; ods trace off; %end; %MEND DO_BRANCH; %DO_BRANCH(replicates=20); THIS CODE RUNS BUT NOT THE ABOVE. Below is what I want the macro to do for each case. data Pattern; Age=0; output; run; ods trace on /listing; ods graphics on; ods output ParameterEstimates=CPHH; proc phreg data=new ; model monset*obs(0)=age; baseline covariates=pattern out=survs survival= survival  ; run; ods trace off; ... View more