Beauchamp:CIMS
Causal Inference of Multisensory Speech Model
NB: This code is in active development
A copy of this help file is included in the code pack as README.txt
If you run into trouble with any step, please contact me: john dot magnotti at gmail dot com. If the model fitting fails to converge for your dataset, you may need to send me at least a portion of the data so I can replicate the error.
If you find this code useful, please cite our work:
Magnotti, J. F., Ma, W. J., & Beauchamp, M. S. (submitted). Causal Inference in Multisensory Speech Perception.
Fitting the CIMS model to your data
System setup
- Install GNU R
- Download the CIMS CODE PACK (link this to Dropbox)
- Extract the zip file to your Desktop/ or other preferred location
Data setup
The model code assumes the data are stored in a matrix format with rows as subjects and each column a separate asynchrony. The first row is used as labels for the columns. Each cell stores the number of times the subject judged the stimulus at the given asynchrony as synchronous. If there are 17 subject and 15 asynchronies, the file will have 18 rows (1st row is header row) and 15 columns. For multi-condition experiments, place each asynchrony/condition combination in a separate column. Asynchronies within a condition must be contiguous. If the study has 4 conditions and 15 asynchronies, then the columns 1-15 will be treated as condition 1, 16-30 as condition 2, and so on. See data.csv in the code pack for the data used in Magnotti, Ma, & Beauchamp.
If you are unfamiliar with R, the best approach is to run through all the model building steps using the included data.csv file, and then try with your own data.
Program setup
We need to ensure R can find the data and code files
- Launch R
- Open the file fit_models.R: File-> Open Document
- We need to make 4 changes before running the code. See the comments in the code file for additional direction
- Line 5: Set the path to be the location of the downloaded files. If you extracted the code pack to your desktop, the path may already be correct
- Line 11: Set the location of the data to be fit
- Line 14: Set the value of max_count to be the total number of trials at each asynchrony for each condition.
- Line 17: Set the value of the asynchronies used in each condition. The order of the asynchronies must match the order in the data file. For multiple condition experiments, list the asynchronies only once.
- Run the setup code to make sure there are no errors
- Highlight lines 1 through 17 using the mouse
- Execute the code by using the R menu: Edit -> Execute
Fitting the model
Highlight and execute each of the following lines in turn
cl = makeCluster(detectCores())
cims.model = cims(n.reps=1000)
gauss.model = gauss()
Model Parameters
- The resulting parameters for each model are saved to cims_out.csv and gauss_out.csv.
- The predicted values for each model are saved to cims_predicted.csv and gauss_predicted.csv.
Model Comparison
Batch comparison within and across conditions
within_condition_tests(cims.model, gauss.model) across_condition_tests(cims.model, gauss.model)
Advanced Functions
We include here helpful functions for further data exploration. The functions assume you have run all the code in fit_models.R. Intrepid users are encouraged to let the the source be their guide.
Extract the log Likelihood from each model
logLik(cims.model)
logLik(cims.model, separate=T)
logLik(gauss.model)
logLik(gauss.model, separate=T)
To plot model fits
ebar_plot(asyncs, count_mat/max_count)
add_fit_line(cims.model, col='orange')
add_fit_line(gauss.model, col='blue')
To obtain the raw predicted values:
cims.predicted = predict(cims.model)
gauss.predicted = predict(gauss.model)
Copyright/Licensing
http://i.creativecommons.org/l/by-sa/3.0/88x31.png
This work is licensed under a Creative Commons Attribution-ShareAlike 3.0 Unported License.