In this paper we consider a heteroscedastic transformation model of the form Λν(Y ) = m(X)+σ(X)ϵ, where Λν belongs to a parametric family of monotone transformations, m(.) and σ(.) are unknown but smooth functions, ϵ is independent of the d-dimensional vector of covariates X, E(ϵ) = 0 and Var(ϵ) = 1. We consider the estimation of the unknown components of the model, ν, m(.), σ(.), and the distribution of ϵ, and we show the asymptotic normality of the proposed estimators. We propose tests for the validity of the model, and establish the limiting distribution of the test statistics under the null hypothesis. A bootstrap procedure is proposed to approximate the critical values of the tests. We carried out a simulation study to verify the small sample behavior of the proposed estimators and tests, and illustrate our method with a dataset.