#(a)
> data = read.csv("https://c...content-available-to-author-only...n.de/book/pizza_delivery.csv")
> data
> a
<- lm
(time ~temperature
+ branch
+ day
+ operator
+ driver
+ bill
+ pizzas
+ discount_customer
,data
= data
) > summary(a)
> #(b)
> lcl
<-coefficients
(a
) - qt
(0.975,1248) * sqrt(diag
(vcov
(a
))) > ucl
<-coefficients
(a
) + qt
(0.975,1248) * sqrt(diag
(vcov
(a
))) > cbind(coefficients(a),lcl,ucl)
> #(c)
> sum(residuals(a)^2)/(a$df.residual)
> #(d)
> t = data[,3]
> t
> st<-sum((t - mean(t))^2)
> se = sum(residuals(a)^2)
> rsuare = 1-(se/st)
> radj = 1-((se/1248)/(st/1265))
> rsuare
> radj
> #(e)
> a2
<- lm
(time~temperature
+ branch
+ day
+ driver
+ bill
+ pizzas
,data
= data
) > predict(a2,data[031,])
> #end
> # your code goes here# your code goes here
IyhhKQo+IGRhdGEgPSByZWFkLmNzdigiaHR0cHM6Ly9jLi4uY29udGVudC1hdmFpbGFibGUtdG8tYXV0aG9yLW9ubHkuLi5uLmRlL2Jvb2svcGl6emFfZGVsaXZlcnkuY3N2IikKPiBkYXRhCj4gYSA8LSBsbSh0aW1lIH50ZW1wZXJhdHVyZSArIGJyYW5jaCArIGRheSArIG9wZXJhdG9yICsgZHJpdmVyICsgYmlsbCArIHBpenphcyArIGRpc2NvdW50X2N1c3RvbWVyLGRhdGEgPSBkYXRhKQo+IHN1bW1hcnkoYSkKCj4gIyhiKQo+IGxjbDwtY29lZmZpY2llbnRzKGEpIC0gcXQoMC45NzUsMTI0OCkgKiBzcXJ0KGRpYWcodmNvdihhKSkpCj4gdWNsPC1jb2VmZmljaWVudHMoYSkgKyBxdCgwLjk3NSwxMjQ4KSAqIHNxcnQoZGlhZyh2Y292KGEpKSkKPiBjYmluZChjb2VmZmljaWVudHMoYSksbGNsLHVjbCkKPiAjKGMpCj4gc3VtKHJlc2lkdWFscyhhKV4yKS8oYSRkZi5yZXNpZHVhbCkKPiAjKGQpCj4gdCA9IGRhdGFbLDNdCj4gdAo+IHN0PC1zdW0oKHQgLSBtZWFuKHQpKV4yKQo+IHNlID0gc3VtKHJlc2lkdWFscyhhKV4yKQo+IHJzdWFyZSA9IDEtKHNlL3N0KQo+IHJhZGogPSAxLSgoc2UvMTI0OCkvKHN0LzEyNjUpKQo+IHJzdWFyZQo+IHJhZGoKPiAjKGUpCj4gYTI8LSBsbSh0aW1lfnRlbXBlcmF0dXJlICsgYnJhbmNoICsgZGF5ICsgZHJpdmVyICsgYmlsbCArIHBpenphcyxkYXRhID0gZGF0YSkKPiBwcmVkaWN0KGEyLGRhdGFbMDMxLF0pCgo+ICNlbmQKPiAjIHlvdXIgY29kZSBnb2VzIGhlcmUjIHlvdXIgY29kZSBnb2VzIGhlcmU=