Selecting Data for mouse measures

Comments

• Pascal

  September 2018 edited September 2019

  Hi!

  that is in principle possible, but the specific implementation depends on the measure you want to use.

  The following code gives an example how you can calculate the Euclidean distance traveled from the beginning of the trajectory until timestamp 600 (ms):

  # Load mousetrap library
  library(mousetrap)
  
  # Calculate mouse-tracking derivatives for each position (distance, velocity, acceleration)
  mt_example <- mt_derivatives(mt_example)
  
  # Calculate standard mouse-tracking measures
  mt_example <- mt_measures(mt_example)
  
  # Demonstrate that sum of distances corresponds to total_dist measure
  table(rowSums(mt_example$trajectories[,,"dist"],na.rm=TRUE)-mt_example$measures$total_dist)
  
  # Set distances for timestamps beyond 600 ms to 0
  mt_example$trajectories[,,"dist"][mt_example$trajectories[,,"timestamps"]>600] <- 0
  
  # Calculate sum of distances for timestamps <= 600 ms as a new measure in mt_example$measures
  mt_example$measures$total_dist_600 <- rowSums(mt_example$trajectories[,,"dist"],na.rm=TRUE)
  

  Hope this helps. If you have any additional questions, feel free to ask!

  Best,

  Pascal
• a.g.leontiev

  September 2018 edited September 2018

  Thank you so much! It really does help!

  I have one additional question: how do I summarize these measures by subjID (and possibly, by some other factors)?

  thank you!
• a.g.leontiev

  September 2018

  This is how I did it:

  average_measures <- mt_aggregate_per_subject(
    mt_data,
    use="measures",
    use_variables=c("total_dist_600"),
    subject_id="subjID"
  )
  
• Pascal

  September 2018

  You can use the mt_aggregate_per_subject function for this. If you specify a condition variable in use2_variables, the measures will be summarized per participant and condition.

  average_measures <- mt_aggregate_per_subject(
    mt_example,
    use="measures",
    use_variables=c("total_dist_600 "),
    use2_variables="Condition",
    subject_id="subjID"
  )
  
  

  Best,

  Pascal
• a.g.leontiev

  September 2018

  thank you!
• a.g.leontiev

  October 2018

  I have another question. Should I want to calculate other measures, like acceleration and velocity, will setting trajectories to zero be sufficient? In other words, if i want to calculate say velocity after certain time period, do I only need to change total_dist to vel_max?

  Thank you!
• Pascal

  October 2018

  You would have to make a few changes:

  # Set velocities for timestamps beyond 600 ms to 0
  mt_example$trajectories[,,"vel"][mt_example$trajectories[,,"timestamps"]>600] <- 0
  
  # Calculate maximum velocity
  mt_example$measures$vel_max_600 <- apply(mt_example$trajectories[,,"vel"],1,max,na.rm=TRUE)
  

  Instead of the second line, you could also run mt_measures afterwards. It will automatically calculate the total_dist and vel_max based on the modified values, so the resulting values for total_dist an vel_max will be the ones for the first 600 ms.