We consider the estimation of smooth regression functions in a class of conditionally parametric co-variate-
response models. Independent and identically distributed observations are available from the distribution of
(Z,X), where
Z is a real-valued co-variate with some unknown distribution, and the
response X conditional on
Z is distributed according to the density
p(·,ψ(Z)), where
p(·,θ) is a one-parameter exponential family. The function
ψ is a smooth
monotone function. Under this formulation, the regression function
E(X|Z) is
monotone in the co-variate
Z (and can be expressed as a one–one function of
ψ); hence the term “
monotone response model”. Using a penalized least squares approach that incorporates both monotonicity and smoothness, we develop a scheme for producing smooth
monotone estimates of the regression function and also the function
ψ across this entire class of
models. Point-wise asymptotic normality of this estimator is established, with the rate of convergence depending on the smoothing parameter. This enables construction of Wald-type (point-wise) as well as pivotal confidence sets for
ψ and also the regression function. The methodology is extended to the general heteroscedastic model, and its asymptotic properties are discussed.