drugrisk {eba} | R Documentation |

## Perceived Health Risk of Drugs

### Description

In summer 2007, a survey was conducted at the Department of Psychology,
University of Tuebingen. Hundred and ninety-two participants were presented
with all 15 unordered pairs of the names of six drugs or substances and
asked to choose the drug they judged as more dangerous for their health.
The six drugs were alcohol (alc), tobacco (tob), cannabis (can), ecstasy
(ecs), heroine (her), and cocaine (coc). Choice frequencies were aggregated
in four groups defined by gender and age.

### Usage

`data(drugrisk)`

### Format

A 3d array consisting of four square matrices of choice
frequencies (row drugs are judged over column drugs):

`drugrisk[, , group = "female30"]`

holds the choices
of the 48 female participants up to 30 years of age.

`drugrisk[, , group = "female31"]`

holds the choices
of the 48 female participants from 31 years of age.

`drugrisk[, , group = "male30"]`

holds the choices
of the 48 male participants up to 30 years of age.

`drugrisk[, , group = "male31"]`

holds the choices
of the 48 male participants from 31 years of age.

### Source

Wickelmaier, F. (2008).
Analyzing paired-comparison data in R using probabilistic choice
models. Presented at the *R User Conference 2008*,
August 12-14, Dortmund, Germany.

### Examples

```
data(drugrisk)
## Bradley-Terry-Luce (BTL) model
btl <- eba(drugrisk[, , group = "male30"])
## Elimination-by-aspects (EBA) model, 1 additional aspect
A1 <- list(c(1), c(2,7), c(3,7), c(4,7), c(5,7), c (6,7))
eba1 <- eba(drugrisk[, , group = "male30"], A1)
## EBA model, 2 additional aspects
A2 <- list(c(1), c(2,7), c(3,7), c(4,7,8), c(5,7,8), c(6,7,8))
eba2 <- eba(drugrisk[, , group = "male30"], A2)
## Model selection
anova(btl, eba1, eba2)
## Utility scale values (BTL for females, EBA for males)
dotchart(cbind(
apply(drugrisk[, , 1:2], 3, function(x) uscale(eba(x), norm = 1)),
apply(drugrisk[, , 3:4], 3, function(x) uscale(eba(x, A2), norm = 1))
), xlab="Utility scale value (BTL and EBA models)",
main="Perceived health risk of drugs",
panel.first = abline(v = 1, col = "gray"), log = "x")
mtext("(Wickelmaier, 2008)", line = 0.5)
```

[Package

*eba* version 1.10-0

Index]