Hello,
I am running error rate analysis. It is my results below. When I compare
aov1 and aov2, X square = 4.05, p = 0.044, which indicates that adding the
factor "Congruity" improved the fitting of model. However, the
following Z
value is less than 1 and p value for Z is 1, which means that
"Congruity" is
not significant at all. Therefore, these two parts are not consistent, one
is sig., the other is not. I cannot figure out why they are so different
and inconsistent? Which p value I should report? Any comments are
appreciated.
> options(digits=2)
> data <- read.delim("strooptaskdata2.txt", header=T)
> data <- data[data$condition != "both_incongru",]
> data <- data[data$condition != "syll_incongru",]
> data <- data[data$condition != "syll_congru",]
> data <- data[data$soa == "0",]
> attach(data)
> error <- numeric()
> error[RT==RT] <- 0
> error[RT < 0] <- 1
> x<-length (RT); x
[1] 2220> RT[RT<0]<-NA
> y<-sum(is.na(RT)); y
[1] 9> perc1<-y/x;perc1
[1] 0.0041>
> RT[RT<200]<-NA
> RT[RT>2000]<-NA
> z<-sum(is.na(RT));z
[1] 29> perc2<-z/x;perc2
[1] 0.013>
> congruity <- numeric()
> congruity[condition=="both_congru"] <- "congru"
> congruity[condition=="neutral"] <- "neutral"
> congruity <- as.factor(congruity)
> tapply(error*100, congruity, mean, na.rm=T)
congru neutral
0.00 0.51 > aov1 <- lmer(error ~ 1 + (1 | subj) + (1 | color),
family="binomial")
> aov2 <- lmer(error ~ congruity + (1 | subj) + (1 | color),
> family="binomial")
> anova(aov1, aov2)
Data:
Models:
aov1: error ~ 1 + (1 | subj) + (1 | color)
aov2: error ~ congruity + (1 | subj) + (1 | color)
Df AIC BIC logLik Chisq Chi Df Pr(>Chisq)
aov1 3 119.8 136.9 -56.9
aov2 4 117.8 140.6 -54.9 4.05 1 0.044 *
---
Signif. codes: 0 ?***? 0.001 ?**? 0.01 ?*? 0.05 ?.? 0.1 ? ? 1
> summary(aov2)
Generalized linear mixed model fit by the Laplace approximation
Formula: error ~ congruity + (1 | subj) + (1 | color)
AIC BIC logLik deviance
118 141 -54.9 110
Random effects:
Groups Name Variance Std.Dev.
subj (Intercept) 1.429 1.20
color (Intercept) 0.533 0.73
Number of obs: 2220, groups: subj, 37; color, 4
Fixed effects:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -22.0 2026.8 -0.01087 1
congruityneutral 15.8 2026.8 0.00781 1
Correlation of Fixed Effects:
(Intr)
congrtyntrl -1.000
--
View this message in context:
http://r.789695.n4.nabble.com/results-of-model-comparisons-and-summary-tp4663423.html
Sent from the R help mailing list archive at Nabble.com.