calls_tonality_analysis.m

% info = inputdlg(["Start day","Interval","End day"], 'double');
% 
% start_day = str2num(cell2mat(info(1)));
% interval = str2num(cell2mat(info(2)));
% end_day = str2num(cell2mat(info(3)));

start_day = 4;  %First recordings on P6
interval = 2;   %Record every 2 days
end_day = 18;   %Last recordings on P18
n_days = length(start_day:interval:end_day);

%% initialize variables
m_tonality = [];
se_tonality = [];

%% Get daily means for all variables

subj_means_tonality = [];
subj_sds = [];
subj_names = unique(T.Var2);

for k=1:n_days
    day = (start_day-interval)+(k*interval);
    indexes = T.Var1==int2str(day);
    subT = T(indexes,:);
    
    m_tonality(end+1) = mean(subT.Tonality);
    for n=1:length(subj_names)
        indiv_indexes = subT.Var2==subj_names(n);
        subsubT = subT(indiv_indexes,:);
        subj_means_tonality(n,k) = mean(subsubT.Tonality);
        subj_sds(n,k) = std(subsubT.Tonality);
    end
    
    m_tonality(end) = mean(subT.Tonality); %Mean tonality for all calls that day
end

se_tonality = std(subj_sds, 'omitnan');
daily_subject_means = mean(subj_means_tonality, 'omitnan');  %Mean across subjects per day
subject_7day = mean(subj_means_tonality, 2, 'omitnan');      %Mean across days per subject
%% Figures

[p,tab,stats] = kruskalwallis(subj_means_tonality);
hold on; box off;
xticklabels([start_day:interval:end_day]);
title('Tonality');
xlabel('Postnatal Day');
ylabel('Mean subject USV tonality (dB/kHz)');

figure(3); hold on;
errorbar(daily_subject_means,se_tonality,'k-','LineWidth',1)
plotSpread(subj_means_tonality)
%plot(daily_subject_means,'k-','LineWidth',1)
box off; zoom out;
title('Tonality');
xlabel('Postnatal Day');
ylabel('Mean subject USV tonality (dB/kHz)');
xticklabels([start_day:interval:end_day]);


figure(6)
hold on;
colormap = parula;
colormap12 = colormap(1:256/12:256,:)
subj_means_tonality_sorted = sortrows(subj_means_tonality,1,'ascend')

for m=1:12
    plot(subj_means_tonality_sorted(m,:),'Color',colormap12(m,:),'LineWidth',1.5);
end
title('Tonality');
xlabel('Postnatal Day');
ylabel('Mean subject USV tonality (dB/kHz)');
xticklabels([start_day:interval:end_day]);