08-08-2016 11:42 AM
Getting a lot of observations that are dropped when using PROC GLIMMIX procedure. The note reads:
08-08-2016 01:53 PM
It would be helpful to know what options you are using for PROC GLIMMIX and the other procedure.
For example, the "not a proportion" note might indicate that you are using the DIST=BINOMINAL option instead of the DIST=BINARY option. There was an post last year in which the OP had a similar note.
08-08-2016 01:57 PM
Doing a quasi-likelihood regression, so here are the options:
PROC GLIMMIX DATA = TEST;
MODEL QALY =
/LINK = LOGIT S DIST = BETA;
output out=fracout pred(ilink)=pred lcl(ilink)=lower ucl(ilink)=upper;
08-08-2016 02:26 PM
With link=logit and dist=beta, I would guess that the values falling into the "not a proportion" are greater than 1 or less than 0. That would lead to this. As far as missing response/missing fixed effects part, it looks like that a large part of these overlap. OLS would not detect that a large part of the values did not fit the specified distribution.
08-08-2016 03:43 PM
You seem to have several observations where your response is larger than 1 or smaller than 0. That is why these are "not proportions". Probably you made an error in calculating the proprotion. And more seriously: the beta conditional distribution is defined between 0 and 1 (0 < prop < 1). That means that all 0s and 1s are converted to missing values. You cannot use the beta distribution if you have 0 and 1 unless you want to throw away data. Many references fail to make this clear. To get around this, you only have ad hoc solutions. For instance: Convert all 0s to a very small number (smaller than the smallest nonzero real value that you could observe). Same idea for the 1s. Of course, this is creating artificial data. Would be OK for occassional 0s and 1s, but it appears that you have many.