fit_TIRT_stan: Fit Thurstonian IRT models in Stan

View source: R/stan.R

fit_TIRT_stanR Documentation

Fit Thurstonian IRT models in Stan

Description

Fit Thurstonian IRT models in Stan

Usage

fit_TIRT_stan(data, init = 0, ...)

Arguments

data

An object of class 'TIRTdata'. see make_TIRT_data for documentation on how to create one.

init

Initial values of the parameters. Defaults to 0 as it proved to be most stable.

...

Further arguments passed to rstan::sampling.

Value

A 'TIRTfit' object.

Examples

# load the data
data("triplets")

# define the blocks of items
blocks <-
  set_block(c("i1", "i2", "i3"), traits = c("t1", "t2", "t3"),
          signs = c(1, 1, 1)) +
  set_block(c("i4", "i5", "i6"), traits = c("t1", "t2", "t3"),
            signs = c(-1, 1, 1)) +
  set_block(c("i7", "i8", "i9"), traits = c("t1", "t2", "t3"),
            signs = c(1, 1, -1)) +
  set_block(c("i10", "i11", "i12"), traits = c("t1", "t2", "t3"),
            signs = c(1, -1, 1))

# generate the data to be understood by 'thurstonianIRT'
triplets_long <- make_TIRT_data(
  data = triplets, blocks = blocks, direction = "larger",
  format = "pairwise", family = "bernoulli", range = c(0, 1)
)


# fit the data using Stan
fit <- fit_TIRT_stan(triplets_long, chains = 1)
print(fit)
predict(fit)



thurstonianIRT documentation built on Aug. 22, 2023, 5:08 p.m.