library(dplyr)
## Warning: package 'dplyr' was built under R version 3.4.2
##
## Attaching package: 'dplyr'
## The following objects are masked from 'package:stats':
##
## filter, lag
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union
library(datasets)
data(ToothGrowth)
ToothGrowth
## len supp dose
## 1 4.2 VC 0.5
## 2 11.5 VC 0.5
## 3 7.3 VC 0.5
## 4 5.8 VC 0.5
## 5 6.4 VC 0.5
## 6 10.0 VC 0.5
## 7 11.2 VC 0.5
## 8 11.2 VC 0.5
## 9 5.2 VC 0.5
## 10 7.0 VC 0.5
## 11 16.5 VC 1.0
## 12 16.5 VC 1.0
## 13 15.2 VC 1.0
## 14 17.3 VC 1.0
## 15 22.5 VC 1.0
## 16 17.3 VC 1.0
## 17 13.6 VC 1.0
## 18 14.5 VC 1.0
## 19 18.8 VC 1.0
## 20 15.5 VC 1.0
## 21 23.6 VC 2.0
## 22 18.5 VC 2.0
## 23 33.9 VC 2.0
## 24 25.5 VC 2.0
## 25 26.4 VC 2.0
## 26 32.5 VC 2.0
## 27 26.7 VC 2.0
## 28 21.5 VC 2.0
## 29 23.3 VC 2.0
## 30 29.5 VC 2.0
## 31 15.2 OJ 0.5
## 32 21.5 OJ 0.5
## 33 17.6 OJ 0.5
## 34 9.7 OJ 0.5
## 35 14.5 OJ 0.5
## 36 10.0 OJ 0.5
## 37 8.2 OJ 0.5
## 38 9.4 OJ 0.5
## 39 16.5 OJ 0.5
## 40 9.7 OJ 0.5
## 41 19.7 OJ 1.0
## 42 23.3 OJ 1.0
## 43 23.6 OJ 1.0
## 44 26.4 OJ 1.0
## 45 20.0 OJ 1.0
## 46 25.2 OJ 1.0
## 47 25.8 OJ 1.0
## 48 21.2 OJ 1.0
## 49 14.5 OJ 1.0
## 50 27.3 OJ 1.0
## 51 25.5 OJ 2.0
## 52 26.4 OJ 2.0
## 53 22.4 OJ 2.0
## 54 24.5 OJ 2.0
## 55 24.8 OJ 2.0
## 56 30.9 OJ 2.0
## 57 26.4 OJ 2.0
## 58 27.3 OJ 2.0
## 59 29.4 OJ 2.0
## 60 23.0 OJ 2.0
teeth <- tbl_df(ToothGrowth)
teeth
## # A tibble: 60 x 3
## len supp dose
## <dbl> <fctr> <dbl>
## 1 4.2 VC 0.5
## 2 11.5 VC 0.5
## 3 7.3 VC 0.5
## 4 5.8 VC 0.5
## 5 6.4 VC 0.5
## 6 10.0 VC 0.5
## 7 11.2 VC 0.5
## 8 11.2 VC 0.5
## 9 5.2 VC 0.5
## 10 7.0 VC 0.5
## # ... with 50 more rows
str(teeth)
## Classes 'tbl_df', 'tbl' and 'data.frame': 60 obs. of 3 variables:
## $ len : num 4.2 11.5 7.3 5.8 6.4 10 11.2 11.2 5.2 7 ...
## $ supp: Factor w/ 2 levels "OJ","VC": 2 2 2 2 2 2 2 2 2 2 ...
## $ dose: num 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 0.5 ...
summarise(teeth)
## data frame with 0 columns and 0 rows
teeth
## # A tibble: 60 x 3
## len supp dose
## <dbl> <fctr> <dbl>
## 1 4.2 VC 0.5
## 2 11.5 VC 0.5
## 3 7.3 VC 0.5
## 4 5.8 VC 0.5
## 5 6.4 VC 0.5
## 6 10.0 VC 0.5
## 7 11.2 VC 0.5
## 8 11.2 VC 0.5
## 9 5.2 VC 0.5
## 10 7.0 VC 0.5
## # ... with 50 more rows
glimpse(teeth)
## Observations: 60
## Variables: 3
## $ len <dbl> 4.2, 11.5, 7.3, 5.8, 6.4, 10.0, 11.2, 11.2, 5.2, 7.0, 16....
## $ supp <fctr> VC, VC, VC, VC, VC, VC, VC, VC, VC, VC, VC, VC, VC, VC, ...
## $ dose <dbl> 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 0.5, 1.0, 1....
teeth_df<- table(teeth)
teeth.freq <- table(teeth$dose)
teeth.freq
##
## 0.5 1 2
## 20 20 20
teeth.freq <- table(teeth$supp)
teeth.freq
##
## OJ VC
## 30 30
teeth.freq <- table(teeth$len)
teeth.freq
##
## 4.2 5.2 5.8 6.4 7 7.3 8.2 9.4 9.7 10 11.2 11.5 13.6 14.5 15.2
## 1 1 1 1 1 1 1 1 2 2 2 1 1 3 2
## 15.5 16.5 17.3 17.6 18.5 18.8 19.7 20 21.2 21.5 22.4 22.5 23 23.3 23.6
## 1 3 2 1 1 1 1 1 1 2 1 1 1 2 2
## 24.5 24.8 25.2 25.5 25.8 26.4 26.7 27.3 29.4 29.5 30.9 32.5 33.9
## 1 1 1 2 1 4 1 2 1 1 1 1 1
teeth.relfreq <- teeth.freq / nrow(teeth)
teeth.relfreq
##
## 4.2 5.2 5.8 6.4 7 7.3
## 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667
## 8.2 9.4 9.7 10 11.2 11.5
## 0.01666667 0.01666667 0.03333333 0.03333333 0.03333333 0.01666667
## 13.6 14.5 15.2 15.5 16.5 17.3
## 0.01666667 0.05000000 0.03333333 0.01666667 0.05000000 0.03333333
## 17.6 18.5 18.8 19.7 20 21.2
## 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667 0.01666667
## 21.5 22.4 22.5 23 23.3 23.6
## 0.03333333 0.01666667 0.01666667 0.01666667 0.03333333 0.03333333
## 24.5 24.8 25.2 25.5 25.8 26.4
## 0.01666667 0.01666667 0.01666667 0.03333333 0.01666667 0.06666667
## 26.7 27.3 29.4 29.5 30.9 32.5
## 0.01666667 0.03333333 0.01666667 0.01666667 0.01666667 0.01666667
## 33.9
## 0.01666667
round.teeth.relfreq <- round (teeth.relfreq, digits = 2)
round.teeth.relfreq
##
## 4.2 5.2 5.8 6.4 7 7.3 8.2 9.4 9.7 10 11.2 11.5 13.6 14.5 15.2
## 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.02 0.03 0.03 0.03 0.02 0.02 0.05 0.03
## 15.5 16.5 17.3 17.6 18.5 18.8 19.7 20 21.2 21.5 22.4 22.5 23 23.3 23.6
## 0.02 0.05 0.03 0.02 0.02 0.02 0.02 0.02 0.02 0.03 0.02 0.02 0.02 0.03 0.03
## 24.5 24.8 25.2 25.5 25.8 26.4 26.7 27.3 29.4 29.5 30.9 32.5 33.9
## 0.02 0.02 0.02 0.03 0.02 0.07 0.02 0.03 0.02 0.02 0.02 0.02 0.02
cbind(round.teeth.relfreq)
## round.teeth.relfreq
## 4.2 0.02
## 5.2 0.02
## 5.8 0.02
## 6.4 0.02
## 7 0.02
## 7.3 0.02
## 8.2 0.02
## 9.4 0.02
## 9.7 0.03
## 10 0.03
## 11.2 0.03
## 11.5 0.02
## 13.6 0.02
## 14.5 0.05
## 15.2 0.03
## 15.5 0.02
## 16.5 0.05
## 17.3 0.03
## 17.6 0.02
## 18.5 0.02
## 18.8 0.02
## 19.7 0.02
## 20 0.02
## 21.2 0.02
## 21.5 0.03
## 22.4 0.02
## 22.5 0.02
## 23 0.02
## 23.3 0.03
## 23.6 0.03
## 24.5 0.02
## 24.8 0.02
## 25.2 0.02
## 25.5 0.03
## 25.8 0.02
## 26.4 0.07
## 26.7 0.02
## 27.3 0.03
## 29.4 0.02
## 29.5 0.02
## 30.9 0.02
## 32.5 0.02
## 33.9 0.02
library(dplyr)
library(ggplot2)
crosstab <- table(teeth$dose, teeth$len)
crosstab
##
## 4.2 5.2 5.8 6.4 7 7.3 8.2 9.4 9.7 10 11.2 11.5 13.6 14.5 15.2 15.5
## 0.5 1 1 1 1 1 1 1 1 2 2 2 1 0 1 1 0
## 1 0 0 0 0 0 0 0 0 0 0 0 0 1 2 1 1
## 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
##
## 16.5 17.3 17.6 18.5 18.8 19.7 20 21.2 21.5 22.4 22.5 23 23.3 23.6
## 0.5 1 0 1 0 0 0 0 0 1 0 0 0 0 0
## 1 2 2 0 0 1 1 1 1 0 0 1 0 1 1
## 2 0 0 0 1 0 0 0 0 1 1 0 1 1 1
##
## 24.5 24.8 25.2 25.5 25.8 26.4 26.7 27.3 29.4 29.5 30.9 32.5 33.9
## 0.5 0 0 0 0 0 0 0 0 0 0 0 0 0
## 1 0 0 1 0 1 1 0 1 0 0 0 0 0
## 2 1 1 0 2 0 3 1 1 1 1 1 1 1
library(ggthemes)
## Warning: package 'ggthemes' was built under R version 3.4.2
crosstab2 <- table(teeth$dose, teeth$supp)
crosstab2
##
## OJ VC
## 0.5 10 10
## 1 10 10
## 2 10 10
ftable(crosstab2)
## OJ VC
##
## 0.5 10 10
## 1 10 10
## 2 10 10
table(teeth$dose, teeth$len)
##
## 4.2 5.2 5.8 6.4 7 7.3 8.2 9.4 9.7 10 11.2 11.5 13.6 14.5 15.2 15.5
## 0.5 1 1 1 1 1 1 1 1 2 2 2 1 0 1 1 0
## 1 0 0 0 0 0 0 0 0 0 0 0 0 1 2 1 1
## 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
##
## 16.5 17.3 17.6 18.5 18.8 19.7 20 21.2 21.5 22.4 22.5 23 23.3 23.6
## 0.5 1 0 1 0 0 0 0 0 1 0 0 0 0 0
## 1 2 2 0 0 1 1 1 1 0 0 1 0 1 1
## 2 0 0 0 1 0 0 0 0 1 1 0 1 1 1
##
## 24.5 24.8 25.2 25.5 25.8 26.4 26.7 27.3 29.4 29.5 30.9 32.5 33.9
## 0.5 0 0 0 0 0 0 0 0 0 0 0 0 0
## 1 0 0 1 0 1 1 0 1 0 0 0 0 0
## 2 1 1 0 2 0 3 1 1 1 1 1 1 1
plot(teeth$dose)
plot(crosstab)
plot(crosstab2)
plot(teeth$dose)
xlab ="dose"
ylab = "supp"
main= "condition of teeth"
type="O"
plot(teeth$len)
xlab ="dose"
ylab = "supp"
main= "condition of teeth"
type="L'"
col= ("orange")
hist(teeth$len)
