Bioequivalence and Bioavailability Forum • Period DFs

Period DFs [Regulatives / Guidelines]

posted by ElMaestro – Denmark, 2011-08-23 18:58 (5425 d 22:27 ago) – Posting: # 7299
Views: 8,775

This post is extensively edited as I did initially not catch that we are only dealing with observations of the Reference formulation.

Dear d_labes and HS,

Took the file provided by HS and and rn it through R as follows:




a=read.table("ema1.csv", header=T)

b=a[grep("R", a$Formulation, ignore.case=T),]


Muddle=lm(log(b$Data)~factor(b$Period)+factor(b$Sequence)+factor(b$Subject))

anova(Muddle)

Analysis of Variance Table



Response: log(b$Data)

                  Df  Sum Sq Mean Sq F value    Pr(>F)    

factor(b$Period)   3   1.570 0.52326  2.6253   0.05703 .  

factor(b$Subject) 75 118.662 1.58216  7.9380 2.454e-16 ***

Residuals         71  14.151 0.19931                      

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

drop1(Muddle, test="F")

Single term deletions



Model:

log(b$Data) ~ factor(b$Period) + factor(b$Sequence) + factor(b$Subject)

                   Df Sum of Sq     RSS      AIC F value     Pr(F)    

<none>                           14.151 -196.125                      

factor(b$Period)    2     0.773  14.924 -192.150  1.9383    0.1515    

factor(b$Sequence)  0     0.000  14.151 -196.125                      

factor(b$Subject)  75   118.662 132.813  -10.254  7.9380 2.454e-16 ***

---

Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

What it really, really strange is this:

M=model.matrix(Muddle)

qr(M)$rank

[1] 79

Nützt nüchts. Ich werde's nie verstehen. Ist auch scheissegal.

—
Pass or fail!
ElMaestro

Complete thread:

EMA replicate surprise d_labes 2011-08-23 13:01
- EMA replicate surprise ElMaestro 2011-08-23 13:15
  - Big bäng CV d_labes 2011-08-23 13:58
    - Big bäng CV ElMaestro 2011-08-23 15:02
      - Period DFs d_labes 2011-08-23 15:28
        
        Period DFs Helmut 2011-08-23 15:36
        Period DFs ElMaestro 2011-08-23 16:26
        Period DFsElMaestro 2011-08-23 16:58
        
        Period DFs Helmut 2011-08-24 00:44
        
        Lessons learned? d_labes 2011-08-24 11:24
        
        Lessons learned? Helmut 2011-08-24 13:51
        Lessons learned? ElMaestro 2011-08-28 21:34
        
        Period DFs - EMA code literally in R d_labes 2011-08-24 10:47
        
        Period DFs - EMA code literally in R ElMaestro 2011-08-24 13:04
        
        call for expert’s opinions Helmut 2011-08-24 14:01
        
        call for expert’s opinions jdetlor 2011-11-01 03:05
        
        The finance sector ElMaestro 2011-11-02 09:48