Welcome!

Sign in with your CogSci, Facebook, Google, or Twitter account.

Or register to create a new account.

We'll use your information only for signing in to this forum.

Supported by

Sample by sample time series analysis with generalTestBF

JasperJasper Posts: 1
edited February 10 in JASP & BayesFactor

Hi all,

I have a question about the use of generalTestBF from the BayesFactor package for R. I believe I followed all the steps provided on http://bayesfactorpcl.r-forge.r-project.org/#generalTestBF correctly, but I am not entirely sure if I am interpreting the results correctly.

For the analysis of the time series of a continuous dependent variable (pupil size), I have constructed general linear mixed models (GLMM) using generalTestBF for each sample in the time series. In the GLMM, the dependent variable is modeled as a function of two 2-level factors with a random intercept for each subject

m <- generalTestBF( 
pupil ~ task * luminance + subject,
data = tmpdata,
whichRandom = 'subject',
whichModels = 'all',
iterations = 10000,
neverExclude = 'subject' 
)

To find evidence for an effect of any of the independent variable (IV) I summed the Bayes factors of all the models containing that IV, and divided that by the sum of Bayes factors of all the models that did not contain that IV.

After running the above model for all samples in the time series, you end up with, what I interpreted as, a Bayes Factor time series. Example plot below.

To interpret this Bayes factor time series, I picked the commonly used BF > 3 to indicate epochs with evidence in favor of an effect of the IV on the dependent variable (in the example plot illustrated by the black horizontal bar).

Is my interpretation of this output correct? Or am I missing an important step in the analysis?

Jasper

Comments

  • EJEJ Posts: 409

    Hi Jasper (great name, BTW)

    This is really a question for Richard Morey. You could try to private-message him perhaps.

    Cheers,
    E.J.

  • sebastiaansebastiaan Posts: 2,856

    There's much bigger issues in the world, I know. But I first have to take care of the world I know.
    cogsci.nl/smathot

Sign In or Register to comment.