Programming the statistical procedures from SAS

Setting up the random effects model in Proc mixed

Reply
Contributor
Posts: 40

Setting up the random effects model in Proc mixed

My experiment design is as follows.

Each sample in the experiment is run on 2 “kit”s. Each "kit" has 2 “operator”s. Each operator performs the experiment in 2 “lab”s. And,each of these combinations uses a different “part” of the sample-that’s a total of 8 different “part”s from the sample run randomly on different “kit-operator-lab” combinations. (another option is to use 4 "part"s from the same sample and run the same combination on the 2 kits)

All the variables "kit, operator, lab" are random effects and I would like to estimate the variance due to each of these.

It would be very helpful to know how to specify these effects in the "random" statement. Something like this?

proc mixed;

class kit operator room part;

model y=;

random kit operator*lab part(kit);

run;

Thanks!

Contributor
Posts: 40

Setting up the random effects model in Proc mixed

I wanted to correct the class statement. It should read:

class kit operator lab part;

Please share your thoughts on what is the best way to setup the model.

Thanks!

Frequent Contributor
Posts: 75

Setting up the random effects model in Proc mixed

would you be able to show us your data (or at least a subset of it)?

Contributor
Posts: 40

Setting up the random effects model in Proc mixed

I am afraid I won't be able to share the data. I realize that I will have to run diagnostics to make sure the model is correct and also modify the mdoel in case there are convergence issues. But, I just want to setup a model that is correct for the design. Thanks

Frequent Contributor
Posts: 75

Setting up the random effects model in Proc mixed

no problem. i am just having a hard time visualizing your data structure. hopefully some of the other users on here can chime in. can you show some made up dummy data (something that shows your data structure)?

Ask a Question
Discussion stats
  • 4 replies
  • 158 views
  • 0 likes
  • 2 in conversation