Re: Calculate overall differences

KatieBachert23 · Posted 07-31-2020 12:17 PM

I am working on this project and need to calculate "overall differences in time spent between periods" with the data I have. I have imported my excel workbook into SAS studio and was wondering what is the next step I should take. Here is a screenshot of what I am working with, there are 5 different time periods in the data set.

PaigeMiller · Posted 07-31-2020 12:26 PM

Are you talking about differences of the means among the 5 time periods? Or something else? What does the word "overall" indicate here?

--
Paige Miller

KatieBachert23 · Posted 08-04-2020 11:41 AM

I am referring to the overall difference of the mean from each of the behaviors (pecking, preening, dustbathing, and shaking) based on the 5 time periods. I want to see if there were differences in times spent per behavior per bird based on the time periods. And then I want to see if there were differences in times spent per behavior per flock as well.

Thank you

PaigeMiller · Posted 08-04-2020 12:51 PM

By bird:

proc glm data=have;
    class time_period bird;
    model pecking preening dustbathing shaking = time_period bird time_period*bird;
    lsmeans time_period bird time_period*bird/lines tdiff;
run;
quit;

By flock and bird

proc glm data=have;
    class time_period bird flock;
    model pecking preening dustbathing shaking = time_period flock bird(flock) time_period*flock time_period*bird(flock);
    lsmeans time_period flock bird(flock) time_period*flock time_period*bird(flock)/lines tdiff;
run;
quit;

Since the second analysis accounts for the effects of both bird and flock, I would simply do that analysis and not the analysis which doesn't use flock.

In either case, if you have an awful lot of birds, maybe this isn't the best way to go.

--
Paige Miller

KatieBachert23 · Posted 08-04-2020 02:05 PM

I put the code in and this was the outcome (see attachments). I see that I need to fix the data= part. How do I get the data I want to the data=_____? I then did proc glm data=mite_sen (which is the table I showed in my first post) and says that " ERROR: File WORK.MITE_SEN.DATA does not exist"

Thank you!

KatieBachert23 · Posted 08-04-2020 03:01 PM

I figured out the data= and fixed the wording so the variable matched up. When I ran the code this came up...do you know what this means?

PaigeMiller · Posted 08-04-2020 03:29 PM

It's trying to do some plots, and it is having problems. How many different birds and flocks do you have?

You might want to try turning off the ODS GRAPHICS before you run PROC GLM, although certainly some of the plots can be helpful. Use this command

ods graphics off;

--
Paige Miller

KatieBachert23 · Posted 08-04-2020 03:59 PM

There are 4 flocks total and each flock has 12 birds

KatieBachert23 · Posted 08-06-2020 08:56 AM

I did the ods graphic off; statement and it still won't run

PaigeMiller · Posted 08-06-2020 09:51 AM

What do you mean by "it still won't run"? What do you see? Can you show us the ENTIRE log for PROC GLM, including the code, NOTEs, WARNINGs and ERRORs? Please copy the log as text and then paste it into the window that appears when you click on the </> icon here. Do not show us the log as a screen capture.

--
Paige Miller

KatieBachert23 · Posted 08-06-2020 02:02 PM

Code:

options validvarname=v7;
libname mite_sen xlsx "/folders/myfolders/myfolders/UC Davis Thesis/Mite Sensor.xlsx";
run;

proc glm data=mite_sen.percentofday;
class time_period bird flock;
model pecking_s preening_s dustbathing_s shaking_s = time_period flock bird(flock) time_period*flock time_period*bird(flock);
lsmeans time_period flock bird(flock) time_period*flock time_period*bird(flock)/lines tdiff;
run;
quit;

Log:

1 OPTIONS NONOTES NOSTIMER NOSOURCE NOSYNTAXCHECK;

72

73 proc glm data=mite_sen.percentofday;

NOTE: Variable Name Change. Time Period -> Time_Period

NOTE: Variable Name Change. Pecking Hz -> Pecking_Hz

NOTE: Variable Name Change. Preening Hz -> Preening_Hz

NOTE: Variable Name Change. Dustbathing Hz -> Dustbathing_Hz

NOTE: Variable Name Change. Shaking Hz -> Shaking_Hz

NOTE: Variable Name Change. Pecking s -> Pecking_s

NOTE: Variable Name Change. Preening s -> Preening_s

NOTE: Variable Name Change. Dustbathing s -> Dustbathing_s

NOTE: Variable Name Change. Shaking s -> Shaking_s

NOTE: Variable Name Change. Pecking PercentOfDay -> Pecking_PercentOfDay

NOTE: Variable Name Change. Preening PercentOfDay -> Preening_PercentOfDay

NOTE: Variable Name Change. Dustbathing PercentOfDay -> Dustbathing_PercentOfDay

NOTE: Variable Name Change. Shaking PercentOfDay -> Shaking_PercentOfDay

74 class time_period bird flock;

75 model pecking_s preening_s dustbathing_s shaking_s = time_period flock bird(flock) time_period*flock

75 ! time_period*bird(flock);

76 lsmeans time_period flock bird(flock) time_period*flock time_period*bird(flock)/lines tdiff;

77 run;

NOTE: The import data set has 1896 observations and 16 variables.