SAS Support Communities

Dear all,   I would like to assess the relationship between time of day of (protein) intake and overall (protein) intake. Therefore, I want to calculate corelations between the proportion ingested during each hour of the day and the total amount ingested over the entire day. In the end I would like to create a graph such as the one shown in the attachment. In this example, however, time is divided into categories, whereas I want to handle time as a continuous variable.   The difficulty is that I am working with repeated measures (2 observations per participant). I do not know exactly how to handle this in my syntax.   My data looks as follows. So, for each participant I have 48 records (24 hours for each day). Since I want to correlate the proportion of protein ingested during each hour of the day, I created the variable (proportion_protein), by dividing protein_hour by protein_day and multiplying it by 100. data test1; input subject:12. day:12. hours:12. protein_day:12. protein_hour:12. proportion_protein:12.; datalines; 1 1 8 99 0 0 1 1 9 99 11 11 1 1 10 99 0 0 ... ... ... 750 2 22 56 0 0 750 2 23 56 0.8 1.5 run;     I've thought about the following syntax (adapted from: http://www2.sas.com/proceedings/sugi29/198-29.pdf). However I do not know if I am correct. Do I have to add both proportion_protein and hour into the model, or is this not needed? I have the feeling that if I do not do it, SAS won't know which proportions belong to which hours of the day. ods graphics on; PROC mixed DATA=test method=reml covtest plots(maxpoints=none); Class subject day (ref=first) hour (ref=first); MODEL prot_day = proportion_protein hour proportion_protein*hour / solution ddfm=kr ; Random prot_hour /type=un subject=subject v vcorr; Repeated day / type=un subject=hour(subject) r rcorr; run; ods graphics off;   After running this syntax, log says the following: NOTE: Convergence criteria met. NOTE: PROCEDURE MIXED used (Total process time): real time 9.96 seconds cpu time 9.82 seconds     I do not exacly know if this is the outcome I needed. Also no graphics are shown. Is it even possible to translate this data into a graph as the one shown in the attachement. Or is this not possible? Can someone maybe help me with this? Many thanks in advance.   ... View more

Dear all,   I would like to calculate time differences between hour-of-wake-up and hours-of-eating (difference in minutes). For each participant I have the hour-of-wake up (e.g. 08:30; character variable $5.) and the hours at which foods are consumed (e.g. 09; character variable $3.). I've recently read a lot about handling time in SAS, however I do not (yet) succeed to calculate the difference in time. I have mainly difficulties with the difference in format statement between the two variables. Can someone please explain me how to calculate the difference?   See attachment for example dataset and the output I would like to have.   Many thanks in advance     ... View more

Hi!   I want to show the distribution of food intake that is consumed over the day, and more specifically... I want to show the distribution of both protein and energy intake (both on the y-axis) that is consumed over the day, per food consumption occasion (fco) (e.g. breakfast, lunch) and per unit of time (both on the x-axis) in a graphical way.   I would like to project energy and protein as a vertical bar, and if possible I want both variables to be in the same bar (protein as part of energy per food consumption occasion). Also, I would like to present the food consumption occasion (fco) in different colours, so it is easier to observe differences per fco.   I've already tried some gplot syntaxes, but these did not produce the desired outcome. PROC GPLOT data=test_1; PLOT sum_prot*hour=fco PLOT sum_enkcal*hour=fco; RUN;   This syntax created two separate plots for protein and energy intake (see attachment for plot energy intake) instead of one in which protein and energy are combined. Also, no bars were created (which is logical when using gplot), and therefore protein and energy intake could not be presented in the same bar.   Can someone maybe help me with this? Is it even possible to present protein and energy intake in the same bar, and if, what kind of syntax do I need for this?   Many thanks in advance.   Best Regards, Teuni ... View more