Hi, I need to estimate parameters a and b as well as R-squared from data for each group (expno=1, expno=2, expno=3 etc.,) separately at once. Could someone help me to modify the below codes? Many thanks FILENAME CSV "/folders/myfolders/umar2021.csv" TERMSTR=CRLF;
PROC IMPORT DATAFILE=CSV
OUT=umar
DBMS=CSV
REPLACE;
RUN;
title 'LWR for shells';
data exp;
set umar;
if expno=1 or expno=2 or expno=3;
run;
proc nlin data=exp method=marquardt;
parms B=3 A=0.01;
model w=A*tl**B;
output out=expp p=pw r=w_residual;
run;
proc summary data=expp;
var w_residual w;
output out=stats css(w)=sstot uss(w_residual)=ssres N=N;
run;
data expp;
set stats;
rsquared=1-(ssres/sstot);
adjrsquared = 1-(1-rsquared)*(N-1) / (N- 3 -1);
run;
proc print data=expp;
run;
... View more