-
Notifications
You must be signed in to change notification settings - Fork 1
Expand file tree
/
Copy pathscurvy_analysis_script.R
More file actions
13 lines (13 loc) · 910 Bytes
/
scurvy_analysis_script.R
File metadata and controls
13 lines (13 loc) · 910 Bytes
1
2
3
4
5
6
7
8
9
10
11
12
13
scurvy <- read.table("scurvy.csv", header = TRUE, sep = ',')
summary(scurvy)
str(scurvy)
scurvy$treatment <- factor(scurvy$treatment)
scurvy$gum_rot_d6 <- factor(scurvy$gum_rot_d6)
scurvy$skin_sores_d6 <- factor(scurvy$skin_sores_d6)
table(scurvy$study_id, scurvy$gum_rot_d6)
table (scurvy$study_id, scurvy$weakness_of_the_knees_d6)
boxplot(scurvy$gum_rot_d6 ~ scurvy$treatment, data = scurvy, cex.axis = 1.2, xlab = "", ylab = "")
stripchart(scurvy$gum_rot_d6 ~ scurvy$treatment, data = scurvy, add = TRUE, vertical = TRUE, method = "jitter", pch = 21, col = "maroon", bg = "bisque")
scurvy$treatment <- factor(scurvy$treatment)
boxplot(scurvy$weakness_of_the_knees_d6 ~ scurvy$treatment, data = scurvy, cex.axis = 1.2, xlab = "", ylab = "")
stripchart(scurvy$weakness_of_the_knees_d6 ~ scurvy$treatment, data = scurvy, add = TRUE, vertical = TRUE, method = "jitter", pch = 21, col = "maroon", bg = "bisque")