Group/Individual Details
- Shishen Chen (s3442140)
- Shipren Jayadev (s3744421)
Executive Statement
Is there a statistical significant difference in prices between Coles and Woolworths?
That is the aim of this investigation. To begin we have chosen to randomly sample products from both supermarkets primarily focusing on the meats and vegetables section. The reason we chose this product range is that they are easily comparable in terms of price and will have the least influence by branding other than the branding from Coles and Woolworths.
We have elected to choose a sample size of 40 products from each supermarket. The method we obtained the prices was through both supermarkets online store obtaining their prices per unit kilogram.
After conducting the hyphothesis test below, we can conclude that there is no statistical difference in the prices between Woolworths and Coles given the variables and samples taken for this investigation.
Load Packages and Data
library(tidyr)
library(dplyr)
library(magrittr)
library(car)
Summary Statistics
Use R to summarise the data from the investigation. Include an appropriate plot to help visualise the data. Describe the trend.
#Summary Stats
price <- read.csv("Price.csv")
price_g <- price %>% gather('Coles','Woolworths',key = "supermarket",value = "Price")
price_g <- price_g %>% mutate(supermarket = as.factor(supermarket))
price_g %>% group_by(supermarket) %>% summarise(Min = min(Price,na.rm = TRUE),
Q1 = quantile(Price,probs = 0.25,
na.rm = TRUE),
Median = median(Price,na.rm = TRUE),
Q3 = quantile(Price,probs = 0.75,
na.rm = TRUE),
Max = max(Price,na.rm = TRUE),
Mean = mean(Price, na.rm = TRUE),
SD = sd(Price,na.rm = TRUE),
n = n(),
Missing = sum(is.na(Price)))
price_g %>% boxplot(Price~supermarket, data = ., ylab = "Price($)")

Hypothesis Test
Use R to perform an appropriate hypothesis test to determine which supermarket is the cheapest. You need to explain your choice of hypothesis test, any assumptions and the significance level.
1.Comparing two independent sample means with unknown population variance.
2.Large sample used (n>30 for both groups) so normality can be assumed.
3.Population homogeneity of variance.
Mc = Mean of Coles Price
Mw = Mean of Woolworths Price
Variance1^2 = Variance of Coles Price
Variance2^2 = Variance of Woolworths Price
H0: Mc = Mw (No difference in prices between Supermarkets)
HA: Mc != Mw (Prices between Supermarkets are different)
HA: Mc > Mw (Coles is more expensive than Woolworths)
HA: Mc < Mw (Woolworths is more expensive than Coles)
# H0 : Variance1^2 = Variance2^2
# HA : Variance1^2 != Variance2^2
leveneTest(Price~supermarket, data = price_g)
Levene's Test for Homogeneity of Variance (center = median)
Df F value Pr(>F)
group 1 0.0129 0.91
78
#Two Sample t-test
t.test(Price~supermarket, data = price_g, var.equal = TRUE, alternative = "two.sided")
Two Sample t-test
data: Price by supermarket
t = 0.13941, df = 78, p-value = 0.8895
alternative hypothesis: true difference in means is not equal to 0
95 percent confidence interval:
-3.130926 3.602426
sample estimates:
mean in group Coles mean in group Woolworths
8.81150 8.57575
Interpretation
After doing the Levene Test, since p > 0.05, we can assume the variances are homogeneous. Samples can also be assumed to be approximately normal since they are both more than 30 in size.
The estimated diffrence in means between the two supermarkets is Mc - Mw = 8.81150 - 8.57575 = 0.23575 The 95% CI is [-3.130926, 3.602426]
p-value = 0.8895.
Since the p-value is more than 0.05, we will not be rejecting the null hyphothesis.
There is no statistical difference in the prices between Coles and Woolworths.
Discussion
Our findings confirmed our preliminary assumption wihch was there is no significant difference in the two supermarkets. The strength of our investigation comes from the product range we sampled from as meats and vegetables tend to be less influenced by brand names and tend to have a similar supplier cost so it limits price influences other than the ones imposed on the products by Coles and Woolworths.
One weakness would be that the sample size is small and a better approach in the future would be to select even larger sample sizes or to even consider stratified sampling.
LS0tDQp0aXRsZTogIk1BVEgxMzI0IEFzc2lnbm1lbnQgMiINCnN1YnRpdGxlOiAiU3VwZXJtYXJrZXQgUHJpY2UgV2FycyINCm91dHB1dDogaHRtbF9ub3RlYm9vaw0KLS0tDQoNCiMjIEdyb3VwL0luZGl2aWR1YWwgRGV0YWlscw0KDQoqIFNoaXNoZW4gQ2hlbiAoczM0NDIxNDApDQoqIFNoaXByZW4gSmF5YWRldiAoczM3NDQ0MjEpDQoNCiMjIEV4ZWN1dGl2ZSBTdGF0ZW1lbnQNCg0KSXMgdGhlcmUgYSBzdGF0aXN0aWNhbCBzaWduaWZpY2FudCBkaWZmZXJlbmNlIGluIHByaWNlcyBiZXR3ZWVuIENvbGVzIGFuZCBXb29sd29ydGhzPw0KDQpUaGF0IGlzIHRoZSBhaW0gb2YgdGhpcyBpbnZlc3RpZ2F0aW9uLiBUbyBiZWdpbiB3ZSBoYXZlIGNob3NlbiB0byByYW5kb21seSBzYW1wbGUgcHJvZHVjdHMgZnJvbSBib3RoIHN1cGVybWFya2V0cyBwcmltYXJpbHkgZm9jdXNpbmcgb24gdGhlIG1lYXRzIGFuZCB2ZWdldGFibGVzIHNlY3Rpb24uIFRoZSByZWFzb24gd2UgY2hvc2UgdGhpcyBwcm9kdWN0IHJhbmdlIGlzIHRoYXQgdGhleSBhcmUgZWFzaWx5IGNvbXBhcmFibGUgaW4gdGVybXMgb2YgcHJpY2UgYW5kIHdpbGwgaGF2ZSB0aGUgbGVhc3QgaW5mbHVlbmNlIGJ5IGJyYW5kaW5nIG90aGVyIHRoYW4gdGhlIGJyYW5kaW5nIGZyb20gQ29sZXMgYW5kIFdvb2x3b3J0aHMuDQoNCldlIGhhdmUgZWxlY3RlZCB0byBjaG9vc2UgYSBzYW1wbGUgc2l6ZSBvZiA0MCBwcm9kdWN0cyBmcm9tIGVhY2ggc3VwZXJtYXJrZXQuIFRoZSBtZXRob2Qgd2Ugb2J0YWluZWQgdGhlIHByaWNlcyB3YXMgdGhyb3VnaCBib3RoIHN1cGVybWFya2V0cyBvbmxpbmUgc3RvcmUgb2J0YWluaW5nIHRoZWlyIHByaWNlcyBwZXIgdW5pdCBraWxvZ3JhbS4NCg0KQWZ0ZXIgY29uZHVjdGluZyB0aGUgaHlwaG90aGVzaXMgdGVzdCBiZWxvdywgd2UgY2FuIGNvbmNsdWRlIHRoYXQgdGhlcmUgaXMgbm8gc3RhdGlzdGljYWwgZGlmZmVyZW5jZSBpbiB0aGUgcHJpY2VzIGJldHdlZW4gV29vbHdvcnRocyBhbmQgQ29sZXMgZ2l2ZW4gdGhlIHZhcmlhYmxlcyBhbmQgc2FtcGxlcyB0YWtlbiBmb3IgdGhpcyBpbnZlc3RpZ2F0aW9uLiANCg0KIyMgTG9hZCBQYWNrYWdlcyBhbmQgRGF0YQ0KDQpgYGB7cn0NCmxpYnJhcnkodGlkeXIpDQpsaWJyYXJ5KGRwbHlyKQ0KbGlicmFyeShtYWdyaXR0cikNCmxpYnJhcnkoY2FyKQ0KYGBgDQoNCiMjIFN1bW1hcnkgU3RhdGlzdGljcw0KDQpVc2UgUiB0byBzdW1tYXJpc2UgdGhlIGRhdGEgZnJvbSB0aGUgaW52ZXN0aWdhdGlvbi4gSW5jbHVkZSBhbiBhcHByb3ByaWF0ZSBwbG90IHRvIGhlbHAgdmlzdWFsaXNlIHRoZSBkYXRhLiBEZXNjcmliZSB0aGUgdHJlbmQuIA0KDQpgYGB7cn0NCiNTdW1tYXJ5IFN0YXRzDQoNCnByaWNlIDwtIHJlYWQuY3N2KCJQcmljZS5jc3YiKQ0KcHJpY2VfZyA8LSBwcmljZSAlPiUgZ2F0aGVyKCdDb2xlcycsJ1dvb2x3b3J0aHMnLGtleSA9ICJzdXBlcm1hcmtldCIsdmFsdWUgPSAiUHJpY2UiKQ0KcHJpY2VfZyA8LSBwcmljZV9nICU+JSBtdXRhdGUoc3VwZXJtYXJrZXQgPSBhcy5mYWN0b3Ioc3VwZXJtYXJrZXQpKQ0KcHJpY2VfZyAlPiUgZ3JvdXBfYnkoc3VwZXJtYXJrZXQpICU+JSBzdW1tYXJpc2UoTWluID0gbWluKFByaWNlLG5hLnJtID0gVFJVRSksDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBRMSA9IHF1YW50aWxlKFByaWNlLHByb2JzID0gMC4yNSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbmEucm0gPSBUUlVFKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIE1lZGlhbiA9IG1lZGlhbihQcmljZSxuYS5ybSA9IFRSVUUpLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgUTMgPSBxdWFudGlsZShQcmljZSxwcm9icyA9IDAuNzUsDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIG5hLnJtID0gVFJVRSksDQogICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICBNYXggPSBtYXgoUHJpY2UsbmEucm0gPSBUUlVFKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIE1lYW4gPSBtZWFuKFByaWNlLCBuYS5ybSA9IFRSVUUpLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgU0QgPSBzZChQcmljZSxuYS5ybSA9IFRSVUUpLA0KICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgbiA9IG4oKSwNCiAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgICAgIE1pc3NpbmcgPSBzdW0oaXMubmEoUHJpY2UpKSkNCg0KcHJpY2VfZyAlPiUgYm94cGxvdChQcmljZX5zdXBlcm1hcmtldCwgZGF0YSA9IC4sIHlsYWIgPSAiUHJpY2UoJCkiKQ0KYGBgDQoNCg0KIyMgSHlwb3RoZXNpcyBUZXN0DQoNClVzZSBSIHRvIHBlcmZvcm0gYW4gYXBwcm9wcmlhdGUgaHlwb3RoZXNpcyB0ZXN0IHRvIGRldGVybWluZSB3aGljaCBzdXBlcm1hcmtldCBpcyB0aGUgY2hlYXBlc3QuIFlvdSBuZWVkIHRvIGV4cGxhaW4geW91ciBjaG9pY2Ugb2YgaHlwb3RoZXNpcyB0ZXN0LCBhbnkgYXNzdW1wdGlvbnMgYW5kIHRoZSBzaWduaWZpY2FuY2UgbGV2ZWwuIA0KDQoxLkNvbXBhcmluZyB0d28gaW5kZXBlbmRlbnQgc2FtcGxlIG1lYW5zIHdpdGggdW5rbm93biBwb3B1bGF0aW9uIHZhcmlhbmNlLg0KDQoyLkxhcmdlIHNhbXBsZSB1c2VkIChuPjMwIGZvciBib3RoIGdyb3Vwcykgc28gbm9ybWFsaXR5IGNhbiBiZSBhc3N1bWVkLiANCg0KMy5Qb3B1bGF0aW9uIGhvbW9nZW5laXR5IG9mIHZhcmlhbmNlLg0KDQpNYyA9IE1lYW4gb2YgQ29sZXMgUHJpY2UNCg0KTXcgPSBNZWFuIG9mIFdvb2x3b3J0aHMgUHJpY2UNCg0KVmFyaWFuY2UxXjIgPSBWYXJpYW5jZSBvZiBDb2xlcyBQcmljZQ0KDQpWYXJpYW5jZTJeMiA9IFZhcmlhbmNlIG9mIFdvb2x3b3J0aHMgUHJpY2UNCg0KSDA6IE1jID0gTXcgKE5vIGRpZmZlcmVuY2UgaW4gcHJpY2VzIGJldHdlZW4gU3VwZXJtYXJrZXRzKQ0KDQpIQTogTWMgIT0gTXcgKFByaWNlcyBiZXR3ZWVuIFN1cGVybWFya2V0cyBhcmUgZGlmZmVyZW50KQ0KDQpIQTogTWMgPiBNdyAoQ29sZXMgaXMgbW9yZSBleHBlbnNpdmUgdGhhbiBXb29sd29ydGhzKQ0KDQpIQTogTWMgPCBNdyAoV29vbHdvcnRocyBpcyBtb3JlIGV4cGVuc2l2ZSB0aGFuIENvbGVzKQ0KDQoNCmBgYHtyfQ0KIyBIMCA6IFZhcmlhbmNlMV4yID0gVmFyaWFuY2UyXjINCiMgSEEgOiBWYXJpYW5jZTFeMiAhPSBWYXJpYW5jZTJeMg0KbGV2ZW5lVGVzdChQcmljZX5zdXBlcm1hcmtldCwgZGF0YSA9IHByaWNlX2cpDQoNCiNUd28gU2FtcGxlIHQtdGVzdA0KdC50ZXN0KFByaWNlfnN1cGVybWFya2V0LCBkYXRhID0gcHJpY2VfZywgdmFyLmVxdWFsID0gVFJVRSwgYWx0ZXJuYXRpdmUgPSAidHdvLnNpZGVkIikNCmBgYA0KDQojIyBJbnRlcnByZXRhdGlvbg0KDQpBZnRlciBkb2luZyB0aGUgTGV2ZW5lIFRlc3QsIHNpbmNlIHAgPiAwLjA1LCB3ZSBjYW4gYXNzdW1lIHRoZSB2YXJpYW5jZXMgYXJlIGhvbW9nZW5lb3VzLiANClNhbXBsZXMgY2FuIGFsc28gYmUgYXNzdW1lZCB0byBiZSBhcHByb3hpbWF0ZWx5IG5vcm1hbCBzaW5jZSB0aGV5IGFyZSBib3RoIG1vcmUgdGhhbiAzMCBpbiBzaXplLiANCg0KVGhlIGVzdGltYXRlZCBkaWZmcmVuY2UgaW4gbWVhbnMgYmV0d2VlbiB0aGUgdHdvIHN1cGVybWFya2V0cyBpcyBNYyAtIE13ID0gOC44MTE1MCAtIDguNTc1NzUgPSAwLjIzNTc1DQpUaGUgOTUlIENJIGlzIFstMy4xMzA5MjYsIDMuNjAyNDI2XQ0KDQpwLXZhbHVlID0gMC44ODk1LiANCg0KU2luY2UgdGhlIHAtdmFsdWUgaXMgbW9yZSB0aGFuIDAuMDUsIHdlIHdpbGwgbm90IGJlIHJlamVjdGluZyB0aGUgbnVsbCBoeXBob3RoZXNpcy4gDQoNClRoZXJlIGlzIG5vIHN0YXRpc3RpY2FsIGRpZmZlcmVuY2UgaW4gdGhlIHByaWNlcyBiZXR3ZWVuIENvbGVzIGFuZCBXb29sd29ydGhzLiANCg0KDQojIyBEaXNjdXNzaW9uIA0KDQpPdXIgZmluZGluZ3MgY29uZmlybWVkIG91ciBwcmVsaW1pbmFyeSBhc3N1bXB0aW9uIHdpaGNoIHdhcyB0aGVyZSBpcyBubyBzaWduaWZpY2FudCBkaWZmZXJlbmNlIGluIHRoZSB0d28gc3VwZXJtYXJrZXRzLiBUaGUgc3RyZW5ndGggb2Ygb3VyIGludmVzdGlnYXRpb24gY29tZXMgZnJvbSB0aGUgcHJvZHVjdCByYW5nZSB3ZSBzYW1wbGVkIGZyb20gYXMgbWVhdHMgYW5kIHZlZ2V0YWJsZXMgdGVuZCB0byBiZSBsZXNzIGluZmx1ZW5jZWQgYnkgYnJhbmQgbmFtZXMgYW5kIHRlbmQgdG8gaGF2ZSBhIHNpbWlsYXIgc3VwcGxpZXIgY29zdCBzbyBpdCBsaW1pdHMgcHJpY2UgaW5mbHVlbmNlcyBvdGhlciB0aGFuIHRoZSBvbmVzIGltcG9zZWQgb24gdGhlIHByb2R1Y3RzIGJ5IENvbGVzIGFuZCBXb29sd29ydGhzLiANCg0KT25lIHdlYWtuZXNzIHdvdWxkIGJlIHRoYXQgdGhlIHNhbXBsZSBzaXplIGlzIHNtYWxsIGFuZCBhIGJldHRlciBhcHByb2FjaCBpbiB0aGUgZnV0dXJlIHdvdWxkIGJlIHRvIHNlbGVjdCBldmVuIGxhcmdlciBzYW1wbGUgc2l6ZXMgb3IgdG8gZXZlbiBjb25zaWRlciBzdHJhdGlmaWVkIHNhbXBsaW5nLiANCg0K