In tmsalab/edm: Exploratory Diagnostic Models

title: "Assessing Spatial Data" author: "JJB" date: "r Sys.Date()" output: rmarkdown::html_vignette vignette: > %\VignetteIndexEntry{Assessing Spatial Data} %\VignetteEngine{knitr::rmarkdown} %\VignetteEncoding{UTF-8}

Install the required packages

# install.packages("devtools")

devtools::install_github("tmsalab/edm", auth_token = "secret")
devtools::install_github("tmsalab/edmdata", auth_token = "secret")

Using the `edm` package

Load the packages

# Load required packages
library("edm")
library("edmdata")

# Load data
data("items_revised_psvtr", package="edmdata")

Perform Estimation Routine

# Version of package used
message("Running simulation with edm v",
        as.character(utils::packageVersion("edm")),
        sep = "")


### Configuration

# Burn in
burnin = 20000

# Length of chain
chain_length = 10000

# Number of attributes
K = 4

# Coerce to matrix
trial_matrix = as.matrix(trial_matrix)

# Launch job
edina_model = edina(trial_matrix, K, burnin, chain_length)

Evaluate Results of Q and IRT Estimation

#### Parameters

## Guessing and Slipping
out$coefficients

# Guessing
out$coefficients[, 1]

# Slipping
out$coefficients[, 2]

# Pi
out$pi_classes

# Q Matrix
out$est_Q

### Print methods
print(edina_model)

### Calculate model heuristic with default alpha = 0.05
summary(edina_model)

### Calculate model heuristic
summary(edina_model, alpha = 0.10)