0.1 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 .

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:

summary (cars)
##      speed           dist       
##  Min.   : 4.0   Min.   :  2.00  
##  1st Qu.:12.0   1st Qu.: 26.00  
##  Median :15.0   Median : 36.00  
##  Mean   :15.4   Mean   : 42.98  
##  3rd Qu.:19.0   3rd Qu.: 56.00  
##  Max.   :25.0   Max.   :120.00

0.2 Including Plots

You can also embed plots, for example:

plot (pressure)

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

# install.packages("tidyverse")
# install.packages("datasets")
library(tidyverse)
## Warning: package 'tidyverse' was built under R version 4.3.3
## Warning: package 'ggplot2' was built under R version 4.3.3
## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.4
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ ggplot2   3.5.0     ✔ tibble    3.2.1
## ✔ lubridate 1.9.3     ✔ tidyr     1.3.0
## ✔ purrr     1.0.2     
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors
library(datasets)

1 chickwts

view(chickwts)
# For the chicken
head(chickwts)
##   weight      feed
## 1    179 horsebean
## 2    160 horsebean
## 3    136 horsebean
## 4    227 horsebean
## 5    217 horsebean
## 6    168 horsebean
tail(chickwts)
##    weight   feed
## 66    352 casein
## 67    359 casein
## 68    216 casein
## 69    222 casein
## 70    283 casein
## 71    332 casein
glimpse(chickwts)
## Rows: 71
## Columns: 2
## $ weight <dbl> 179, 160, 136, 227, 217, 168, 108, 124, 143, 140, 309, 229, 181…
## $ feed   <fct> horsebean, horsebean, horsebean, horsebean, horsebean, horsebea…
str(chickwts)
## 'data.frame':    71 obs. of  2 variables:
##  $ weight: num  179 160 136 227 217 168 108 124 143 140 ...
##  $ feed  : Factor w/ 6 levels "casein","horsebean",..: 2 2 2 2 2 2 2 2 2 2 ...
# Get the mean and standard deviation for the chicken weight
mean(chickwts$weight)
## [1] 261.3099
median(chickwts$weight)
## [1] 258
sd(chickwts$weight)
## [1] 78.0737
?(chickwts)
boxplot(chickwts$weight~chickwts$feed)

table(chickwts$feed)
## 
##    casein horsebean   linseed  meatmeal   soybean sunflower 
##        12        10        12        11        14        12
pie(table(chickwts$feed))

`

LS0tDQp0aXRsZTogIlNlc2lvbiAiDQphdXRob3I6ICJBbG1hIFNhbnRpYWdvIChBMDA4MzY2MzYpIg0KZGF0ZTogIjI0LTA2LTIwMjQiDQpvdXRwdXQ6DQogIGh0bWxfZG9jdW1lbnQ6DQogICAgdG9jOiB0cnVlIA0KICAgIHRvY19mbG9hdDogdHJ1ZQ0KICAgIGNvZGVfZG93bmxvYWQ6IHRydWUNCiAgICB0aGVtZTogam91cm5hbA0KICAgIG51bWJlcl9zZWN0aW9uczogdHJ1ZQ0KZWRpdG9yX3M6IA0KICBjaHVua19vdXRwdXRfdHlwZTogY29uc29sZQ0KZWRpdG9yX29wdGlvbnM6IA0KICBjaHVua19vdXRwdXRfdHlwZTogY29uc29sZQ0KLS0tDQoNCiMjIFIgTWFya2Rvd24gDQoNClRoaXMgaXMgYW4gUiBNYXJrZG93biBkb2N1bWVudC4gTWFya2Rvd24gaXMgYSBzaW1wbGUgZm9ybWF0dGluZyBzeW50YXggZm9yIGF1dGhvcmluZyBIVE1MLCBQREYsIGFuZCBNUyBXb3JkIGRvY3VtZW50cy4gRm9yIG1vcmUgZGV0YWlscyBvbiB1c2luZyBSIE1hcmtkb3duIHNlZSAuDQoNCldoZW4geW91IGNsaWNrIHRoZSAqKktuaXQqKiBidXR0b24gYSBkb2N1bWVudCB3aWxsIGJlIGdlbmVyYXRlZCB0aGF0IGluY2x1ZGVzIGJvdGggY29udGVudCBhcyB3ZWxsIGFzIHRoZSBvdXRwdXQgb2YgYW55IGVtYmVkZGVkIFIgY29kZSBjaHVua3Mgd2l0aGluIHRoZSBkb2N1bWVudC4gWW91IGNhbiBlbWJlZCBhbiBSIGNvZGUgY2h1bmsgbGlrZSB0aGlzOg0KDQpgYGB7ciBjYXJzfQ0Kc3VtbWFyeSAoY2FycykNCmBgYA0KDQojIyBJbmNsdWRpbmcgUGxvdHMNCg0KWW91IGNhbiBhbHNvIGVtYmVkIHBsb3RzLCBmb3IgZXhhbXBsZTogDQoNCmBgYHtyfQ0KcGxvdCAocHJlc3N1cmUpDQpgYGANCg0KTm90ZSB0aGF0IHRoZSBgZWNobyA9IEZBTFNFYCBwYXJhbWV0ZXIgd2FzIGFkZGVkIHRvIHRoZSBjb2RlIGNodW5rIHRvIHByZXZlbnQgcHJpbnRpbmcgb2YgdGhlIFIgY29kZSB0aGF0IGdlbmVyYXRlZCB0aGUgcGxvdC4NCg0KYGBge3J9DQojIGluc3RhbGwucGFja2FnZXMoInRpZHl2ZXJzZSIpDQojIGluc3RhbGwucGFja2FnZXMoImRhdGFzZXRzIikNCg0KYGBgDQoNCmBgYHtyfQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpsaWJyYXJ5KGRhdGFzZXRzKQ0KYGBgDQoNCiMgPHNwYW4gc3R5bGU9ImNvbG9yOiM4QjNBNjI7Ij5jaGlja3d0czwvc3Bhbj4NCg0KYGBge3IgbWVzc2FnZT1UUlVFLCB3YXJuaW5nPUZBTFNFfQ0KdmlldyhjaGlja3d0cykNCiMgRm9yIHRoZSBjaGlja2VuDQpoZWFkKGNoaWNrd3RzKQ0KdGFpbChjaGlja3d0cykNCmdsaW1wc2UoY2hpY2t3dHMpDQpzdHIoY2hpY2t3dHMpDQojIEdldCB0aGUgbWVhbiBhbmQgc3RhbmRhcmQgZGV2aWF0aW9uIGZvciB0aGUgY2hpY2tlbiB3ZWlnaHQNCm1lYW4oY2hpY2t3dHMkd2VpZ2h0KQ0KbWVkaWFuKGNoaWNrd3RzJHdlaWdodCkNCnNkKGNoaWNrd3RzJHdlaWdodCkNCj8oY2hpY2t3dHMpDQpib3hwbG90KGNoaWNrd3RzJHdlaWdodH5jaGlja3d0cyRmZWVkKQ0KdGFibGUoY2hpY2t3dHMkZmVlZCkNCnBpZSh0YWJsZShjaGlja3d0cyRmZWVkKSkNCg0KDQoNCmBgYA0KDQpg