Displaying 2 results from an estimated 2 matches for "lm21".
Did you mean:
lm01
2007 Oct 29
1
lm design matrix bug?
...y[idx] = 30*i+10*j+100*i*j+30*k-60
x1[idx] = i
x2[idx] = j
idx = idx+1
}
}
}
lm11 = lm(y ~ factor(x1)*factor(x2), x=1)
summary(lm11)
unique(predict(lm11))
X = lm11$x; X
P = solve(t(X)%*%X) %*% t(X); round(P,3)
y[3] = NA
y[6] = NA
y[12] = NA
y[18] = NA
y[24] = NA
y[27] = NA
lm21 = lm(y ~ factor(x1)*factor(x2), x=1)
summary(lm21)
unique(predict(lm21))
X = lm21$x; X
P = solve(t(X)%*%X) %*% t(X); round(P,3)
2002 Sep 14
0
p.s. regarding stripchart missing-data report (PR#2019)
...e+minrisk.numbers.sm+QOL.numbers.sm
#+worry.numbers.sm+trust.numbers.sm+mss.m.score)
#anova(lm18)
#lm19 <- lm(toohigh.sm ~ stai.score+educ+sex+truste.score+minrisk.numbers.sm+QOL.numbers.sm
#+worry.numbers.sm+comfort.numbers.sm)
#anova(lm19)
lm20 <- lm(toohigh.sm ~ stai.score)
anova(lm20)
lm21 <- lm(toohigh.sm ~ educ)
anova(lm21)
t.test(toohigh.sm[sex==1], toohigh.sm[sex==2])
lm22 <- lm(toohigh.sm ~ educ)
anova(lm22)
lm23 <- lm(toohigh.sm ~ lot.score)
anova(lm23)
lm24 <- lm(toohigh.sm ~ denial.score)
anova(lm24)
lm25 <- lm(toohigh.sm ~ truste.score)
anova(lm25)
lm26 &...