In this paper we develop statistical models for bankruptcy prediction of Italian firms in the limited liability sector, using annual balance sheet information. Several issues involved in default risk analysis are investigated, such as the structure of the data-base, the sampling procedure and the influence of predictors. In particular we focus on the variable selection problem, comparing innovative techniques based on shrinkage with traditional stepwise methods. The predictive performance of the proposed default risk model has been evaluated by means of different accuracy measures. The results of the analysis, carried out on a data-set of financial ratios expressly created from a sample of industrial firms annual reports, give evid...