We use type Ia supernovae (SN Ia) data in combination with recent baryonic acoustic oscillations (BAO) and cosmic microwave background (CMB) observations to constrain a kink-like parametrization of the deceleration parameter ($q$). This $q$-parametrization can be written in terms of the initial ($q_i$) and present ($q_0$) values of the deceleration parameter, the redshift of the cosmic transition from deceleration to acceleration ($z_t$) and the redshift width of such transition ($\tau$). By assuming a flat space geometry, $q_i=1/2$ and adopting a likelihood approach to deal with the SN Ia data we obtain, at the 68% confidence level (C.L.), that: $z_t=0.56^{+0.13}_{-0.10}$, $\tau=0.47^{+0.16}_{-0.20}$ and $q_0=-0.31^{+0.11}_{-0.11}$ when we combine BAO/CMB observations with SN Ia data processed with the MLCS2k2 light-curve fitter. When in this combination we use the SALT2 fitter we get instead, at the same C.L.: $z_t=0.64^{+0.13}_{-0.07}$, $\tau=0.36^{+0.11}_{-0.17}$ and $q_0=-0.53^{+0.17}_{-0.13}$. Our results indicate, with a quite general and model independent approach, that MLCS2k2 favors Dvali-Gabadadze-Porrati-like cosmological models, while SALT2 favors $\Lambda$CDM-like ones. Progress in determining the transition redshift and/or the present value of the deceleration parameter depends crucially on solving the issue of the difference obtained when using these two light-curve fitters.

By using the Sloan Digital Sky Survey (SDSS) first year type Ia supernova (SN Ia) compilation, we compare two different approaches (traditional \chi^2 and complete likelihood) to determine parameter constraints when the magnitude dispersion is to be estimated as well. We consider cosmological constant + Cold Dark Matter (\Lambda CDM) and spatially flat, constant w Dark Energy + Cold Dark Matter (FwCDM) cosmological models and show that, for current data, there is a small difference in the best fit values and $\sim$ 30% difference in confidence contour areas in case the MLCS2k2 light-curve fitter is adopted. For the SALT2 light-curve fitter the differences are less significant ($\lesssim$ 13% difference in areas). In both cases the likelihood approach gives more restrictive constraints. We argue for the importance of using the complete likelihood instead of the \chi^2 approach when dealing with parameters in the expression for the variance.

The type Ia supernovae observational data is one of the most important in observational cosmology nowadays. Here we present the first public version of BETOCS (BayEsian Tools for Observational Cosmology using SNe Ia), which is a powerful and high productivity tool aimed to help the theoretical physicist community investigate cosmological models using type Ia supernovae (SNe Ia) observational data. BETOCS is applied to the generalized Chaplygin gas model (GCGM), traditional Chaplygin gas model (CGM) and $\Lambda$CDM, ranging from 5 to 3 free parameters, respectively. The ``gold sample'' of 157 supernovae data is used. It is shown that the Chaplygin gas scenario is viable (in most cases the $\Lambda$CDM is disfavoured) and the quartessence scenario (that unifies the description for dark matter and dark energy) is favoured. The Hubble parameter ($H_0$) is important and should not be fixed and it can be estimated or marginalized with or without the Hubble Space Telescope prior.