Hi @Choit and welcome to the SAS Support Communities!
PROC NLIN can estimate the parameters:
data have;
input x y;
cards;
1 0.000516
2 0.000847
3 0.001459
4 0.001939
5 0.002075
6 0.001749
;
ods output ParameterEstimates=est;
proc nlin data=have;
parms a=3 b=2 c=.01;
model y=c*pdf('gamma',x,a,b);
run;
proc print data=est noobs;
var parameter estimate;
run;
Result:
Parameter Estimate a 3.2866 b 2.1503 c 0.0165
The estimates for a and b differ slightly from your Excel values, which might be due to rounding error in variable Y. If the Y values are taken as exact values, the sum of squares is smaller for the above estimates. This is actually not maximum likelihood estimation. The Gauss-Newton method with initial values a=3, b=2, c=.01 was used to minimize the sum of squares.
Hi @Choit and welcome to the SAS Support Communities!
PROC NLIN can estimate the parameters:
data have;
input x y;
cards;
1 0.000516
2 0.000847
3 0.001459
4 0.001939
5 0.002075
6 0.001749
;
ods output ParameterEstimates=est;
proc nlin data=have;
parms a=3 b=2 c=.01;
model y=c*pdf('gamma',x,a,b);
run;
proc print data=est noobs;
var parameter estimate;
run;
Result:
Parameter Estimate a 3.2866 b 2.1503 c 0.0165
The estimates for a and b differ slightly from your Excel values, which might be due to rounding error in variable Y. If the Y values are taken as exact values, the sum of squares is smaller for the above estimates. This is actually not maximum likelihood estimation. The Gauss-Newton method with initial values a=3, b=2, c=.01 was used to minimize the sum of squares.
Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.
If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website.
Learn the difference between classical and Bayesian statistical approaches and see a few PROC examples to perform Bayesian analysis in this video.
Find more tutorials on the SAS Users YouTube channel.