Heliconia <- read.csv("Heliconia.csv")
Heliconia
## variety length
## 1 bihai 47.12
## 2 bihai 46.75
## 3 bihai NA
## 4 bihai 47.12
## 5 bihai 46.67
## 6 bihai 47.43
## 7 <NA> 46.44
## 8 bihai 46.64
## 9 bihai 48.07
## 10 bihai 48.34
## 11 bihai 48.15
## 12 bihai 50.26
## 13 bihai 50.12
## 14 bihai 46.34
## 15 bihai 46.94
## 16 bihai 48.36
## 17 red 41.90
## 18 red 42.01
## 19 red 41.93
## 20 red 43.09
## 21 red 41.47
## 22 red 41.69
## 23 red 39.78
## 24 red 40.57
## 25 red 39.63
## 26 red 42.18
## 27 red 40.66
## 28 red 37.87
## 29 red 39.16
## 30 red 37.40
## 31 red 38.20
## 32 red 38.07
## 33 red 38.10
## 34 red 37.97
## 35 red 38.79
## 36 red 38.23
## 37 red 38.87
## 38 red 37.78
## 39 red 38.01
## 40 yellow 36.78
## 41 yellow 37.02
## 42 yellow 36.52
## 43 yellow 36.11
## 44 yellow 36.03
## 45 yellow 35.45
## 46 yellow 38.13
## 47 yellow 37.10
## 48 yellow 35.17
## 49 yellow 36.82
## 50 yellow 36.66
## 51 yellow 35.68
## 52 yellow 36.03
## 53 yellow 34.57
## 54 yellow 34.63
summary(Heliconia)
## variety length
## Length:54 Min. :34.57
## Class :character 1st Qu.:37.10
## Mode :character Median :39.16
## Mean :40.96
## 3rd Qu.:46.44
## Max. :50.26
## NA's :1
mean(Heliconia$length)
## [1] NA
sd(Heliconia$length)
## [1] NA
mean(Heliconia$length, na.rm=TRUE)
## [1] 40.95868
sd(Heliconia$length, na.rm=TRUE)
## [1] 4.714614
library(RcmdrMisc)
## Loading required package: car
## Loading required package: carData
## Loading required package: sandwich
numSummary(Heliconia$length, groups=Heliconia$variety)
## mean sd IQR 0% 25% 50% 75% 100% data:n
## bihai 47.73643 1.2352490 1.495 46.34 46.7975 47.275 48.2925 50.26 14
## red 39.71130 1.7987630 3.495 37.40 38.0850 39.160 41.5800 43.09 23
## yellow 36.18000 0.9753241 1.235 34.57 35.5650 36.110 36.8000 38.13 15
## data:NA
## bihai 1
## red 0
## yellow 0
with(Heliconia, numSummary(length, groups=variety))
## mean sd IQR 0% 25% 50% 75% 100% data:n
## bihai 47.73643 1.2352490 1.495 46.34 46.7975 47.275 48.2925 50.26 14
## red 39.71130 1.7987630 3.495 37.40 38.0850 39.160 41.5800 43.09 23
## yellow 36.18000 0.9753241 1.235 34.57 35.5650 36.110 36.8000 38.13 15
## data:NA
## bihai 1
## red 0
## yellow 0
hist(Heliconia$length)
hist(Heliconia$length, col = "red")
boxplot(length ~ variety, data=Heliconia)
