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
[[alternative HTML version deleted]]
______________________________________________
[email protected] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
