Bioequivalence and Bioavailability Forum • Period DFs

Period DFs - EMA code literally in R [Regulatives / Guidelines]

posted by ElMaestro – Denmark, 2011-08-24 15:04 (5012 d 15:24 ago) – Posting: # 7303
Views: 7,190

Dear d_labes

❝ factor(b$sequence) vanished in nirvana :confused:

❝ Clever or bug?

I have no idea. I am one big question mark.

❝ muddle2 <- lm(log(Cmax)~ sequence + subject%in%sequence + period, data=refdata)

My eyes, my eyes, it hurts it hurts! Subject%in%sequence when all subjects have unique identifiers... Initial muddle matrix will contain 77 columns of pure zeros. Makes no practical difference for the results, though.

❝ drop1(muddle2)

drop1 does not work like in SAS when it comes to sequence. When R drops Sequence to get the deltaSS for this factor it does not automatically recognise that Sequence is already in the model because it is contained in Subject. SAS does that (or at least: should do it, but might obviously be failing here).

Your example with first the sequential muddle with period last and then a drop1 is truly brilliant. Here I would expect exactly the same SS for period between the two anovas. It just isn't the case.

Where's hrotter and jdetlor? There must be an explanation.

—
Pass or fail!
ElMaestro

Complete thread:

EMA replicate surprise d_labes 2011-08-23 13:01 [Regulatives / Guidelines]
- 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 DFs ElMaestro 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 RElMaestro 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