Suppose $N$ experiments can be made in varying conditions. Each of them yields an estimate $f_i$ of a continuous (and, if necessary, positive) function of x over some interval. Experiment $i$ is repeated a number of times $n_i$, and each time it yields a function $f_i^j$. I am given the following data $$\forall x \quad \forall i \quad {n_i,f_i(x),\sigma_i(x)}$$ where $$f_i(x) = \frac{1}{n_i}\sum_{j=1}^{n_i} f_i^j(x) $$ $$\sigma_i^2(x) = \frac{1}{n_i}\sum_{j=1}^{n_i} (f_i^j(x)-f_i(x))^2 $$
Suppose experiment $1$ is the reference experiment. I would like to know whether some or all of the $f_i$ agree with function $f_1$. If they do not agree, I would like to know the region in which there is disagreement. Note however that since the experiments were made in different conditions, it is expected that the variance functions do not agree, and have to be treated as different for each experiment.
For now, I was planning on doing a two-way ANOVA on each $x$ value followed by a two-sample $t$-test if necessary to determine the disagreeing experiments. However, I was looking for a method that treats the function as a whole, heteroscedastic, and that is bayesian if possible. I could not find anything. Literature references appreciated.