function grid = Stats(X, sex, age, col, grid)

    X = X(strcmp(X.ACTUALSEX, sex), :);
    if strcmp(age, 'AHY') || strcmp(age, 'HY')
        X = X(strcmp(X.Age, age), :);
    end
    
    grid{1, col} = age;
    grid{2, col} = sex;
    grid{3, col} = sprintf(' n = %d', size(X,1));

    grid{4, 1} = 'Culmen length (mm)';
    grid{4, col} = sprintf('%.2f ± %.2f', mean(X.Cul), std(X.Cul));
    grid{5, 1} = '';
    grid{5, col} = sprintf('(%.1f–%.1f)', min(X.Cul), max(X.Cul));

    grid{6, 1} = 'Tarsometatarsus length (mm)';
    grid{6, col} = sprintf('%.2f ± %.2f', mean(X.Tar), std(X.Tar));
    grid{7, 1} = '';
    grid{7, col} = sprintf('(%.1f–%.1f)', min(X.Tar), max(X.Tar));

    grid{8, 1} = 'Toe length (mm)';
    grid{8, col} = sprintf('%.2f ± %.2f', mean(X.Toe), std(X.Toe));
    grid{9, 1} = '';
    grid{9, col} = sprintf('(%.1f–%.1f)', min(X.Toe), max(X.Toe));

    grid{10, 1} = 'Weight (gm)';
    grid{10, col} = sprintf('%.2f ± %.2f', mean(X.Wt), std(X.Wt));
    grid{11, 1} = '';
    grid{11, col} = sprintf('(%.1f–%.1f)', min(X.Wt), max(X.Wt));
end

