The error you're encountering suggests that the groupsummary function is expecting a table or a dataset array, but it's receiving a standard numeric matrix instead. To resolve this issue, try converting matrix into a table format before using groupsummary.
To do this follow the given below steps:
- Convert the filtered data matrix into a table format.
- Apply "groupsummary" on the table.
The revised code for the same is given below:
% Sample data
data = [
8 4 340; 8 4 328; 8 3 218; 8 4 338; 8 3 213;
8 4 328; 8 3 254; 8 4 323; 8 4 340; 8 3 273;
9 3 580; 9 4 363; 9 4 371; 9 3 374; 9 3 383;
9 3 302; 9 4 406; 9 3 390; 9 3 380; 9 3 366;
9 4 468
];
% Extract columns
subjects = data(:, 1);
trialcodes = data(:, 2);
latencies = data(:, 3);
% Define a function to remove outliers
removeOutliers = @(latencies) latencies(abs(latencies - mean(latencies)) <= 0.5 * std(latencies));
% Group by subject and trialcode
[G, ~] = findgroups(subjects, trialcodes);
% Apply the function to each group
cleanedData = splitapply(@(latencies) {removeOutliers(latencies)}, latencies, G);
% Reconstruct the data matrix without outliers
filteredData = [];
for i = 1:length(cleanedData)
if ~isempty(cleanedData{i})
group = unique(data(G == i, 1:2), 'rows');
filteredData = [filteredData; repmat(group, size(cleanedData{i}, 1), 1), cleanedData{i}];
end
end
% Convert filtered data to a table
filteredTable = array2table(filteredData, 'VariableNames', {'Subject', 'TrialCode', 'Latency'});
% Calculate the mean latency for each subject and trialcode
averageLatencies = groupsummary(filteredTable, {'Subject', 'TrialCode'}, 'mean', 'Latency');
% Display the result
disp(averageLatencies);
To understand more about "groupsummary" refer to the following documentation:
- https://www.mathworks.com/help/matlab/ref/double.groupsummary.html
Hope this helps!
