Let Γ(·, λ) be smooth, i.e. C ∞ , embeddings from Ω onto Ω λ , where Ω and Ω λ are bounded domains with smooth boundary in the complex plane and λ varies in I = [0, 1]. Suppose that Γ is smooth on Ω × I and f is a smooth function on ∂Ω × I. Let u(·, λ) be the harmonic functions on Ω λ with boundary values f (·, λ). We show that u(Γ(z, λ), λ) is smooth on Ω × I. Our main result is proved for suitable Hölder spaces for the Dirichlet and Neumann problems with parameter. By observing that the regularity of solutions of the two problems with parameter is not local, we show the existence of smooth embeddings Γ(·, λ) from D, the closure of the unit disc, onto Ω λ such that Γ is smooth on D × I and real analytic at ( √ −1, 0) ∈ D × I, but for every family of Riemann mappings R(·, λ) from Ω λ onto D, the function R(Γ(z, λ), λ) is not real analytic at ( √ −1, 0) ∈ D × I.2010 Mathematics Subject Classification. 31A10, 45B05, 30C35, 35B30. Key words and phrases. Dirichlet and Neumann problems, Kellogg's theorem with parameter, domains with parameter, integral equations with parameter.