Thanks PG Stats, I run like that: FILENAME CSV "/folders/myfolders/Sir nadir.csv" TERMSTR=CRLF; PROC IMPORT DATAFILE=CSV OUT=geumar2 DBMS=CSV REPLACE; RUN; title 'geumar'; data group; set geumar2; run; data group; set geumar2; rename x=t y=S_t; label y="S(t)" x="t"; run; proc nlin data=exp; parms slope_0 = 1.0 intercept_1 = 1.75 intercept_2 = 0.75 exponent = 1.0 ; slope = slope_0*fish_size**exponent; group_1 = (group = 1); group_2 = (group = 2); model S_t = intercept_1*group_1 + intercept_2*group_2 - slope*t; output out=fishPred predicted=p_S_t; run; it doesn't give any error but it doesn't converge. WARNING: Step size shows no improvement. WARNING: PROC NLIN failed to converge. I am sharing a file which has total 6 groups, I edited your constricted proc to include others group too but the mass exponent is -48. According to my data analysis done with hands, it should be near to 0.55.
... View more