Why is there no matching of timing data when trying to match partic... (2024)

% Load participant data

participant_data = readtable('945386.xlsx');

keyTap_rt = participant_data.keyTap_rt;

Rhythmsequence = participant_data.Rhythmsequence;

trials_thisN = participant_data.trials_thisN;

% Display Rhythmsequence for debugging

disp('Rhythmsequence:');

disp(Rhythmsequence);

% Load note onset data

note_onset_data = readtable('RPTiming.xls');

note_onset_Rhythmsequence = note_onset_data.Rhythmsequence;

note_onset_Notetiming = note_onset_data.Notetiming;

% Define durations for different audio stimuli

file_duration_data = { ...

'44_neutral_eighths.mp3', 33; ...

'44_prosody_eighths.mp3', 33; ...

'44_reverse_eighths.mp3', 33; ...

'78_neutral_eighths.mp3', 30; ...

'78_prosody_eighths.mp3', 30; ...

'78_reverse_eighths.mp3', 30 ...

};

% Initialize duration_map

duration_map = containers.Map();

% Populate duration_map

for i = 1:size(file_duration_data, 1)

filename = file_duration_data{i, 1};

duration = file_duration_data{i, 2};

duration_map(filename) = duration;

end

% Initialize matched data arrays

matched_keyTap_rt = [];

matched_Rhythmsequence = [];

matched_Notetiming = [];

% Store results for each trial

trial_results = struct('filename', {}, 'circular_mean', {}, 'circular_variance', {});

% Initialize a map to store results per audio file

audiofile_results = containers.Map('KeyType', 'char', 'ValueType', 'any');

% Define the list of relevant audio files

relevant_files = {'44_neutral_eighths.mp3', '44_prosody_eighths.mp3', '44_reverse_eighths.mp3', '78_neutral_eighths.mp3', '78_prosody_eighths.mp3', '78_reverse_eighths.mp3'};

% Iterate through each trial in the participant data from rows 7 to 30

for i = 7:30

% Get the full filename from the Rhythmsequence

full_filename = Rhythmsequence{i};

% Check if full_filename or keyTap_rt is empty

if isempty(full_filename) || isempty(keyTap_rt{i})

warning('Empty full_filename or keyTap_rt at index %d', i);

continue;

end

% Debugging statement to check the full_filename

disp('Full filename:');

disp(full_filename);

% Remove 'resources/' prefix if it exists

filename_without_prefix = strrep(full_filename, 'resources/', '');

% Debugging statement to check the filename without the prefix

disp('Filename without prefix:');

disp(filename_without_prefix);

% Check if the filename without prefix is in the list of relevant files

if ~ismember(filename_without_prefix, relevant_files)

continue; % Skip this trial if the file is not relevant

end

% Check if the filename without prefix is present in the duration map

if isKey(duration_map, filename_without_prefix)

% Retrieve duration from duration_map

duration_of_audio_stimulus = duration_map(filename_without_prefix);

% Debugging statement to check the duration

disp('Duration of audio stimulus:');

disp(duration_of_audio_stimulus);

else

warning('Duration not found for audio file: %s', filename_without_prefix);

continue; % Skip this trial if duration is not found

end

% Find the index of the matching Rhythmsequence in the note onset data

matching_idx = find(strcmp(note_onset_Rhythmsequence, full_filename));

% Match keyTap.rt with corresponding Notetiming based on Rhythmsequence

for j = 1:numel(matching_idx)

% Extract note onset timings for the current Rhythmsequence

note_onset_times = note_onset_Notetiming{matching_idx(j)};

% Convert note onset timings to numeric array

note_onset_times = str2double(strsplit(note_onset_times, ','));

% Check if the trial index is within the range of available keyTap_rt values

if i <= numel(keyTap_rt)

% Iterate over each note onset and find the closest tapping event

for k = 1:numel(note_onset_times)

% Find the index of the closest tapping timestamp to the current note onset timestamp

[~, closest_idx] = min(abs(keyTap_rt{i} - note_onset_times(k)));

% Match participant's tapping timestamp with note onset timestamp

matched_keyTap_rt = [matched_keyTap_rt; keyTap_rt{i}(closest_idx)];

matched_Rhythmsequence = [matched_Rhythmsequence; full_filename];

% For the first note, set the timestamp to 0.0

if k == 1

matched_Notetiming = [matched_Notetiming; 0.0];

else

matched_Notetiming = [matched_Notetiming; note_onset_times(k)];

end

end

% Convert timestamps to angles for current trial

angles = mod(2*pi * matched_keyTap_rt / duration_of_audio_stimulus, 2*pi);

% Calculate circular mean and variance for current trial

trial_circular_mean = circ_mean(angles);

trial_circular_variance = circ_var(angles);

% Perform Rayleigh test for circular mean

p_mean = circ_rtest(angles);

% Perform Rayleigh test for circular variance

p_variance = circ_rtest(angles);

% Store results for current trial

trial_results(end+1).filename = full_filename;

trial_results(end).circular_mean = trial_circular_mean;

trial_results(end).circular_variance = trial_circular_variance;

% Retrieve or initialize results for the current audio file

if isKey(audiofile_results, filename_without_prefix)

results = audiofile_results(filename_without_prefix);

else

results = struct('means', [], 'variances', []);

end

% Append current trial results to the audio file results

results.means(end+1) = trial_circular_mean;

results.variances(end+1) = trial_circular_variance;

% Update the audio file results in the map

audiofile_results(filename_without_prefix) = results;

end

end

end