R code for lm.mod does not make sense [🇷 for BE/BA]
The R code for lm.mod is (obtained by typing lm.mod at a prompt):
There are several instances of
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]
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 [🇷 for BE/BA]
- 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
- ANOVA F-test results and regulators d_labes 2015-06-12 15:44
- ANOVA F-test results and regulators ElMaestro 2015-06-12 15:29
- ANOVA F-test results and regulators d_labes 2015-06-12 15:23
- Type III again ElMaestro 2015-06-12 12:18
- SAS OT d_labes 2015-06-12 15:32
- Type III kuddel muddel Helmut 2015-06-12 11:27
- Reference datasets yjlee168 2015-06-11 19:56
- Test hypothesis of using Type I/III MS... yjlee168 2015-06-12 09:01
- Reference datasets ElMaestro 2015-06-11 12:01
- Reference datasets Helmut 2015-06-11 11:50
- R code for lm.mod does not make sense yjlee168 2015-06-11 09:44
- 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-10 11:34
- R code for lm.mod does not make sense dmfisher 2015-06-10 01:38
- 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-09 17:57