-
Notifications
You must be signed in to change notification settings - Fork 0
/
dexter.qmd
42 lines (36 loc) · 1.65 KB
/
dexter.qmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
## Analysing with dexter
Dexter is an R (R Development Core Team 2005) package intended as a robust and fairly comprehensive system for managing and analyzing test data organized in booklets. It includes facilities for importing and managing test data, assessing and improving the quality of data through basic test-and-item analysis, fitting an IRT model, and computing various estimates of ability.
Work through the examples here: <https://cran.r-project.org/web/packages/dexter/vignettes/dexter.html> and in this paper: <https://dexter-psychometrics.github.io/dexter/articles/blog/2023-05-20-dexter-paper>
## Extension exercise
When you have finished, work through the data on sentence structure and produce a report in Quarto.
```{r}
library(dplyr)
library(dexter)
head(verbAggrRules, 10)
db = start_new_project(verbAggrRules, "verbAggression.db", person_properties=list(gender="unknown"))
add_booklet(db, verbAggrData, "agg")
add_item_properties(db, verbAggrProperties)
head(verbAggrProperties)
get_booklets(db)
head(get_items(db))
get_persons(db) %>%
glimpse()
tt = tia_tables(db)
tt$booklets
tt$items
distractor_plot(db, 'S1DoShout')
m = fit_inter(db, booklet_id=='agg')
plot(m, "S1DoScold", show.observed=TRUE)
plot(m, 'S1DoCurse', summate=FALSE)
mSit = fit_domains(db, item_property= "situation")
plot(mSit)
parms = fit_enorm(db)
parms_gibbs = fit_enorm(db, method='Bayes')
head(coef(parms_gibbs))
pv = plausible_values(db, parms)
plot(density(pv$PV1), bty='l', main='verbal aggression', xlab='plausible value')
pv = merge(pv, get_persons(db))
boxplot(PV1~gender, data=pv, border='black')
profile_plot(db, item_property='mode', covariate='gender')
close_project(db)
```