eba (version 1.5-2)

thurstone: Thurstone-Mosteller model (Case V)

Description

Fits a Thurstone-Mosteller model (Case V) by maximum likelihood.

Usage

thurstone(M)

Arguments

M

a square matrix or a data frame consisting of absolute choice frequencies; row stimuli are chosen over column stimuli

Value

• estimate
  a vector of parameter estimates (scale values), first element is set to zero
• goodness.of.fit
  the goodness of fit statistic including the likelihood ratio fitted vs. saturated model (-2logL), the degrees of freedom, and the p-value of the corresponding chi-square distribution
• tm.glm
  the output from a call to glm

Details

thurstone provides a front end for glm.

See Critchlow & Fligner (1991) for more details.

References

Critchlow, D.E., & Fligner, M.A. (1991). Paired comparison, triple comparison, and ranking experiments as generalized linear models, and their implementation in GLIM. Psychometrika, 56, 517--533.

See Also

eba, strans, pcX, kendall.u, circular, glm.

Examples

## Taste data (David, 1988, p. 116)
dat <- matrix(c( 0,  3,  2,  2,
                12,  0, 11,  3,
                13,  4,  0,  5,
                13, 12, 10,  0), 4, 4, byrow=TRUE)
thurstone(dat)  # Thurstone-Mosteller model fits OK