Original Values

  • The original values from the combine data set
1
max.value 56.55
PHDD.mean 0.18
DPDD.mean 4.95
no.drinks.mean 190.95
no.days.drinking.mean 28.31
Cor.DrInC.PHD 0.42
Abstinence.W1 0.24
Abstinence.W1.1 0.24
Low.Risk.W1 0.11
Medium.Risk.W1 0.15
High.Risk.W1 0.18
Very.High.Risk.W1 0.31
Abstinence.W2 0.24
Low.Risk.W2 0.05
Medium.Risk.W2 0.09
High.Risk.W2 0.14
Very.High.Risk.W2 0.48
Abstinence.W3 0.24
Low.Risk.W3 0.45
Medium.Risk.W3 0.11
High.Risk.W3 0.10
Very.High.Risk.W3 0.09
tt.N_PHD 0.15
tt.Tx_PHD 0.23
tt.N_DPD 0.76
tt.Tx_DPD 0.95
per.NO.HeavyDD 0.43

No. of redistributions

  • Because the gamma distribution created some extremely high numbers we end up redistributin anyone who scored higher than 80 Standardard drinks. These two plots represent the number of mean times this occured at a given level of a simulation and the max times it occurred. The total number of cells in the data set was 80640 therefore even though this occurred it was only for a small number of the total data set.

Simulation Plots

  • I simulated 6 different levels of mode (range = 1.0 to 1.5) and 10 different levels of SD (range = 0.0 to 0.14)both have an impact on the mean of the drinking bias (X axis). There are 1000 simulations at each combination of levels of mode and SD (each point on the graph) for a total of 90000 simulations.
  • The black horizontal line on each point represents the value of the original COMBINE data set (N = 896; after removing 487 subjects due to missigness)

Drink Bias Plots

  • These are profiles of the drink bias distributions we used for each level of the simulation. These are based on a gamma distribution which is parameterized by two different parameters (shape and rate). However, one can convert a mode and SD of a distribution into a shape and rate parameters.
  • \(rate = ( mode + sqrt( mode^2 + 4*sd^2 ) ) / ( 2 * sd^2 )\)
  • \(shape = 1 + mode * rate\)

2017-04-24