R Dataset / Package lme4 / VerbAgg


On this Picostat.com statistics page, you will find information about the VerbAgg data set which pertains to Verbal Aggression item responses. The VerbAgg data set is found in the lme4 R package. You can load the VerbAgg data set in R by issuing the following command at the console data("VerbAgg"). This will load the data into a variable called VerbAgg. If R says the VerbAgg data set is not found, you can try installing the package by issuing this command install.packages("lme4") and then attempt to reload the data. If you need to download R, you can go to the R project website. You can download a CSV (comma separated values) version of the VerbAgg R data set. The size of this file is about 437,902 bytes.

Verbal Aggression item responses


These are the item responses to a questionaire on verbal aggression. These data are used throughout De Boeck and Wilson, Explanatory Item Response Models (Springer, 2004) to illustrate various forms of item response models.


A data frame with 7584 observations on the following 13 variables.


the subject's Trait Anger score as measured on the State-Trait Anger Expression Inventory (STAXI)


the subject's gender - a factor with levels M and F


the item on the questionaire, as a factor


the subject's response to the item - an ordered factor with levels no < perhaps < yes


the subject identifier, as a factor


behavior type - a factor with levels curse, scold and shout


situation type - a factor with levels other and self indicating other-to-blame and self-to-blame


behavior mode - a factor with levels want and do


dichotomous version of the response - a factor with levels N and Y




De Boeck and Wilson (2004), Explanatory Item Response Models, Springer.


## Show how  r2 := h(resp) is defined:
with(VerbAgg, stopifnot( identical(r2, {
     r <- factor(resp, ordered=FALSE); levels(r) <- c("N","Y","Y"); r})))xtabs(~ item + resp, VerbAgg)
xtabs(~ btype + resp, VerbAgg)
round(100 * ftable(prop.table(xtabs(~ situ + mode + resp, VerbAgg), 1:2), 1))
person <- unique(subset(VerbAgg, select = c(id, Gender, Anger)))
densityplot(~ Anger, person, groups = Gender, auto.key = list(columns = 2),
            xlab = "Trait Anger score (STAXI)")if(lme4:::testLevel() >= 3) { ## takes about 15 sec
print(fmVA <- glmer(r2 ~ (Anger + Gender + btype + situ)^2 +
 		   (1|id) + (1|item), family = binomial, data =
		   VerbAgg), corr=FALSE)
                       ## much faster but less accurate
print(fmVA0 <- glmer(r2 ~ (Anger + Gender + btype + situ)^2 +
                    (1|id) + (1|item), family = binomial, data =
                    VerbAgg, nAGQ=0L), corr=FALSE)

Dataset imported from https://www.r-project.org.

Title Authored on Content type
R Dataset / Package psych / bfi March 9, 2018 - 1:06 PM Dataset
OpenIntro Statistics Dataset - scotus_healthcare August 9, 2020 - 2:38 PM Dataset
R Dataset / Package psych / withinBetween March 9, 2018 - 1:06 PM Dataset
R Dataset / Package Stat2Data / Kids198 March 9, 2018 - 1:06 PM Dataset
R Dataset / Package Ecdat / Wages1 March 9, 2018 - 1:06 PM Dataset
Attachment Size
dataset-97081.csv 427.64 KB
Dataset License
GNU General Public License v2.0
Documentation License
GNU General Public License v2.0