Homework1

Quarto

Quarto enables you to weave together content and executable code into a finished document. To learn more about Quarto see https://quarto.org.

Running Code

When you click the Render button a document will be generated that includes both content and the output of embedded code. You can embed code like this:

library(tidyverse)

Warning: package 'tidyverse' was built under R version 4.1.3

Warning: package 'ggplot2' was built under R version 4.1.3

Warning: package 'tibble' was built under R version 4.1.3

Warning: package 'tidyr' was built under R version 4.1.3

Warning: package 'readr' was built under R version 4.1.3

Warning: package 'purrr' was built under R version 4.1.3

Warning: package 'dplyr' was built under R version 4.1.3

Warning: package 'stringr' was built under R version 4.1.3

Warning: package 'forcats' was built under R version 4.1.3

Warning: package 'lubridate' was built under R version 4.1.3

-- Attaching core tidyverse packages ------------------------ tidyverse 2.0.0 --
v dplyr     1.1.2     v readr     2.1.4
v forcats   1.0.0     v stringr   1.5.0
v ggplot2   3.4.2     v tibble    3.2.1
v lubridate 1.9.2     v tidyr     1.3.0
v purrr     1.0.1     
-- Conflicts ------------------------------------------ tidyverse_conflicts() --
x dplyr::filter() masks stats::filter()
x dplyr::lag()    masks stats::lag()
i Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

library(car)

Warning: package 'car' was built under R version 4.1.3

Loading required package: carData

Attaching package: 'car'

The following object is masked from 'package:dplyr':

    recode

The following object is masked from 'package:purrr':

    some

data1 <- read.csv("Milk_Production.csv", colClasses = c(NA, "numeric"))
head(data1)

  EweType Diet Week MilkProd
1       2    1    0     1260
2       2    1    0     1900
3       2    1    0     1380
4       2    1    0     3910
5       2    1    0     2900
6       2    2    0     2210

str(data1)

'data.frame':   240 obs. of  4 variables:
 $ EweType : int  2 2 2 2 2 2 2 2 2 2 ...
 $ Diet    : num  1 1 1 1 1 2 2 2 2 2 ...
 $ Week    : int  0 0 0 0 0 0 0 0 0 0 ...
 $ MilkProd: num  1260 1900 1380 3910 2900 2210 1850 1750 1900 1950 ...

glimpse(data1)

Rows: 240
Columns: 4
$ EweType  <int> 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 1, 1, 1, 1, 1, 1~
$ Diet     <dbl> 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1, 2~
$ Week     <int> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0~
$ MilkProd <dbl> 1260, 1900, 1380, 3910, 2900, 2210, 1850, 1750, 1900, 1950, 1~

data1$EweType<- as.factor(data1$EweType)
data1$Diet<- as.factor(data1$Diet)
data1$Week <- as.character(data1$Week)
data1$MilkProd<- as.numeric(data1$MilkProd)
tail(data1)

    EweType Diet Week MilkProd
235       1    2    7     1000
236       1    3    7      550
237       1    3    7      820
238       1    3    7      520
239       1    3    7      630
240       1    3    7      500

data1

    EweType Diet Week MilkProd
1         2    1    0     1260
2         2    1    0     1900
3         2    1    0     1380
4         2    1    0     3910
5         2    1    0     2900
6         2    2    0     2210
7         2    2    0     1850
8         2    2    0     1750
9         2    2    0     1900
10        2    2    0     1950
11        2    3    0     1190
12        2    3    0     1190
13        2    3    0      590
14        2    3    0     1410
15        2    3    0     1320
16        1    1    0     4120
17        1    1    0     3530
18        1    1    0     1940
19        1    1    0     3130
20        1    1    0     3230
21        1    2    0     1480
22        1    2    0     2180
23        1    2    0     2610
24        1    2    0     2320
25        1    2    0     1740
26        1    3    0     1740
27        1    3    0     1680
28        1    3    0     1890
29        1    3    0     1970
30        1    3    0     2770
31        2    1    1     1640
32        2    1    1     2570
33        2    1    1     2760
34        2    1    1     3190
35        2    1    1     2440
36        2    2    1     2330
37        2    2    1      930
38        2    2    1     1810
39        2    2    1     1500
40        2    2    1     1800
41        2    3    1      760
42        2    3    1     1030
43        2    3    1      910
44        2    3    1      990
45        2    3    1      950
46        1    1    1     4810
47        1    1    1     3040
48        1    1    1     2000
49        1    1    1     2440
50        1    1    1     4040
51        1    2    1     1870
52        1    2    1     1700
53        1    2    1     1850
54        1    2    1     1420
55        1    2    1     1260
56        1    3    1     1190
57        1    3    1     1240
58        1    3    1     1380
59        1    3    1     1170
60        1    3    1     1790
61        2    1    2     1080
62        2    1    2     1100
63        2    1    2     1550
64        2    1    2     1050
65        2    1    2     1280
66        2    2    2     1700
67        2    2    2      850
68        2    2    2     1540
69        2    2    2      950
70        2    2    2     1350
71        2    3    2      540
72        2    3    2      700
73        2    3    2      600
74        2    3    2      750
75        2    3    2      700
76        1    1    2     1850
77        1    1    2     1380
78        1    1    2     1650
79        1    1    2     1230
80        1    1    2     1580
81        1    2    2     1500
82        1    2    2     1380
83        1    2    2     1120
84        1    2    2     1130
85        1    2    2      860
86        1    3    2      960
87        1    3    2     1080
88        1    3    2     1120
89        1    3    2      960
90        1    3    2     1020
91        2    1    3     1000
92        2    1    3     1520
93        2    1    3      980
94        2    1    3     1180
95        2    1    3     1540
96        2    2    3     1180
97        2    2    3      780
98        2    2    3     1210
99        2    2    3      500
100       2    2    3     1000
101       2    3    3      500
102       2    3    3      500
103       2    3    3      320
104       2    3    3      240
105       2    3    3      460
106       1    1    3     1600
107       1    1    3     1250
108       1    1    3     1420
109       1    1    3     1020
110       1    1    3     1450
111       1    2    3     1300
112       1    2    3     1100
113       1    2    3      750
114       1    2    3      900
115       1    2    3      640
116       1    3    3      600
117       1    3    3      980
118       1    3    3      940
119       1    3    3      750
120       1    3    3      570
121       2    1    4      950
122       2    1    4     1200
123       2    1    4     1380
124       2    1    4      980
125       2    1    4      740
126       2    2    4      890
127       2    2    4      380
128       2    2    4      810
129       2    2    4      260
130       2    2    4      800
131       2    3    4      370
132       2    3    4      350
133       2    3    4      250
134       2    3    4      170
135       2    3    4      380
136       1    1    4     1680
137       1    1    4     1200
138       1    1    4      960
139       1    1    4     1100
140       1    1    4      620
141       1    2    4      900
142       1    2    4      720
143       1    2    4      680
144       1    2    4      480
145       1    2    4      800
146       1    3    4      700
147       1    3    4      800
148       1    3    4      680
149       1    3    4      400
150       1    3    4      420
151       2    1    5     1000
152       2    1    5     1300
153       2    1    5     1000
154       2    1    5      850
155       2    1    5      800
156       2    2    5      850
157       2    2    5      300
158       2    2    5      600
159       2    2    5      300
160       2    2    5      650
161       2    3    5      250
162       2    3    5      250
163       2    3    5      200
164       2    3    5      100
165       2    3    5      270
166       1    1    5     1380
167       1    1    5      920
168       1    1    5     1280
169       1    1    5     1120
170       1    1    5      800
171       1    2    5      640
172       1    2    5      800
173       1    2    5      420
174       1    2    5      380
175       1    2    5      900
176       1    3    5      660
177       1    3    5      760
178       1    3    5      640
179       1    3    5      480
180       1    3    5      460
181       2    1    6      900
182       2    1    6     1350
183       2    1    6     1300
184       2    1    6     1100
185       2    1    6      600
186       2    2    6      750
187       2    2    6      560
188       2    2    6      750
189       2    2    6      250
190       2    2    6      600
191       2    3    6      360
192       2    3    6      300
193       2    3    6      150
194       2    3    6      200
195       2    3    6      200
196       1    1    6     1100
197       1    1    6      850
198       1    1    6     1100
199       1    1    6     1200
200       1    1    6      750
201       1    2    6      400
202       1    2    6      520
203       1    2    6      500
204       1    2    6      400
205       1    2    6      950
206       1    3    6      650
207       1    3    6      900
208       1    3    6      850
209       1    3    6      500
210       1    3    6      400
211       2    1    7      800
212       2    1    7     1250
213       2    1    7     1050
214       2    1    7      950
215       2    1    7      500
216       2    2    7      600
217       2    2    7      350
218       2    2    7      500
219       2    2    7      200
220       2    2    7      500
221       2    3    7      300
222       2    3    7      100
223       2    3    7      100
224       2    3    7      200
225       2    3    7      100
226       1    1    7     1000
227       1    1    7     1150
228       1    1    7     1120
229       1    1    7     1450
230       1    1    7      660
231       1    2    7      780
232       1    2    7      620
233       1    2    7      550
234       1    2    7      450
235       1    2    7     1000
236       1    3    7      550
237       1    3    7      820
238       1    3    7      520
239       1    3    7      630
240       1    3    7      500

You can add options to executable code like this

library(tidyverse)
library(car)
data2 <- read.csv("Liveweight.csv", colClasses = c(NA, "numeric"))
head(data2)

  Ewetype Diet Week LiveWeigth
1       1    1    0       67.5
2       1    1    0       86.0
3       1    1    0       85.5
4       1    1    0       74.5
5       1    1    0       84.0
6       1    2    0       80.0

str(data2)

'data.frame':   240 obs. of  4 variables:
 $ Ewetype   : int  1 1 1 1 1 1 1 1 1 1 ...
 $ Diet      : num  1 1 1 1 1 2 2 2 2 2 ...
 $ Week      : int  0 0 0 0 0 0 0 0 0 0 ...
 $ LiveWeigth: num  67.5 86 85.5 74.5 84 80 76.5 71.5 75 78.5 ...

glimpse(data2)

Rows: 240
Columns: 4
$ Ewetype    <int> 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2,~
$ Diet       <dbl> 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 1, 1, 1, 1, 1,~
$ Week       <int> 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,~
$ LiveWeigth <dbl> 67.5, 86.0, 85.5, 74.5, 84.0, 80.0, 76.5, 71.5, 75.0, 78.5,~

data2$Ewetype<- as.factor(data2$Ewetype)
data2$Diet<- as.factor(data2$Diet)
data2$Week <- as.character(data2$Week)
data2$LiveWeigth<- as.numeric(data2$LiveWeigth)
tail(data2)

    Ewetype Diet Week LiveWeigth
235       2    2    7       63.0
236       2    3    7       57.5
237       2    3    7       62.5
238       2    3    7       62.5
239       2    3    7       83.5
240       2    3    7       80.0

data2

    Ewetype Diet Week LiveWeigth
1         1    1    0       67.5
2         1    1    0       86.0
3         1    1    0       85.5
4         1    1    0       74.5
5         1    1    0       84.0
6         1    2    0       80.0
7         1    2    0       76.5
8         1    2    0       71.5
9         1    2    0       75.0
10        1    2    0       78.5
11        1    3    0       70.0
12        1    3    0       93.5
13        1    3    0       77.5
14        1    3    0       91.5
15        1    3    0       83.5
16        2    1    0       79.5
17        2    1    0       88.5
18        2    1    0       76.5
19        2    1    0       59.5
20        2    1    0       90.5
21        2    2    0       87.5
22        2    2    0       70.0
23        2    2    0       72.0
24        2    2    0       64.0
25        2    2    0       62.5
26        2    3    0       55.5
27        2    3    0       65.0
28        2    3    0       64.0
29        2    3    0       81.0
30        2    3    0       81.5
31        1    1    1       68.0
32        1    1    1       88.0
33        1    1    1       87.0
34        1    1    1       77.5
35        1    1    1       85.5
36        1    2    1       81.0
37        1    2    1       71.0
38        1    2    1       72.0
39        1    2    1       75.5
40        1    2    1       81.0
41        1    3    1       66.5
42        1    3    1       89.0
43        1    3    1       74.0
44        1    3    1       86.5
45        1    3    1       75.5
46        2    1    1       82.5
47        2    1    1       89.5
48        2    1    1       77.0
49        2    1    1       62.5
50        2    1    1       91.0
51        2    2    1       85.5
52        2    2    1       69.5
53        2    2    1       70.5
54        2    2    1       63.0
55        2    2    1       62.0
56        2    3    1       54.5
57        2    3    1       64.0
58        2    3    1       63.5
59        2    3    1       78.5
60        2    3    1       81.0
61        1    1    2       68.0
62        1    1    2       90.0
63        1    1    2       88.0
64        1    1    2       80.5
65        1    1    2       87.5
66        1    2    2       81.0
67        1    2    2       67.5
68        1    2    2       72.0
69        1    2    2       75.5
70        1    2    2       82.5
71        1    3    2       62.0
72        1    3    2       85.5
73        1    3    2       70.0
74        1    3    2       78.5
75        1    3    2       69.5
76        2    1    2       86.0
77        2    1    2       88.0
78        2    1    2       77.0
79        2    1    2       64.0
80        2    1    2       91.0
81        2    2    2       84.5
82        2    2    2       68.0
83        2    2    2       67.0
84        2    2    2       62.0
85        2    2    2       60.5
86        2    3    2       55.0
87        2    3    2       62.5
88        2    3    2       64.0
89        2    3    2       77.0
90        2    3    2       80.0
91        1    1    3       66.5
92        1    1    3       89.5
93        1    1    3       87.0
94        1    1    3       84.0
95        1    1    3       89.5
96        1    2    3       80.0
97        1    2    3       67.5
98        1    2    3       69.0
99        1    2    3       77.0
100       1    2    3       82.0
101       1    3    3       63.0
102       1    3    3       83.5
103       1    3    3       73.5
104       1    3    3       76.5
105       1    3    3       71.0
106       2    1    3       94.5
107       2    1    3       88.0
108       2    1    3       77.5
109       2    1    3       66.0
110       2    1    3       92.0
111       2    2    3       83.0
112       2    2    3       67.0
113       2    2    3       65.5
114       2    2    3       61.5
115       2    2    3       61.0
116       2    3    3       54.5
117       2    3    3       61.0
118       2    3    3       62.5
119       2    3    3       75.5
120       2    3    3       80.5
121       1    1    4       69.5
122       1    1    4       90.0
123       1    1    4       82.0
124       1    1    4       86.0
125       1    1    4       89.0
126       1    2    4       82.5
127       1    2    4       71.0
128       1    2    4       70.5
129       1    2    4       79.5
130       1    2    4       80.5
131       1    3    4       63.5
132       1    3    4       84.0
133       1    3    4       75.5
134       1    3    4       75.0
135       1    3    4       76.0
136       2    1    4       83.5
137       2    1    4       86.0
138       2    1    4       78.0
139       2    1    4       64.5
140       2    1    4       90.5
141       2    2    4       82.5
142       2    2    4       68.0
143       2    2    4       64.5
144       2    2    4       62.5
145       2    2    4       61.5
146       2    3    4       53.5
147       2    3    4       61.0
148       2    3    4       62.5
149       2    3    4       75.5
150       2    3    4       80.5
151       1    1    5       69.5
152       1    1    5       87.5
153       1    1    5       80.5
154       1    1    5       83.0
155       1    1    5       83.0
156       1    2    5       80.0
157       1    2    5       71.5
158       1    2    5       71.5
159       1    2    5       80.0
160       1    2    5       79.5
161       1    3    5       61.0
162       1    3    5       81.5
163       1    3    5       71.0
164       1    3    5       71.0
165       1    3    5       70.0
166       2    1    5       81.5
167       2    1    5       83.5
168       2    1    5       73.5
169       2    1    5       61.5
170       2    1    5       87.5
171       2    2    5       80.5
172       2    2    5       65.0
173       2    2    5       59.5
174       2    2    5       60.0
175       2    2    5       60.0
176       2    3    5       52.0
177       2    3    5       59.0
178       2    3    5       61.5
179       2    3    5       75.0
180       2    3    5       79.0
181       1    1    6       72.0
182       1    1    6       92.5
183       1    1    6       84.0
184       1    1    6       88.0
185       1    1    6       87.5
186       1    2    6       82.0
187       1    2    6       74.0
188       1    2    6       74.0
189       1    2    6       83.0
190       1    2    6       81.0
191       1    3    6       63.0
192       1    3    6       82.5
193       1    3    6       75.0
194       1    3    6       72.0
195       1    3    6       75.0
196       2    1    6       85.5
197       2    1    6       89.5
198       2    1    6       79.5
199       2    1    6       63.0
200       2    1    6       94.5
201       2    2    6       85.5
202       2    2    6       68.5
203       2    2    6       60.5
204       2    2    6       65.5
205       2    2    6       62.0
206       2    3    6       54.5
207       2    3    6       61.0
208       2    3    6       62.0
209       2    3    6       80.0
210       2    3    6       79.0
211       1    1    7       77.5
212       1    1    7       98.5
213       1    1    7       88.5
214       1    1    7       92.5
215       1    1    7       92.0
216       1    2    7       85.0
217       1    2    7       77.5
218       1    2    7       77.5
219       1    2    7       85.0
220       1    2    7       82.0
221       1    3    7       67.5
222       1    3    7       83.5
223       1    3    7       81.0
224       1    3    7       72.5
225       1    3    7       76.5
226       2    1    7       93.5
227       2    1    7       94.5
228       2    1    7       85.0
229       2    1    7       64.5
230       2    1    7       99.5
231       2    2    7       90.0
232       2    2    7       72.0
233       2    2    7       62.0
234       2    2    7       70.0
235       2    2    7       63.0
236       2    3    7       57.5
237       2    3    7       62.5
238       2    3    7       62.5
239       2    3    7       83.5
240       2    3    7       80.0

The echo: false option disables the printing of code (only output is displayed).