Bioequivalence and Bioavailability Forum • R code quest

R code quest [Regulatives / Guidelines]

posted by d_labes – Berlin, Germany, 2010-02-25 17:34 (5969 d 14:18 ago) – Posting: # 4819
Views: 6,772

Dear ElMaestro,

❝ Lund <- function (Linmo, alpha, Verbose)

❝ {

❝ ...

❝ q=qr(model.matrix(Linmo))$rank

❝ n=length(residuals(Linmo))

❝ ...

why did you use this here? :ponder:

:ponder:

Is there any hidden secret behind?

I would do:

   summary(Linmo)$df[2]  # error degrees of freedom



   ...

   t <- qt((1-alpha/n/2),df-1)    

   isr2 <- df*(t^2)/(t^2+df-1)  

   ...

—
Regards,

Detlew

Complete thread:

Canada: Draft Guidances published Helmut 2010-01-28 00:52
- Outlier earlybird 2010-02-09 08:35
  - Outlier (tests?) Helmut 2010-02-09 14:55
    - Outlier (tests?) earlybird 2010-02-09 16:45
      - Outlier (tests?) Helmut 2010-02-09 17:05
        
        Outlier (tests?) earlybird 2010-02-09 17:18
        Outlier "im Quadrat" d_labes 2010-02-10 13:51
        
        Hotelling's T² Helmut 2010-02-11 17:56
    - Lund! d_labes 2010-02-10 13:36
      - Lund! Helmut 2010-02-12 20:04
        
        intercept ElMaestro 2010-02-12 20:29
        
        intercept Helmut 2010-02-13 00:03
        
        intercept ElMaestro 2010-02-13 01:04
        
        Lund's critical values, technically d_labes 2010-02-22 12:52
        
        R code ElMaestro 2010-02-25 14:22
        
        R code questd_labes 2010-02-25 16:34
        
        R code quest ElMaestro 2010-02-25 17:06

Admin contact
23,656 posts in 4,994 threads, 1,571 registered users;
200 visitors (0 registered, 200 guests [including 16 identified bots]).
Forum time: 08:52 CEST (Europe/Vienna)

Try to learn something about everything
and everything about something. Thomas Henry Huxley

The Bioequivalence and Bioavailability Forum is hosted by
BEBAC

Ing. Helmut Schütz