Bioequivalence and Bioavailability Forum • R code for lm.mod does not make sense

R code for lm.mod does not make sense [🇷 for BE/BA]

posted by dmfisher – USA, 2015-06-09 05:55 (3664 d 04:33 ago) – Posting: # 14930
Views: 52,598
(edited on 2015-06-09 07:46)

The R code for lm.mod is (obtained by typing lm.mod at a prompt):

    cat("\n")

    target.lm <- lm(target ~ seq + subj:seq + prd + drug, data = inputdata)

    BearAnova = anova(target.lm)

    row.names(BearAnova)[4] = "subj(seq)"

    show(BearAnova[2:5, ])

    cat("\n")

    colnames(BearAnova) <- c(" DF", "   Type I SS", " Mean Square", 

        " F Value", " Pr > F")

    show(BearAnova[2:4, ])

    cat("\n")

    colnames(BearAnova) <- c(" DF", " Type III SS", " Mean Square", 

        " F Value", " Pr > F")

    show(BearAnova[2:4, ])

    cat("-----------------------------------------------------------\n")

    cat("\n")

    cat("Tests of Hypothesis using the Type I MS for \n")

    cat("SUBJECT(SEQUENCE) as an error term\n\n")

    show(summary(aov(target ~ seq, data = inputdata)))

    cat("\n")

    cat("Tests of Hypothesis using the Type III MS for\n")

    cat("SUBJECT(SEQUENCE) as an error term\n\n")

    show(summary(aov(target ~ seq, data = inputdata)))

    cat("---\n")

    cat("Sum Sq. = Type III SS\n\n")

    cat(" Tests for carryover effects and direct formulation effects\n")

    cat("------------------------------------------------------------")

    show(summary(aov(target ~ seq * drug + Error(subj), data = inputdata)))

    cat("------------------------------------------------------------\n")

There are several instances of
show(BearAnova ...)
The column names are changed between these commands but the contents of the object are unchanged. This makes no sense — how could the values be the same for Type I and Type III SS?

Edit: category changed [Ohlbe]
Edit: Please follow the Forum’s Policy. [Helmut]

Complete thread:

R code for lm.mod does not make sensedmfisher 2015-06-09 03:55
- R code for lm.mod does not make sense yjlee168 2015-06-09 17:57
  - R code for lm.mod does not make sense ElMaestro 2015-06-09 21:07
    - R code for lm.mod does not make sense yjlee168 2015-06-10 11:34
      - R code for lm.mod does not make sense ElMaestro 2015-06-10 14:22
        
        R code for lm.mod does not make sense yjlee168 2015-06-11 09:44
        
        Reference datasets Helmut 2015-06-11 11:50
        
        Reference datasets ElMaestro 2015-06-11 12:01
        
        Reference datasets yjlee168 2015-06-11 22:16
        Type III kuddel muddel d_labes 2015-06-12 09:04
        
        Type III kuddel muddel Helmut 2015-06-12 11:27
        
        ANOVA F-test results and regulators d_labes 2015-06-12 15:23
        
        ANOVA F-test results and regulators ElMaestro 2015-06-12 15:29
        
        ANOVA F-test results and regulators d_labes 2015-06-12 15:44
        
        ANOVA F-test results and regulators ElMaestro 2015-06-12 15:52
        
        Type III again ElMaestro 2015-06-12 12:18
        
        SAS OT d_labes 2015-06-12 15:32
        
        Reference datasets yjlee168 2015-06-11 19:56
        Test hypothesis of using Type I/III MS... yjlee168 2015-06-12 09:01
  - R code for lm.mod does not make sense dmfisher 2015-06-10 01:38