**Wiekvoet**, and kindly contributed to R-bloggers)

A while ago I had a post: ‘Should I use premium Diesel? Setup. Since that time the data has been acquired. This post describes the results.

### Data

Data is registered by me in 2014 and 2015. 2014 has standard Diesel, while 2015 has premium. Both are from the same months of the year, so should have approximately the same weather and traffic.

### Results

The figure shows the main result: While the 2015 has some better results, it is not by much. It is not in the ballpark of 5% which was originally defined.

Results from distribution fitting and t.test confirm the visual observation.

#### fitdistr

Standard

mean sd

3.59517971 0.19314598

(0.04828649) (0.03414371)

Premium

mean sd

3.54015573 0.23905047

(0.05976262) (0.04225855)

#### t.test

Welch Two Sample t-test

data: usage by kind

t = -0.6934, df = 28.732, p-value = 0.4936

alternative hypothesis: true difference in means is not equal to 0

95 percent confidence interval:

-0.2173822 0.1073343

sample estimates:

mean in group Premium mean in group Standard

3.540156 3.595180

#### Posterior density for premium Diesel

Previously I defined a prior for the analysis of the data. This was a mixed distribution from three parts. Like Standard fuel, 5% better and something else in the neighborhood of 3.6. Given the simplicity of the model and the non-standard prior I used MCMCpack’s MCMCmetrop1R function. This allowed me to write it all as R code.

The result is a density for mean usage quite close to 3.6.

Iterations = 501:20500

Thinning interval = 1

Number of chains = 1

Sample size per chain = 20000

1. Empirical mean and standard deviation for each variable,

plus standard error of the mean:

Mean SD Naive SE Time-series SE

[1,] 3.5409 0.06413 0.0004535 0.001722

[2,] 0.2699 0.05353 0.0003785 0.001448

2. Quantiles for each variable:

2.5% 25% 50% 75% 97.5%

var1 3.4147 3.494 3.5465 3.5884 3.6523

var2 0.1873 0.232 0.2629 0.2998 0.3953

### Conclusion

### Data and code

**leave a comment**for the author, please follow the link and comment on their blog:

**Wiekvoet**.

R-bloggers.com offers

**daily e-mail updates**about R news and tutorials on topics such as: Data science, Big Data, R jobs, visualization (ggplot2, Boxplots, maps, animation), programming (RStudio, Sweave, LaTeX, SQL, Eclipse, git, hadoop, Web Scraping) statistics (regression, PCA, time series, trading) and more...