Displaying 1 result from an estimated 1 matches for "aptest".
Did you mean:
aaptest
2013 Nov 08
1
Different output from lm() and lmPerm lmp() if categorical variables are included in the analysis
...in the case of including categorical variables:
require(lmPerm)
set.seed(42)
testx1 <- rnorm(100,10,5)
testx2 <- c(rep("a",50),rep("b",50))
testy <- 5*testx1 + 3 + runif(100,-20,20)
test <- data.frame(x1=testx1,x2=
testx2,y=testy)
atest <- lm(y ~ x1*x2,data=test)
aptest <- lmp(y ~ x1*x2,data=test,perm = "", seqs = TRUE, center = FALSE)
summary(atest)
Call:
lm(formula = y ~ x1 * x2, data = test)
Residuals:
Min 1Q Median 3Q Max
-17.1777 -9.5306 -0.9733 7.6840 22.2728
Coefficients:
Estimate Std. Error t value Pr(>...