|
Found out alone:
data t;
input id $ week y x;
datalines;
id1 1 165 80
id1 5 168 82
id1 8 169 83
id2 2 155 60
id2 3 155 60.5
id2 7 156 61
;
/* linear regression: y = a * x + b */
proc reg data = t outest = params outsscp = sscp noprint;
model y = x / selection=rsquare;
by id;
*plot y * x / name='linReg';
Title "Linear regression";
output out = linReg predicted = yhat;
run; quit;
proc sql;
create table final as
select
p.id,
p.x as a,
p.Intercept as b,
p._RSQ_ as r2
from
params p;
quit;
proc sql;
create table final2 as
select
f.*,
s.Intercept as N
from
final f, sscp s
where
f.id = s.id
and s._TYPE_ = 'N';
quit;
| 
