Hi, Let's say I have a dataset of 1000 observation. I want to construct a model using the first 900 observation, and to build predictions for the 901th -1000th observation-in orde to test my model and compute the residual. How do I get predictions to the rest of my dataset? Thanks! ... View more

Hi all, I'm running several stepwise regression and I'm trying to find out how to 1) output only the chosen final model of a stepwise regression. 2) efficiently get the list of the chosen variables of the final regression, how do you recomend to do it? My code: proc reg data=&data outest=&out tableout ; model &Var.=v1-v365 a1-a12 t1-t6 /AIC selection=stepwise SLENTRY=0.05 ; OUTPUT out=&out p=predicted r=residuals STUDENT=STUDENT_residuals; where station=&stat.; run; Thanks in advance, Tal ... View more

Hi, I couldn't find an appropreate forum so... when I want to put a line or a few lines in a note ,I select the line and then press cnrl+? - this adds the /**/. for some reason it no longer works. anybody knows how to restore this option? Thanks! ... View more

Hi , I'm trying to write a macro that will include a "call symput" part. I have a data set of 1 observation , that has the following variables: intercept, &var1 and &var2. I want to assign 3 macro variables: intercept21 var1 var2 that will store these values. How can I do it? Thanks! %macro reg_fill(indat,var1,var2,stat_id); data _null_; set &indat; call symput ("Intercept&stat_id.",Intercept); call symput ("&var1.",&var1); call symput ("&var2.",&var2); run; %put &var1. &var2.; %mend reg_fill; %reg_fill(data,O3,NOx,21); ... View more

Hi! Here's my code- for the 3 values of stat_id I get 3 overlay plots. When I run this i get an automatic legend of the 3 colors of stat_id. I can cancel that automatic legend by creating my own legend, BUT how can I have a graph with no legends at all? Thanks, Tal proc gplot gout=hh.no_legend data=hh.jan_03_excel; plot O3*d=stat_id /vaxis = axis2 haxis=axis1 /*LEGEND=LEGEND1*/ name="O3 January 2003.gif"; where day in (5,6,7,8,9,10,11,12); run; quit; ... View more

Here's my code. I specified Justify and Height options, but I get a centerd title in a fix height (I tried to change the height, but it doesn't affect the graph) What's wrong with my syntax? Thanks in advance! Tal goptions reset=all; legend1 label=('Station') value=(height=1 justify=left 'A' 'B' 'C' )shape=symbol(3,0.5) position=(top right inside ) mode=protect down=3 across=1 frame; title '(a)' JUSTIFY=left HEIGHT=3 ; symbol1 i =join v=dot c=green h=0.1 cm; symbol2 i =join v=dot c=blue h=0.1 cm; symbol3 i =join v=dot c=red h=0.1 cm; axis1 order=('5jan03:00:00'dt to '12jan03:00:00'dt by day) label=(f=swiss h=0.5 cm 'Date and Time') minor=(number=24); axis2 order = (0 to 60 by 10) label=(f=swiss h=0.5 cm a=90 'O3 [ppb]'); proc gplot gout=w data=hh.jan_03_excel; plot O3*d=stat_id /vaxis = axis2 haxis=axis1 LEGEND=LEGEND1 /*legend=legend2 */ name="O3 January 2003.gif"; where day in (5,6,7,8,9,10,11,12); run; quit; ... View more