a GYM is a new variable created by putting (Group || Test Year || Test Month), Test year  (Y) and Test Month (M) were taken from TD variable by Y = year (TD) and M = month (TD)   ... View more

Yes I tray several joints, but I can't find any solution. can you make me your proposal in SAS code!!   ... View more

Next code is based on @Patrick code with sligh changes in steps 2 and 3 I hope this will fit to what you want: data factors; input Effect $ fact1 fact2 fact3 fact4 fact5; datalines; HP11 0.1 -0.2 0.1 0.1 -0.3 HP12 0.2 -0.3 0.3 0.4 0.1 HP13 0.1 0.2 0.3 0.4 0.5 HP21 0.3 -0.4 0.7 0.2 -0.3 HP22 0.1 -0.7 0.2 0.2 -0.2 HP23 0.2 0.2 0.2 0.2 0.5 HP32 0.1 0.5 0.1 0.1 0.2 HP34 0.1 0.3 0.1 0.1 0.3 HP35 0.7 -0.1 0.6 0.3 0.1 SP11 0.3 -0.4 0.7 0.7 -0.3 SP12 0.1 0.2 0.2 0.2 -0.1 SP13 0.9 -0.1 0.5 0.1 0.3 SP21 0.4 -0.1 0.7 0.1 -0.3 SP22 0.5 -0.2 0.1 0.7 -0.3 SP23 0.1 -0.2 0.1 0.3 -0.2 SP32 0.6 -0.2 0.5 0.3 -0.1 SP34 0.7 -0.4 0.5 0.2 0.3 SP35 0.9 -0.2 0.1 0.1 0.1 ; run; data horizons; retain n; drop n; infile datalines truncover; input time val1 val2 val3 val4; n+1; if time=99 then call symput('_nobs',left(n-1)); else output; datalines; 1 0.01 0.02 0.05 0.08 2 0.02 0.03 0.02 0.08 3 0.07 0.02 0.03 0.03 4 0.02 0.02 0.04 0.08 5 0.04 0.07 0.01 0.08 6 0.01 0.07 0.01 0.02 7 0.02 0.05 0.02 0.07 8 0.06 0.02 0.02 0.04 9 0.01 0.06 0.02 0.08 10 0.03 0.06 0.06 0.08 11 0.02 0.04 0.05 0.08 12 0.02 0.01 0.05 0.08 13 0.07 0.01 0.04 0.08 14 0.09 0.01 0.03 0.02 15 0.02 0.06 0.02 0.08 16 0.01 0.07 0.03 0.04 17 0.02 0.09 0.06 0.02 18 0.02 0.02 0.02 0.02 19 0.08 0.01 0.07 0.03 20 0.06 0.03 0.02 0.02 21 0.01 0.02 0.03 0.05 22 0.06 0.04 0.01 0.06 23 0.02 0.04 0.02 0.05 24 0.04 0.01 0.05 0.09 25 0.01 0.01 0.03 0.01 99 ; run; %put No of Values = &_nobs; data want; set factors; retain tar1-tar&_nobs; array tar {&_nobs} tar1-tar&_nobs; do _point=1 to &_nobs; set horizons point=_point; tar(time) = fact1+ (fact2* val1)+(fact3* val2)+(fact4*val3)+ (fact5*val4); end; drop time val1-val4; output; run; ... View more

You have one bar, not two as in code provided. I believe one bar is OR not the concatenare operator.      ... View more

This may get you started data parm ; input id a b c month; datalines; 1 5 0.2 0.003 1 2 8 0.5 0.01 2 ; run; data plot; set parm; plotstart= mdy(month,1,1960); do t=5 to 305; Plotdate= plotstart+t; P=a*(t**b)*exp(-c*t); if t=5 then Ny=p; else NY=.; output; end; /* to use the BREAK option in series plot so each "months" series isn't connected to the other we need a missing Y value*/ plotdate=plotdate+1; p=.; output; run; ods graphics on; proc sgplot data=plot; series x=plotdate y=p/break name='L' ; needle x=plotdate y=ny/lineattrs=(pattern=4) name='N'; format Plotdate monname3.; keylegend 'L'; run; There are additional attribute data sets to set things like line colors associated with values of variables but you'll have to investigate those on your own.   ... View more