Anything wrong with the code?

Accepted Solution Solved
Reply
Occasional Contributor
Posts: 7
Accepted Solution

Anything wrong with the code?

Hi everyone! I have a quick question. I'm using this code to check for if there is a correlation between edad*musc, edad*sexo, musc*sexo and edad*musc*sexo. My question is, in the random statement, is it correct to write toro(edad) like I'm doing? My only random would be toro, which is my independent variable, but edad, as well as musc and sexo are all dependent variables.  

 

 

proc glimmix data=cruza;
class edad musc sexo toro;
model tern_co = edad | musc| sexo;
random toro(edad);
lsmeans edad | musc | sexo / pdiff lines;
proc glimmix data=cruza;
class edad musc sexo toro;
model Tern_cr = edad | musc| sexo;
random Toro(Edad);
lsmeans edad | musc| sexo / pdiff lines;
run;

 

 

Thanks. 


Accepted Solutions
Solution
‎03-20-2018 11:44 AM
Respected Advisor
Posts: 2,521

Re: Anything wrong with the code?

Posted in reply to ValeriaRodz

One concern is that you may be running low on degrees of freedom for error, and possibly eliminating the three-way interaction might help. Generally, I don't put three-way interactions into a model without having some subject matter reason for thinking it belongs. Other than that, it seems like the right model.

--
Paige Miller

View solution in original post


All Replies
Respected Advisor
Posts: 2,521

Re: Anything wrong with the code?

Posted in reply to ValeriaRodz

Impossible to answer

 

We don't know what the variables edad musc sexo and toro represent.

 

We don't know the design of your study either, which would also inform us about what the proper model and random statements are.

--
Paige Miller
Occasional Contributor
Posts: 7

Re: Anything wrong with the code?

Posted in reply to PaigeMiller
Hi! Sorry for not being specific. I'm studying the effect of castration
(sexo) and slaughter age (edad) in meat quality in two different muscles
(musc). We used 30 animals in total, 15 of them were castrated and 15
remained intact. So each animal is represented by the (toro) variable.
Super User
Posts: 22,579

Re: Anything wrong with the code?

Posted in reply to ValeriaRodz

That's a fairly complex model for only 30 observations...

Solution
‎03-20-2018 11:44 AM
Respected Advisor
Posts: 2,521

Re: Anything wrong with the code?

Posted in reply to ValeriaRodz

One concern is that you may be running low on degrees of freedom for error, and possibly eliminating the three-way interaction might help. Generally, I don't put three-way interactions into a model without having some subject matter reason for thinking it belongs. Other than that, it seems like the right model.

--
Paige Miller
☑ This topic is solved.

Need further help from the community? Please ask a new question.

Discussion stats
  • 4 replies
  • 85 views
  • 2 likes
  • 3 in conversation