I am thinking that my confusion come from the 2 following responses : y= (#event / Person-year) and y=#event. In your last post, you propose the model proc glimmix data=count_data; class study; model LOGtox = / offset=logPA2 solution; random intercept / subject=study ; * to test heterogeneity between study; weight inv_varln; * weight = inverse variance; parms (1) (1) / hold=2; run; But I have an approximation of the variance of tox/person-year or log(tox/person-year). So I have to peform a model where the response is log(tox/person-year) and delete the offset. proc glimmix data=count_data; class study; model LOG(tox/person-years) = / solution; random intercept / subject=study ; * to test heterogeneity between study; weight inv_varln; * weight = inverse variance; parms (1) (1) / hold=2; run; In my last posts I was wondering whether I have to consider the denominator of the response (offset) in the function variance determined by the software or whether the software takes into account only the numerator of incidence rate i.e #events.
... View more