library(readxl)
## Warning: package 'readxl' was built under R version 4.0.5
RestaurantGrades <-read_excel('RestaurantGrades.xlsx', sheet = "ads_analysis%202_1fd2e6d3-bf38-")
names(RestaurantGrades)
## [1] "treatment" "pageviews" "calls" "reservations"
## [5] "business_id" "restaurant_type"
RestaurantGrades <- RestaurantGrades[-5]
names(RestaurantGrades)
## [1] "treatment" "pageviews" "calls" "reservations"
## [5] "restaurant_type"
RestaurantGrades$treatment <- as.factor (RestaurantGrades$treatment)
levels (RestaurantGrades$treatment)
## [1] "0" "1" "2"
RestaurantGrades$restaurant_type <- as.factor (RestaurantGrades$restaurant_type)
levels (RestaurantGrades$restaurant_type)
## [1] "chain" "independent"
model.pageviews <- lm(pageviews~treatment, data = RestaurantGrades)
model.calls <- lm(calls~treatment, data = RestaurantGrades)
model.reservations <- lm(reservations~treatment, data = RestaurantGrades)
summary(model.pageviews)
##
## Call:
## lm(formula = pageviews ~ treatment, data = RestaurantGrades)
##
## Residuals:
## Min 1Q Median 3Q Max
## -295.21 -137.19 -86.78 172.81 434.79
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 419.779 1.645 255.18 <2e-16 ***
## treatment1 81.411 2.326 34.99 <2e-16 ***
## treatment2 63.432 2.326 27.27 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 164.5 on 29997 degrees of freedom
## Multiple R-squared: 0.04312, Adjusted R-squared: 0.04306
## F-statistic: 675.9 on 2 and 29997 DF, p-value: < 2.2e-16
summary(model.calls)
##
## Call:
## lm(formula = calls ~ treatment, data = RestaurantGrades)
##
## Residuals:
## Min 1Q Median 3Q Max
## -22.714 -5.388 -1.388 4.980 35.286
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 34.01960 0.07316 464.98 <2e-16 ***
## treatment1 3.36890 0.10347 32.56 <2e-16 ***
## treatment2 7.69490 0.10347 74.37 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 7.316 on 29997 degrees of freedom
## Multiple R-squared: 0.1564, Adjusted R-squared: 0.1563
## F-statistic: 2780 on 2 and 29997 DF, p-value: < 2.2e-16
summary(model.reservations)
##
## Call:
## lm(formula = reservations ~ treatment, data = RestaurantGrades)
##
## Residuals:
## Min 1Q Median 3Q Max
## -23.680 -5.021 -1.021 4.320 37.320
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 33.96040 0.07123 476.762 <2e-16 ***
## treatment1 0.06080 0.10074 0.604 0.546
## treatment2 7.72010 0.10074 76.637 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 7.123 on 29997 degrees of freedom
## Multiple R-squared: 0.2057, Adjusted R-squared: 0.2057
## F-statistic: 3885 on 2 and 29997 DF, p-value: < 2.2e-16
model.pageviews_restaurant_type <- lm(pageviews~treatment+restaurant_type, data = RestaurantGrades)
model.calls_restaurant_type <- lm(calls~treatment+restaurant_type, data = RestaurantGrades)
model.reservations_restaurant_type <- lm(reservations~treatment+restaurant_type, data = RestaurantGrades)
summary(model.pageviews_restaurant_type)
##
## Call:
## lm(formula = pageviews ~ treatment + restaurant_type, data = RestaurantGrades)
##
## Residuals:
## Min 1Q Median 3Q Max
## -249.31 -32.31 -0.70 32.30 242.67
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 611.8987 0.6065 1008.98 <2e-16 ***
## treatment1 81.4114 0.7003 116.26 <2e-16 ***
## treatment2 63.4316 0.7003 90.58 <2e-16 ***
## restaurant_typeindependent -320.1988 0.5836 -548.70 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 49.52 on 29996 degrees of freedom
## Multiple R-squared: 0.9133, Adjusted R-squared: 0.9133
## F-statistic: 1.053e+05 on 3 and 29996 DF, p-value: < 2.2e-16
summary(model.calls_restaurant_type)
##
## Call:
## lm(formula = calls ~ treatment + restaurant_type, data = RestaurantGrades)
##
## Residuals:
## Min 1Q Median 3Q Max
## -29.0634 -3.1559 -0.0634 3.2130 28.9366
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 40.36848 0.06323 638.43 <2e-16 ***
## treatment1 3.36890 0.07301 46.14 <2e-16 ***
## treatment2 7.69490 0.07301 105.39 <2e-16 ***
## restaurant_typeindependent -10.58147 0.06084 -173.91 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 5.163 on 29996 degrees of freedom
## Multiple R-squared: 0.5799, Adjusted R-squared: 0.5799
## F-statistic: 1.38e+04 on 3 and 29996 DF, p-value: < 2.2e-16
summary(model.reservations_restaurant_type)
##
## Call:
## lm(formula = reservations ~ treatment + restaurant_type, data = RestaurantGrades)
##
## Residuals:
## Min 1Q Median 3Q Max
## -29.8025 -2.9399 0.0601 3.0601 31.1975
##
## Coefficients:
## Estimate Std. Error t value Pr(>|t|)
## (Intercept) 40.08237 0.06215 644.918 <2e-16 ***
## treatment1 0.06080 0.07177 0.847 0.397
## treatment2 7.72010 0.07177 107.573 <2e-16 ***
## restaurant_typeindependent -10.20328 0.05980 -170.609 <2e-16 ***
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
##
## Residual standard error: 5.075 on 29996 degrees of freedom
## Multiple R-squared: 0.5969, Adjusted R-squared: 0.5969
## F-statistic: 1.481e+04 on 3 and 29996 DF, p-value: < 2.2e-16
library(ggplot2)
## Warning: package 'ggplot2' was built under R version 4.0.5
library(dplyr)
## Warning: package 'dplyr' was built under R version 4.0.5
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
RestaurantGrades %>%
ggplot(aes(x=restaurant_type, y=pageviews, col = restaurant_type))+
geom_boxplot(alpha = 1)+
geom_line(method = lm)+
facet_wrap(~treatment)
## Warning: Ignoring unknown parameters: method
RestaurantGrades %>%
ggplot(aes(x=restaurant_type, y=calls, col = restaurant_type))+
geom_boxplot()+
geom_line(method = lm)+
facet_wrap(~treatment)
## Warning: Ignoring unknown parameters: method
RestaurantGrades %>%
ggplot(aes(x=restaurant_type, y=reservations, col = restaurant_type))+
geom_boxplot(alpha = 1)+
geom_line(method = lm)+
facet_wrap(~treatment)
## Warning: Ignoring unknown parameters: method
