turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Find a Community

- Home
- /
- SAS Programming
- /
- General Programming
- /
- [need help] Split plot design with repeated measur...

Topic Options

- Subscribe to RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Printer Friendly Page

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

12-30-2013 03:00 PM

Hi all,

I am a graduate student at Auburn university. My experiment design is kind of complicated and I have no idea how to do the statistics. Please help me.

Here I attached a picture of my experiment design. The field experiment was a split plot design: the main plot is different Ca treatments, and the sub-plot is two peanut cultivars. The experiment was repeated 4 times (4 blocks). Peanut samples were collected from each sub-plot three times: 4 weeks before harvest, 2 weeks before harvest and at harvest. Each peanut sample was further divided into 5 maturity classes: yellow 1, yellow 2, orange, brown and black. The response variable is Ca content in peanut seed.

So my experiment includes 4 blocks (Block), 4 treatments (Main plot), 2 peanut cultivars (Sub-plot), 3 sampling times (timing) and 5 maturity classes (M_class). I'd like to know whether there is a difference between sampling times, peanut cultivars, treatments and maturity classes. I know I can use PROC MIXED to do this but I don't know how to write the SAS program. Anyone can help me? Thank you!

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

12-31-2013 08:04 AM

Hi D,

Do you have access to Littell et al.'s *SAS System for Mixed Models, 2nd. ed.*? I'm sure somebody in the agronomy dept. has a copy. Anyway, look at Chapter 16 - Case Studies, for some examples that approach your design (essentially, everything except the repeated measures part).

Here is my UNTESTED code:

proc mixed data=yourdata;

class block trt cult time mclass;

model trt|cult|mclass|time/ddfm=kenwardrogers;

random intercept/subject=block;

repeated time/subject=block(mclass*trt*cult) type=csh; /* A lot of possible choices for covariance type. If you suspect time dependence, you may want to try ar(1), arh(1) or ante(1) */

lsmeans trt|cult|mclass|time/diff; /* You'll probably want to look at sliced effects here, but for now all possible comparisons */

run;

Steve Denham

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

12-31-2013 12:57 PM

Thank you very much !

I actually have the book in hand, and I will read the chapter you recommended.

Rui

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-14-2014 06:04 PM

Hi

I have two questions:

1.

Can you please explain why the subject is block(mclass*trt*cult) in the REPEATED statement?

2.

I don't quite understand the meaning of "random intercept/subject=block". does this treat all the random effects associate with block?

Thank you!

- Mark as New
- Bookmark
- Subscribe
- Subscribe to RSS Feed
- Highlight
- Email to a Friend
- Report Inappropriate Content

01-15-2014 02:12 PM

1. The subject is set up so that there is exactly one observation per timepoint per all of the effects in the model. If you had something other than that, you would get an error message referring to an infinite likelihood.

2. The "random intercept/subject=block" models only a block effect. If you have, or suspect, random effects due to fixed effect by block interactions, those can be modeled by adding the fixed effects that you expect (or need) to the statement BEFORE the slash.

Steve Denham