get_combined_weights_downsample.m

function [ws,tr] = get_combined_weights_downsample(stats,covariate,factor)
    ws=[stats.ws];
    if ~isfield(ws,covariate)
        error('%s is not a field of the weight structure.',covariate);
    end
    tmp = [ws.(covariate)];
    ws = tmp(1);
    offset=0.01;
    ex=-0.5;
    if contains(covariate,'click')
        idx = logspace(ex,log10(ws.tr.lims(2)+offset+10^ex),round(size(ws.data,2)/factor));
        idx= idx-(idx(1)-offset);
        [~,idx ] = min (abs(idx-buildGLM.get_tr(ws.tr)'));
    else
        idx = linspace(ws.tr.lims(1),ws.tr.lims(2),round(numel(ws.data)/factor));
        [~,idx ] = min (abs(idx-buildGLM.get_tr(ws.tr)'));        
    end
    ws.data = cat(1,tmp.data);
    ws.data = ws.data(:,idx);
    tr = buildGLM.get_tr(ws.tr);
    tr = tr(idx);
    if contains(covariate,'click')
        %ws.data = ws.data*8;
    end

    nx = numel(tmp);
    ws.mean_data = mean(ws.data);
    ws.se = std(ws.data)./sqrt(nx);
end