Mohamed Yehia
★

Egypt,
2018-01-23 10:48
(2336 d 19:16 ago)

Posting: # 18261
Views: 4,917

## Which is better for missing periods PROC GLM or PROC MIXED [RSABE / ABEL]

Hi All,

Hope everything is fine and going well.

For Partial and full replicate studies, if we have some subjects with missing periods, i mean for example subjects who attended only 2 periods in partial replicate

I have only one question:

1) Which method to use following the EMA guidelines; method A "GLM" or method B "Mixed"?.

I am not an expert in statistics, however I read some papers stated "PROC GLM uses least squares or method of moments to fit general linear models. On the other hand, PROC MIXED uses Restricted (or residual) Maximum Likelihood (REML). PROC MIXED is recommended to avoid pitfalls of PROC GLM."

Thanks.
Helmut
★★★

Vienna, Austria,
2018-01-23 12:25
(2336 d 17:40 ago)

@ Mohamed Yehia
Posting: # 18263
Views: 4,171

## EMA: Method A or B

Hi Mohamed,

❝ For Partial and full replicate studies, if we have some subjects with missing periods, i mean for example subjects who attended only 2 periods in partial replicate

First answering a question you didn’t ask. In partial replicate designs (TRR|RTR|RRT) I suggest to specify two analysis sets if data of the third period in sequence RRT is missing:
1. Exclude subjects from the assessment of BE since the GL requires at least one treatment of T and R.
2. Keep subjects for the estimation of CVwR.

❝ 1) Which method to use following the EMA guidelines; method A "GLM" or method B "Mixed"?.

❝ I want a clarified answer please with guideline reference if applicable.

Good question. The Q&A document states:

The analysis […] show that this approach (Method A) is feasible even for unbalanced replicate design studies. The advantage of this approach is that it is straightforward and that it appears to be software and software option independent. A simple linear mixed model, which assumes identical within-subject variability (Method B), may be acceptable as long as results obtained with the two methods do not lead to different regulatory decisions. However, in borderline cases and when there are many included subjects who only provide data for a subset of the treatment periods, additional analysis using Method A might be required.

• Not only Method B assumes identical within-subject variability – the same is true for Method A. A major advantage of a full replicate design is that it allows to separately estimate s²wR and s²wT. IMHO, that’s a bizarre assumption which rarely holds true in practice (in most cases I have seen s²wR > s²wT). The estimation is possible either directly by Method C (which seemingly is not “compatible with [the] CHMP guideline”) or modifying the model given in Section 3.4. by including data of T only. The Q&A also states:
“An advantage of Method C is that it directly calculates s²wR. However, sometimes the algorithm fails to converge.”
This observation is correct only for partial replicate designs, since T is administered only once and the model is over-specified. If the algo converges (generally it does), the estimate of s²wR is correct and the estimate of s²wT is nonsense. However, Method C always converges for full replicate designs.
• Since one doesn’t know beforehand whether the outcome will be a “borderline case” (which “border”?) I guess one has to perform both Methods even if Method B was aimed at. The wording gives me (!) the impression that the PKWP assumes that Method A is more conservative than Method B. That’s not even true for the Q&A’s data set I: The CI obtained by Method B is wider (17.80%) than the one obtained by Method A (17.78%). We are working on a paper supporting validation of software for the EMA’s ABEL. In 57% of our reference data sets Method B lead to wider confidence intervals. However, in 96% of cases the BE decision (pass/fail) agreed. Interesting that one of our data sets would fail by Method B but pass by Method A.
From a purely statistical perspective my preferences are: Method C ≫ Method B > Method A. I don’t like the idea to treat subjects as a fixed effect. Given the observations from above, a sponsor probably would fair best with Method A only.

Dif-tor heh smusma 🖖🏼 Довге життя Україна!
Helmut Schütz

The quality of responses received is directly proportional to the quality of the question asked. 🚮
Science Quotes
Mohamed Yehia
★

Egypt,
2018-01-23 12:56
(2336 d 17:09 ago)

@ Helmut
Posting: # 18264
Views: 4,065

## EMA: Method A or B

❝ […] In partial replicate designs (TRR|RTR|RRT) I suggest to specify two analysis sets if data of the third period in sequence RRT is missing:

1. Exclude subjects from the assessment of BE since the GL requires at least one treatment of T and R.

2. Keep subjects for the estimation of CVwR.

Yes

❝ ❝ I want a clarified answer please with guideline reference if applicable.

❝ The Q&A document

❝ From a purely statistical perspective my preferences are: Method C ≫ Method B > Method A. I don’t like the idea to treat subjects as a fixed effect.

That's what I believe too

❝ Given the observations from above, a sponsor probably would fair best with Method A only.

Of course

Edit: Full quote removed. Please delete everything from the text of the original poster which is not necessary in understanding your answer; see also this post #5! doi corrected. [Helmut]