dexter.R
In dexter: Data Management and Analysis of Tests

## ----setup, include=FALSE-----------------------------------------------------
library(knitr)
opts_chunk$set(echo = TRUE,fig.align='center', fig.width=6, fig.height=5)

if (requireNamespace("Cairo", quietly = TRUE)) 
{
   opts_chunk$set(dev='CairoPNG')
}

par_hook = function(before, options, envir)
{
  if(before)
  {
    do.call(par, options$par)
  }
}
knit_hooks$set(par = par_hook)


## ---- message=FALSE-----------------------------------------------------------
library(dplyr)
library(dexter)
head(verbAggrRules, 10)

## ---- eval=FALSE--------------------------------------------------------------
#  db = start_new_project(verbAggrRules, "verbAggression.db", person_properties=list(gender="unknown"))

## ---- include=FALSE-----------------------------------------------------------
db = start_new_project(verbAggrRules, ":memory:", 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)

## ---- eval=FALSE--------------------------------------------------------------
#  tt$booklets

## ---- echo=FALSE--------------------------------------------------------------
kable(tt$booklets, digits=3)

## ---- eval=FALSE--------------------------------------------------------------
#  tt$items

## ---- echo=FALSE--------------------------------------------------------------
kable(tt$items, digits=3)

## -----------------------------------------------------------------------------
distractor_plot(db, 'S1DoShout')

## -----------------------------------------------------------------------------
m = fit_inter(db, booklet_id=='agg')
plot(m, "S1DoScold", show.observed=TRUE)

## -----------------------------------------------------------------------------
plot(m, 'S1DoCurse', summate=FALSE)

## ----par=list(mfrow=c(2,2))---------------------------------------------------
mSit = fit_domains(db, item_property= "situation")
plot(mSit)

## ---- results='hide'----------------------------------------------------------
parms = fit_enorm(db)

## ---- results='hide'----------------------------------------------------------
parms_gibbs = fit_enorm(db, method='Bayes')

## ---- eval=F------------------------------------------------------------------
#  head(coef(parms_gibbs))

## ---- echo=F------------------------------------------------------------------
kable(head(coef(parms_gibbs)), digits=3)

## -----------------------------------------------------------------------------
pv = plausible_values(db, parms)
plot(density(pv$PV1), bty='l', main='verbal aggression', xlab='plausible value')

## ---- par=list(bty='n', fg='white')-------------------------------------------
pv = merge(pv, get_persons(db))

boxplot(PV1~gender, data=pv, border='black')

## ---- fig.width=5, fig.height=5-----------------------------------------------
profile_plot(db, item_property='mode', covariate='gender')

## -----------------------------------------------------------------------------
close_project(db)