Bioequivalence and Bioavailability Forum • R vs. Phoenix & SAS?

R vs. Phoenix & SAS? [🇷 for BE/BA]

posted by Helmut – Vienna, Austria, 2015-04-21 03:02 (3286 d 12:47 ago) – Posting: # 14725
Views: 22,815

Hi Yung-jin,

❝ ❝ For the EMA we would need an “all effects fixed” model (Method A of the Q&A).

❝ if EMA needs an "all effects fixed" model, why do we use a mixed model to do analysis?

Both the BE-GL and the Q&A state that all effect should be fixed. EMA generally accepts a mixed model, but may ask for recalculation in borderline cases (e.g., extreme imbalance, CI close to the AR). Sometimes assessors ask for recalculation even if the CI was very narrow and the PE was close to 1. Happened to me once. The Q&A specifically mentions SAS’ Proc GLM, not Proc Mixed.

❝ I don't remember lme() right now,…

In lme() it should be possible to specify all effects fixed. Actually the random effect(s) are optional. If is doesn’t work, lm() should do.

❝ Does EMA imply that we should use lm() to analyze replicate crossover dataset? :confused:

No. Both the FDA and the EMA do not mandate any software or routine within. Quote from the Q&A:

SAS (version 9.1, SAS Institute Inc, NC) was used in the previous computations. Results obtained by alternative, validated statistical programs are also acceptable except spreadsheets because outputs of spreadsheets are not suitable for secondary assessment.

❝ You should try to get funded from EMA or FDA for these projects. These are really great projects, indeed.

Gimme a break. I’m a lousy R-coder. Agencies want to get results from validated software. They don’t pay for its development.

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

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

Complete thread:

Bear vs. Phoenix & SAS Helmut 2015-04-20 17:34
- R vs. Phoenix & SAS? yjlee168 2015-04-20 19:36
  - R vs. Phoenix & SAS?Helmut 2015-04-21 01:02
    - lme() does not work with all fixed effects yjlee168 2015-04-21 23:41
      - lme() does not work with all fixed effects Astea 2016-11-04 00:13
        
        lmer: Method B (PE catched for imbalanced dataset!!!) and Method C mittyri 2016-11-05 17:38
        
        lmer: Method B (PE catched for imbalanced dataset!!!) and Method C Astea 2016-11-05 19:27
        
        lmer: Method B is ready for scaling mittyri 2016-11-05 20:01
        
        lmer: Method B is ready for scaling Astea 2016-11-06 11:50
        
        lmer: Method B is ready for scaling mittyri 2016-11-07 06:07
- info for lsmeans yjlee168 2015-04-20 21:34
  - info for lsmeans Helmut 2015-04-21 01:15
    - once more about R and replicate designes Astea 2016-11-02 23:43
      - once more about R and replicate designes VStus 2016-11-06 11:34
        
        Getting variance components from the lmer output StatR 2017-02-03 13:53
        
        Getting variance components from the lmer output VStus 2017-02-03 15:47
        
        Getting variance components from the lmer output StatR 2017-02-03 17:12
        
        Getting variance components d_labes 2017-02-07 11:16
        
        Getting variance components StatR 2017-02-07 11:36
        Getting variance components StatR 2017-02-08 08:41
        
        Getting variance components d_labes 2017-02-08 10:13
        
        Getting variance components StatR 2017-02-08 10:19
        
        Data structure Helmut 2017-02-08 10:33
        
        Data structure StatR 2017-02-08 10:49