# Appropriate technique to extract the response function of a time series to two (or more) irregularly recurring stimuli

Say I poll a teenager every day about his overall happiness. Perhaps this has some long term trend, perhaps it has some short term trend. For sure, it is a noisy variable. Over the relevant time period, the teenager experiences two kinds of stimuli at known, but irregularly spaced, times: he has to take statistics exams and he gets to go to school dances.

I would hypothesize that the lead up to an exam should have a continuous impact on his happiness. Then, perhaps, there is a jump in happiness at the instant of the exam. And then there is a continuous impact after the exam. (I would expect the lead-in and exit patterns to be reasonably well contained to the time period around the "exam event").

And I would hypothesize that there is an analogous (but different in detail) kind of pattern around school dances.

So after a year or two, I'll have:

• Happiness measurements at each time
• A list of times at which exams are administered
• A list of times at which dances occurred
• Some beliefs about the continuity and localization of the expected response of happiness to the given types of stimuli as described above.

What would be an appropriate statistical technique to estimate the "shape" of the happiness response function in the vicinity of an exam and, separately, in the vicinity of a dance?

I am willing to assume that the observed impact on happiness of exams and dances should be the sum of the impacts that would be observed in the presence of each separately. However, in practice, I expect that a substantial portion of my observations will be simultaneously impacted by more than one stimulus. (I.e. dances and exams are not sufficiently separated in time to guarantee that any given happiness measurement is impacted by only a single event.)