Skip to contents

Compare the average cross power distribution of two SplicedViews using sampling on each segment

Usage

compare_avg_cross_power2(
  sv1,
  sv2,
  name1,
  name2,
  num_samples,
  columns,
  show_plot = TRUE
)

Arguments

sv1

SplicedView object.

sv2

SplicedView object.

name1

name for first object.

name2

name for second object.

num_samples

number of samples to draw from segments.

columns

column names in the data e.g. c('Nose_x', 'Nose_y').

show_plot

show the plot?

Value

list of two data.frames containing the sampled data.

Examples

r <- get_sample_recording()
d1 <- get_duration_annotation_data(r)
fv_list <- get_filtered_views(r, data_points = "Nose", n = 41, p = 3)
jv <- get_joined_view(fv_list)

# only one relevant section for sample data
splicing_smile_df <- splice_time(d1, tier ='INTERACTION',
  comments = 'Mutual look and smile')
sv_duration_smile <- get_spliced_view(jv, splicing_df = splicing_smile_df)

splicing_alap_df <- splice_time(
  d1, tier = 'FORM', comments = 'Alap'
)
sv_duration_alap <- get_spliced_view(jv, splicing_df = splicing_alap_df)

sample_list <- compare_avg_cross_power2(
sv_duration_smile, sv_duration_alap, 'Smile', 'Alap', num_samples = 100,
  columns = c("Nose_x_Central_Sitar", "Nose_y_Central_Sitar"))
#> Warning: NaNs produced
#> 
  |                                                                            
  |                                                                      |   0%
#> Warning: NaNs produced
#> 
  |                                                                            
  |======================================================================| 100%
#> 
  |                                                                            
  |                                                                      |   0%
  |                                                                            
  |======================================================================| 100%