lecture_11_live

R Markdown

This is an R Markdown document. Markdown is a simple formatting syntax for authoring HTML, PDF, and MS Word documents. For more details on using R Markdown see http://rmarkdown.rstudio.com.

When you click the Knit button a document will be generated that includes both content as well as the output of any embedded R code chunks within the document. You can embed an R code chunk like this:

\(x^2 + 5x\)

\[ \operatorname{Var}(X) = E[(X-\mu)^2] \] A list:

1. item 1
2. item 2
3. item 3

plot with base R:

plot(pressure, type = 'l', pch = 19, col = 'blue',
     main = 'Vapor pressure',
     xlab = 'Temp',
     ylab = 'Pressure')

Including Plots

You can also embed plots, for example:

Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.

Including tables (knitr::kable)

Head of first 6 mtcars

	mpg	cyl	disp	hp	drat	wt	qsec	vs	am	gear	carb
Mazda RX4	21.0	6	160	110	3.90	2.620	16.46	0	1	4	4
Mazda RX4 Wag	21.0	6	160	110	3.90	2.875	17.02	0	1	4	4
Datsun 710	22.8	4	108	93	3.85	2.320	18.61	1	1	4	1
Hornet 4 Drive	21.4	6	258	110	3.08	3.215	19.44	1	0	3	1
Hornet Sportabout	18.7	8	360	175	3.15	3.440	17.02	0	0	3	2
Valiant	18.1	6	225	105	2.76	3.460	20.22	1	0	3	1

Figures

x <- rnorm(1000, mean = 0, sd = 1)
hist(x, breaks = 30, col = 'grey', border = 'white')
rug(x, col = 'red')

Histogram with Base R

knitr::include_graphics('lecture11_sample.jpg')

Figure 3. Included sample image.

Lecture 11 sample

plot(cars$speed, cars$dist,
     pch = 19, col = 'tomato')
abline(lm(dist ~ speed, data = cars), col = 'grey', lwd = 2, lty = 2)

my_model <- lm(dist ~ speed, data = cars)
summary(my_model)

## 
## Call:
## lm(formula = dist ~ speed, data = cars)
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -29.069  -9.525  -2.272   9.215  43.201 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) -17.5791     6.7584  -2.601   0.0123 *  
## speed         3.9324     0.4155   9.464 1.49e-12 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 15.38 on 48 degrees of freedom
## Multiple R-squared:  0.6511, Adjusted R-squared:  0.6438 
## F-statistic: 89.57 on 1 and 48 DF,  p-value: 1.49e-12