Occupants interact with buildings in various ways via their presence (passive effects) and control actions (active effects). Therefore, understanding the influence of occupants is essential if we are to evaluate the performance of a building. In this paper, we model the mean profiles and variability of occupancy variables (presence and actions) separately. We will use a multi-variate Gaussian distribution to generate mean profiles of occupancy variables, while the variability will be represented by a multi-dimensional time series model, within a framework for a meta-analysis that synthesizes occupancy data gathered from a pool of buildings. We then discuss variants of occupancy models with respect to various outcomes of interest such as HVAC energy consumption and peak demand behavior via a sensitivity analysis. Results show that our approach is able to generate stochastic occupancy profiles, requiring minimum additional input from the energy modeler other than standard diversity profiles. Along with the meta analysis, we enable the generalization of previous research results and statistical inferences to choose occupancy variables for future buildings. The sensitivity analysis shows that for aggregated building energy consumption, occupant presence has a smaller impact compared to lighting and appliance usage. Specifically, being accumulatively 55% wrong with regard to presence, only translates to 2% error in aggregated cooling energy in July and 3.6% error in heating energy in January. Such a finding redirects focus to the accurate estimation of lighting and appliance usage for a better prediction of aggregated energy consumption. Furthermore, it proves that accurate knowledge of the mean profiles is sufficient, that is, stochastic occupancy models do not play a significant role in the prediction of aggregated consumption in a conventional office building where the interaction between the operation of building systems and the spatial and temporal variability of occupancy is weak. When it comes to peak demand behavior, occupancy variability should be taken into account, as static profiles are not able to produce adequate estimates of power duration probabilities close to the power peak. (C) 2016 Elsevier Ltd. All rights reserved.