Could anyone help me on how to include offset parameter in NLMIXED model with negative binomial? Is the below inclusion of log_TotalWordsInSentence (log of number of words in a sentence) variable correct?
proc nlmixed data=myData;
xb = b0 + b1a*(Gender = 'female') + u;
mu = exp(xb + log_TotalWordsInSentence);
m = 1/alpha;
ll = lgamma(CountPronounsInSentence+m)-lgamma(CountPronounsInSentence+1)-lgamma(m)
+ CountPronounsInSentence*log(alpha*mu)-(CountPronounsInSentence+m)*log(1+alpha*mu);
model CountPronounsInSentence ~ general(ll);
random u ~ normal(0,s2u) subject=person_id;
run;
sorry, i dont visit the message boards as frequently as i'd like to.... The code looks ok to me. Personally for clarity i might write it as:
bounds b2<=1, b2>=1;
xb = b0 + b1a*(Gender = 'female') + b2*log_TotalWordsInSentence + u;
mu = exp(xb);
and a minor thing: I'd reduce the length of the variable name 'CountPronounsInSentence' to make it readable, i didn't bother to read the loglikelihood equation ("ll = ")... because it's hard to read
@pau13rown is probably a better person to answer this, but your formula looks okay to me. Alternatively, instead of
mu = exp(xb + log_TotalWordsInSentence);
you could use
mu = TotalWordsInSentence*(exp(xb));
which is equivalent (take log of both sides, bring log_TotalWordsInSentence to the left, combine logs, then exponentiate.)
sorry, i dont visit the message boards as frequently as i'd like to.... The code looks ok to me. Personally for clarity i might write it as:
bounds b2<=1, b2>=1;
xb = b0 + b1a*(Gender = 'female') + b2*log_TotalWordsInSentence + u;
mu = exp(xb);
and a minor thing: I'd reduce the length of the variable name 'CountPronounsInSentence' to make it readable, i didn't bother to read the loglikelihood equation ("ll = ")... because it's hard to read
Are you ready for the spotlight? We're accepting content ideas for SAS Innovate 2025 to be held May 6-9 in Orlando, FL. The call is open until September 25. Read more here about why you should contribute and what is in it for you!
ANOVA, or Analysis Of Variance, is used to compare the averages or means of two or more populations to better understand how they differ. Watch this tutorial for more.
Find more tutorials on the SAS Users YouTube channel.