The package currently supports “case weights” - weights on individual observations. I wonder if there is a way to tweak the code so that I can put weights on clusters?
Specifically, the package currently can take a vector of observation weights w , and calculate the total log-likelihood as \sum log(w_ij * P(y_ij)) , where w_ij and y_ij are the weight and response for observation j from cluster i .
However, instead of weighting the probabilities on the observation level, I want to weight the log-likelihood on the cluster level, i.e., I want to calculate \sum w_i l_i ? Here w_i and l_i are the weight and the log-likelihood for cluster i .
Can anyone give a pointer on where I should get started? Thank you.