How to fit my data to Gamma, Weibull and Lognomal distributions?

[t75wez1]

For the sake of simplicity, I just made another data set and append them together with 62 observations.

By using the code below, I get 62 predicted values with duplication by day.

Is that right?

data customer2;set customers;
count=round(count+rand("weibull", 7,3.5),1.);
run;

data customers;set customers customer2;
run;

proc sort;by day;
run;

proc gampl data=customers;
model count=spline(day) / dist=POISSON; ******Final Pick***************;
id _all_;
output out=preds pred=predicted;
run;

[StatDave · Accepted solution]

Yeah, I guess that'll work for adding one other set of data.

[t75wez1]

Many thanks for your helps. You made my day!