Skip to main content

Advertisement

Log in

Identifying Berlin’s land value map using adaptive weights smoothing

  • Original Paper
  • Published:
Computational Statistics Aims and scope Submit manuscript

Abstract

We use adaptive weights smoothing (AWS) of Polzehl and Spokoiny (J R Stat Soc Ser B 62:335–354, 2000; Ann Stat 31:30–57, 2003; Probab Theory Relat Fields 135:335–362, 2006) to estimate a map of land values for Berlin, Germany. Our data are prices of undeveloped land that was transacted between 1996 and 2009. Even though the observed land price is an indicator of the respective land value, it is influenced by transaction noise. The iterative AWS applies piecewise constant regression to reduce this noise and tests at each location for constancy at the margin. If not rejected, further observations are included in the local regression. The estimated land value map conforms overall well with expert-based land values. Our application suggests that the transparent AWS could prove a useful tool for researchers and real estate practitioners alike.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11

Similar content being viewed by others

Notes

  1. See Polzehl and Spokoiny (2006). For simplicity, we nonetheless refer to the approach as AWS throughout the paper.

  2. Fitting a third degree polynomial for the distance to the CBD gives such a decreasing function with an \(R^2=0.2345\).

  3. Scatterplots of the estimated AWS land values, \(\hat{\theta }(\mathbf {x}_i)\), against the land prices, \(y_i\), (not reported) as well as kernel density estimates of the estimated AWS residuals, \(\hat{\epsilon }_i\), (not reported) indicate that the assumption of normally distributed and homoscedastic (log) land prices is approximately satisified in our application below.

  4. Obviously, if only the observation itself was used in step 0, \(A_i^0=1\), and Eq. 11 becomes Eq. 9.

  5. While the exact sampling distribution of \(T^k_{ij}\) can only be derived in iteration step \(k=1\) if the bandwidth \(h^0\) is very small, the \(\chi ^2_1\)-distribution may still be a good approximation in every iteration step \(k\) under the assumptions of normally distributed, homoscedastic errors.

  6. The condition involves the probability that the Kullback Leibler divergence between the adaptive AWS estimate and the globally constant \(\theta \) is bounded and that the bound does not increase in the iteration process for a given propagation level \(\epsilon .\) Becker and Mathé (2013) also propose a method for estimating this probability from the simulated data.

  7. We use the sequence \(\mathbf {h}=\left\{ 1, 2, 3, 4, 5, 7, 9, 11, 14, 18, 22, 28, 35, 44, 55, 69, 86, 108, 135\right\} \) to increase \(h^{k}\) in the \(k^{*}=19\) iteration steps.

  8. \(\lambda =3.8415\) corresponds to the \(5\,\%\) percentile of the \(\chi ^2_1\)-distribution.

References

  • Ahlfeldt GM, Maennig W (2010) Impact of sport arenas on land values: evidence from Berlin. Ann Reg Sci 44:205–227

    Article  Google Scholar 

  • Becker S (2014) The propagation–separation approach. Theoretical study and application to magnetic resonance imaging, PhD thesis, Humboldt-Universität zu Berlin

  • Becker S, Mathé P (2013) A different perspective on the propagation–separation approach. Electr J Stat 7:2702–2736

    Article  MATH  Google Scholar 

  • Bostic RW, Longhofer SD, Redfearn CL (2007) Land leverage: decomposing home price dynamics. Real Estate Econ 35:183–208

    Article  Google Scholar 

  • Bryan KA, Satre P-DG (2009) Semiparametric estimation of land price gradients using large data sets. Fed Reserve Bank Richmond Econom Q 95:53–74

    Google Scholar 

  • Diamond DB (1980) The relationship between amenities and urban land prices. Land Econom 56:21–32

    Article  Google Scholar 

  • Fan J, Marron JS (1994) Fast implementations of nonparametric curve estimators. J Comput Graph Stat 3:35–56

    Google Scholar 

  • Kennedy P (1983) Logarithmic dependent variables and prediction bias. Oxf Bull Econ Stat 45:389–392

    Article  Google Scholar 

  • Polzehl J (2014) Adaptive weights smoothing. R package version 1.9-4, http://CRAN.R-project.org/package=aws

  • Polzehl J, Spokoiny V (2000) Adaptive weights smoothing with applications to image restoration. J R Stat Soc Ser B 62:335–354

    Article  MathSciNet  Google Scholar 

  • Polzehl J, Spokoiny V (2003) Image denoising: pointwise adaptive approach. Ann Stat 31:30–57

    Article  MATH  MathSciNet  Google Scholar 

  • Polzehl J, Spokoiny V (2006) Propagation–separation approach for local likelihood estimation. Probab Theory Relat Fields 135:335–362

    Article  MATH  MathSciNet  Google Scholar 

  • Polzehl J, Spokoiny V (2008) Structural adaptive smoothing by propagation–separation methods. In: Chen C, Härdle W, Unwin A (eds) Handbook of data visualization. Springer, Berlin, pp 471–492

    Chapter  Google Scholar 

  • Stone CJ (1985) Additive regression and other nonparametric models. Ann Stat 13:689–705

    Article  MATH  Google Scholar 

Download references

Acknowledgments

We have benefited from comments received at the Applicable Semiparametrics Conference in Berlin 2013 and from two anonymous referees. Financial support from the Deutsche Forschungsgemeinschaft, CRC 649 Economic Risk, is gratefully acknowledged. The usual disclaimer applies.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Axel Werwatz.

Appendix

Appendix

We adjust the land prices as follows. First, we run the regression

$$\begin{aligned} p_{j,t} = \mathbf {d}_j\varvec{\alpha } + \mathbf {x}_{j}\varvec{\beta } + \mathbf {z}_{j}\varvec{\gamma } + \epsilon _{j,t}, \end{aligned}$$
(15)

where \(p_{j,t}\) is the log price of site \(j\) transacted in quarter \(t\). The column vector \(\mathbf {d}_{j}\) has \(T\) elements: the first element for the overall constant is one, the period \(t\) element is one if site \(j\) was transacted in this period, zero otherwise. 1996Q1 is the omitted reference period. The vector \(\mathbf {x}_{j}\) contains binary indicators for unusual features of the site and for unusual aspects of the business dealings. The vector \(\mathbf {z}_{j}\) contains binary indicators for Berlin’s 96 administrative sub-districts (Ortsteile). The site will be located in one of these sub-districts. The vector \(\mathbf {z}_j\) also contains a binary indicator for site location adjacent to a lake or the bank of a river. Finally, the vector contains binary indicators for site’s location rating. This rating comes from Berlin’s Senate Department for Urban Development and rates natural amenities, the quality of existing buildings, access to public transport and shopping facilities with in the neighborhood. The rating for a site takes one of four values: low, medium, high, very high. The variables in \(\mathbf {z}_j\) control crudely for location effects. Without the inclusion, the estimates of \(\varvec{\alpha }\) and \(\varvec{\beta }\) may suffer from omitted variable bias. Table 4 presents least squares estimates of the model in Eq. 15. The in-sample fit, as measured by the \(R^2\), is reasonably good. Except for the coefficient for ground monument, all coefficients are statistically significant at the 5 % significance level. The signs of the point estimates, as well as their magnitude, are plausible.

Table 4 Parameter estimates for land price adjustment model

Second, given the coefficient estimates, we compute the adjusted log real land price as

$$\begin{aligned} p_{j} =p_{j,t} - (\mathbf {d}_j-\mathbf {d}_{b})\widehat{\varvec{\alpha }} - \mathbf {x}_{j,t}\widehat{\varvec{\beta }}. \end{aligned}$$
(16)

The first entry of \(\mathbf {d}_{b}\) is one, the entries for the four quarters of the year 2009 are 0.25 each, the remaining entries are zero. The term in brackets in Eq. 16 converts prices to the base year 2009. The estimated value of \(\mathbf {z}_{j}\varvec{\gamma } \) is not considered for \(p_j\), because it enters Eq. 15 only to prevent bias. The resulting \(p_{j}\) is in real terms and adjusted for unusual circumstances of the site. Using it in our analysis puts us on an equal footing with the land price information used by local surveyors to produce the BRW. The summary statistics for prices in natural scale in Table 1 are computed using \(P_j=\exp \!\left\{ p_{j} + 0.5 \widehat{\sigma }^2_{\epsilon }\right\} \), where \(\widehat{\sigma }^2_{\epsilon }\) is the estimated variance of the error term in Eq. 15 (Kennedy 1983).

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kolbe, J., Schulz, R., Wersing, M. et al. Identifying Berlin’s land value map using adaptive weights smoothing. Comput Stat 30, 767–790 (2015). https://doi.org/10.1007/s00180-015-0559-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s00180-015-0559-9

Keywords

Navigation