SUMMARYIn this paper, we will deal with interesting example of natural risk modeling -namely of indoor radon concentration, needed for proper exposure assessment and appreciation of quality of preventive measures taken. First, we will illustrate, how the traditional view based on fixed (i.e., time-invariant) coefficient models can be misleading. Then we formulate a flexible (nonparametric) regression dynamic model offering a doable alternative. Situation is not entirely standard here, as the smoothed and weight-giving variables are different. Nevertheless, estimation for our local regression formulation is easily doable even for rather large data via extension of standard local smoothing that we describe. We illustrate on data coming from a rare intensive measurement campaign in an occupied house with simultaneous measurements taken in different rooms. Our model has nice physical interpretation of its parts. We also illustrate how even such a simple model can produce behavior that throws some light on radon experts' discussions about natural radon concentration circadial movement phase. This is because our model acts as a linear but time-varying filter that can change not only amplitude but also phase between different rooms of the same house, suggesting that there might not be a universal phase of (e.g., daily) radon variation. Instead, the phase might depend on how close a particular indoor location is to the radon sources. Finally, we present some ideas about how our model can be expanded in future to cover more complicated situations and settings.