Session 14

adv_data=read.csv("AdvertisingDataV2.csv")
library(psych)
describe(adv_data)

                vars     n     mean      sd  median  trimmed      mad min
adType*            1 30000     2.00    0.82     2.0     2.00     1.48   1
pageViews          2 30000   468.06  168.16   391.0   458.18   149.74 145
phoneCalls         3 30000    37.71    7.97    37.0    37.20     7.41  17
reservations       4 30000    36.55    7.99    36.0    35.97     7.41  15
businessID         5 30000 15000.50 8660.40 15000.5 15000.50 11119.50   1
restaurantType*    6 30000     1.60    0.49     2.0     1.62     0.00   1
                  max range  skew kurtosis    se
adType*             3     2  0.00    -1.50  0.00
pageViews         929   784  0.45    -1.29  0.97
phoneCalls         77    60  0.65     0.57  0.05
reservations       79    64  0.78     0.88  0.05
businessID      30000 29999  0.00    -1.20 50.00
restaurantType*     2     1 -0.41    -1.83  0.00

attach(adv_data) 

oneWayfit <- aov(reservations ~ adType, data = adv_data)
# summary of the ANOVA model
summary(oneWayfit)

               Df  Sum Sq Mean Sq F value Pr(>F)    
adType          2  394228  197114    3885 <2e-16 ***
Residuals   29997 1522018      51                   
---
Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

Error in unique.default(x, nmax = nmax) : 
  unique() applies only to vectors