topic Re: ploting data in SAS Programming

ploting data

always-good — Tue, 29 Nov 2022 20:02:05 GMT

Hello SAS community!

I'm working on fitting the lactation curve of dairy cows according to the calving season (1=winter, 2= spring and 3= summer) using the model of Wood for fitting.

My question here is that I want to draw the three curves together as one curve according to the three seasons and not each curve separately.

Can anyone please correct the code for me. I really need help!

Here is the data and the tried code:

data curves;

input season days Milk;

datalines;

1 22 7.63

1 54 7.25

1 85 8.40

1 113 8.10

1 142 5.75

1 174 6.00

1 208 5.00

1 230 4.00

1 263

1 299

2 17 8.54

2 46 8.71

2 75 7.27

2 107 5.46

2 140 5.82

2 167 6.22

2 187 6.60

2 219 3.00

2 258

2 285

3 16 7.00

3 50 9.00

3 84 7.00

3 107 9.50

3 140 8.00

3 163

;

proc nlin data = curves plots = fit;

parms A = 15 B = 0.19 C = -0.0012 ;

bounds A B C > 0;

by season;

model PL = A * days **b * exp(-C*days);

output out = Fit predicted = Pred;

by season;

run;

title1' lactation curve ';

proc gplot data=fit; plot pred*days milk*days/overlay;

by season;

run;

Thank you, the community!

Re: ploting data

PaigeMiller — Tue, 29 Nov 2022 20:05:31 GMT

Please look at the log for this code and fix the error. (That could be why you can't get a plot; or it may be there are additional errors)

Re: ploting data

always-good — Tue, 29 Nov 2022 20:12:36 GMT

Thank you for your reply!
The code is running well, but that's not exactly the result I'm asking for.
I want to plot three curves in the same curve and I'm not familiar with that.

Re: ploting data

PaigeMiller — Tue, 29 Nov 2022 20:31:30 GMT

Maybe you are running different code than what you showed above? When I copy and paste that code into SAS, there is an error in PROC NLIN.

Re: ploting data

always-good — Tue, 29 Nov 2022 20:58:00 GMT

Please try this code. It is running:
data curves;
input saison Jours PL;
datalines;
1 22 7.63
1 54 7.25
1 85 8.40
1 113 8.10
1 142 5.75
1 174 6.00
1 208 5.00
1 230 4.00
1 263
1 299
2 17 8.54
2 46 8.71
2 75 7.27
2 107 5.46
2 140 5.82
2 167 6.22
2 187 6.60
2 219 3.00
2 258
2 285
3 16 7.00
3 50 9.00
3 84 7.00
3 107 9.50
3 140 8.00
3 163
;
proc nlin data = curves plots = fit;
parms A = 15 B = 0.19 C = -0.0012 ;
bounds A B C > 0;
by saison;
model PL = A * Jours **b * exp(-C*Jours);
output out = Fit predicted = Pred ;
by saison;
run;
title1' courbe de lactation ';
proc gplot data=fit;
plot pred*Jours PL*Jours/overlay;
by saison;
run;

Re: ploting data

PaigeMiller — Tue, 29 Nov 2022 22:22:49 GMT

Thank you, your code now works. However, I don't really understand this part:

I want to draw the three curves together as one curve according to the three seasons and not each curve separately.

Does this mean you want values 1 2 and 3 on the horizontal axis? If so, what should go on the vertical axis? Or do you want something else? I really need a more clear explanation.

Re: ploting data

always-good — Tue, 29 Nov 2022 22:42:58 GMT

Yes this is what I mean, values 1,2 and 3 on the horizontal axis and the same vertical axis for all curves like in the code, I means "Jours".
I want to fit the curve of Milk production "PL" (horizontal axis) for the three seasons 1,2 and 3, according to the days "Jours" (vertical axis).

Re: ploting data

PaigeMiller — Tue, 29 Nov 2022 22:54:09 GMT

@always-good wrote:
Yes this is what I mean, values 1,2 and 3 on the horizontal axis and the same vertical axis for all curves like in the code, I means "Jours".
I want to fit the curve of Milk production "PL" (horizontal axis) for the three seasons 1,2 and 3, according to the days "Jours" (vertical axis).

I'm still not grasping what you want. What "curves"? What is the value on the Y-axis? Can you show me a similar plot that is on the internet somewhere?

Re: ploting data

always-good — Tue, 29 Nov 2022 23:06:37 GMT

Please check this out

Re: ploting data

PaigeMiller — Wed, 30 Nov 2022 11:24:10 GMT

Please include a screen capture of this plot by clicking on the "Insert Photos" icon. I refuse to download files; and furthermore Microsoft Office documents can be a security threat.

Re: ploting data

always-good — Wed, 30 Nov 2022 11:57:09 GMT

lactation curves according to different season of calving

Re: ploting data

PaigeMiller — Wed, 30 Nov 2022 12:24:48 GMT

This does not look anything like what you said earlier "Yes this is what I mean, values 1,2 and 3 on the horizontal axis"

This doesn't seem to match your original description of the plot "My question here is that I want to draw the three curves together as one curve", this plot shows separate curves and does not draw the three curves together as one curve.

Re: ploting data

always-good — Wed, 30 Nov 2022 12:32:12 GMT

Please if you understand what I mean, you can help me! I think it is the same thing meant:
The "PL" for the three seasons 1,2 and 3 are in the horizontal axis. I think it is clear with the picture.
what is the thing that you didn't understand? we can find a solution?

Re: ploting data

PaigeMiller — Wed, 30 Nov 2022 12:38:49 GMT

Not clear. The plot shows a horizontal axis going from 5 to 305. Do you want a horizontal axis that has only values 1, 2 and 3? (Because that's what you said earlier) And do you want ""My question here is that I want to draw the three curves together as one curve", or do you want three curves?

Re: ploting data

always-good — Wed, 30 Nov 2022 13:14:46 GMT

I'm sorry you are right.
the PL values for season1,2 and 3 are in the vertical axis.

Re: ploting data

PaigeMiller — Wed, 30 Nov 2022 13:22:49 GMT

The example you show does not have 1,2,3 on the vertical axis either.

It's time for a re-write. Please start from scratch, and write what you want in the plot in detail, correctly, clearly, completely. Put some effort into this, make sure everything is clear, complete, correct. Please explain how the data in the output data set named FIT corresponds to the plot that you show. I want this all in one place so I don't have to keep scrolling up and down to read the various parts of your message and see your example, and so I don't have to ask five more questions, so I don't even have to ask one more question.

Re: ploting data

Rick_SAS — Wed, 30 Nov 2022 13:55:16 GMT

I suggest you use the "missing value trick" to create an evenly spaced set of Jours for scoring. Then run the model and use PROC SGPLOT to overlay the data and curves. You do this by using the GROUP=saison option:

/* to improve the visualization, add a sequence of evenly spaced
values to the data. This is the scoring set. */
data Score;
Score = 1;
do saison = 1 to 3;
   do Jours = 15 to 250 by 5;
      output;
   end;
end;
run;

data All;
set curves Score;
run;

proc sort data=All;
by saison Jours;
run;

proc nlin data = All plots = fit;
parms A = 15 B = 0.19 C = -0.0012 ;
bounds A B C > 0;
by saison;
model PL = A * Jours **b * exp(-C*Jours);
output out = Fit predicted = Pred ;
run;

title1' courbe de lactation ';
proc sgplot data=fit;
scatter y=PL x=Jours / group=saison;
series y=pred x=Jours / group=saison;
run;

Re: ploting data

always-good — Wed, 30 Nov 2022 14:39:11 GMT

Dear Rick,
This is exactly what I'm talking about. This is what I'm asking for!
The code is running well.
Thanks a lot, and I'm so grateful!
Best regards.