function [rms, nees] = compute_rms_nees_ave( rmsall, neesall )
% compute rms and nees from multiple tests

nsteps = size(rmsall.position, 2);
ntests = size(rmsall.position, 1);

% compute rms
rmspos_sum = sum(rmsall.position.^2, 1);
rms.position = sqrt(rmspos_sum/ntests);

rmsori_sum = sum(rmsall.orientation.^2, 1);
rms.orientation = sqrt(rmsori_sum/ntests);

% compute nees
neesall_pose = neesall.pose(:, 2:end);
neesall_orientation = neesall.orientation(:, 2:end);

neespose_sum = sum(neesall_pose, 1);
nees.pose = neespose_sum/ntests;

neesorientation_sum = sum(neesall_orientation, 1);
nees.orientation = neesorientation_sum/ntests;