# NOT RUN {
## load verbal aggression data
data("VerbalAggression", package = "psychotools")
## fit a rasch model to dichotomized verbal aggression data
rmod <- raschmodel(VerbalAggression$resp2)
## extract the discrimination parameters
dp <- discrpar(rmod)
## extract the variance-covariance matrix
head(vcov(dp))
# }
Run the code above in your browser using DataLab