Hi there! I am looking to develop a growth curve model for aggressive behavior over time using a 2 level multilevel model. Each individual was assessed up to 6 time points, so Level 1 = Wave and Level 2 = Participant. The nature of the outcome (aggression) is a count of behaviors that is heavily skewed, so I've been told I should use a negative binomial model. I would like to first evaluate whether there is significant variability in slopes (i.e., test for random slopes) in a null model, and then ultimately include other variables to potentially explain/account for that variability in slopes in subsequent steps. I think I know how to include random intercepts in the null model using PROC GLIMMIX with the following code: proc glimmix data=AggData method=quad noclprint; title 'Time Null model'; class PARTID; model Aggression = Wave / dist=nb link=log s; random intercept / TYPE=UN subject= PARTID; run; When I run this, I get all of the fit statistics, a table of Covariance Parameter Estimates (which I * think * relate to the random intercept?), and the solutions for Fixed Effects. My questions are: 1) How do I incorporate random slope into this model? I want to allow for each participant to have their own slope over time, but I'm not sure of how to code for that. 2) Is there a way to obtain significance tests or a way to evaluate if the random intercept and slope are significant (i.e., that there is in fact sufficient variability in the participant intercepts and slopes) in the output? I'm relatively new to SAS so have been struggling to make sense of the online resources and programming documentation to find exactly what I'm looking for. Hoping someone here may be able to help me. Thank you!
... View more