Update hmrS_GLM.m

        end

        % update time on stim object

        stimRuns{i+1}(j).data(:,1) = stimRuns{i+1}(j).data(:,1) + max_t(i) + dt_foo;

        stimRuns{i+1}(j).states(:,1) = stimRuns{i+1}(j).states(:,1) + max_t(i) + dt_foo;

        if stimRuns{i+1}(j).data(1,2) ~= 0

            foo1 = stimRuns{i+1}(j).data;

            foo2 = stimRuns{i+1}(j).states;

            foo1(1,1) = stimRuns{i+1}(j).data(:,1) + max_t(i) + dt_foo;

            foo2(1,1) = stimRuns{i+1}(j).states(:,1) + max_t(i) + dt_foo;

            % concatenate

        foo_stim_data = cat(1,foo_stim_data,stimRuns{i+1}(j).data);

        foo_stim_states = cat(1,foo_stim_states,stimRuns{i+1}(j).states);

            foo_stim_data = cat(1,foo_stim_data,foo1);

            foo_stim_states = cat(1,foo_stim_states,foo2);

            clear foo1 foo2

        end

    end