function plots(data_boys, data_girls, label) %1. age 2-4. kth moment of syllables
  
    boys=data_boys;
    girls=data_girls;
    
    for k=1:3
        [r,p]=corr([boys(:,1),boys(:,k+1)],'Type','Pearson'); 
        rp_m(1,k)=r(1,2);
        rp_m(2,k)=p(1,2);
        [r,p]=corr([girls(:,1),girls(:,k+1)],'Type','Pearson'); 
        rp_f(1,k)=r(1,2);
        rp_f(2,k)=p(1,2);
    end
    
    figure;
    subplot 311;
    plot(boys(:,1),boys(:,2),'bx'); hold on; plot(girls(:,1),girls(:,2), 'ro'); 
    p_pab=polyfit(boys(:,1),boys(:,2),1); y_pab=polyval(p_pab,boys(:,1));
    plot(boys(:,1),y_pab,'b');
    p_pag=polyfit(girls(:,1),girls(:,2),1); y_pag=polyval(p_pag,girls(:,1));
    plot(girls(:,1),y_pag,'r'); hold off;
    title(['PA: Corr boys: ', num2str(rp_m(1,1), '%.3f'), '; Corr girls: ', num2str(rp_f(1,1),'%.3f')]); 
    xlabel('age'); ylabel(label);

    subplot 312;
    plot(boys(:,1),boys(:,3),'bx'); hold on; plot(girls(:,1),girls(:,3), 'ro'); 
    p_tab=polyfit(boys(:,1),boys(:,3),1); y_tab=polyval(p_tab,boys(:,1));
    plot(boys(:,1),y_tab,'b');
    p_tag=polyfit(girls(:,1),girls(:,3),1); y_tag=polyval(p_tag,girls(:,1));
    plot(girls(:,1),y_tag,'r'); hold off;
    title(['TA: Corr boys: ', num2str(rp_m(1,2), '%.3f'), '; Corr girls: ', num2str(rp_f(1,2),'%.3f')]); 
    xlabel('age'); ylabel(label);

    subplot 313;
    plot(boys(:,1),boys(:,4),'bx'); hold on; plot(girls(:,1),girls(:,4), 'ro'); 
    p_kab=polyfit(boys(:,1),boys(:,4),1); y_kab=polyval(p_kab,boys(:,1));
    plot(boys(:,1),y_kab,'b');
    p_kag=polyfit(girls(:,1),girls(:,4),1); y_kag=polyval(p_kag,girls(:,1));
    plot(girls(:,1),y_kag,'r'); hold off;
    title(['KA: Corr boys: ', num2str(rp_m(1,3), '%.3f'), '; Corr girls: ', num2str(rp_f(1,3),'%.3f')]); 
    xlabel('age'); ylabel(label);
end
