wants <- c("DescTools", "robustbase")
has <- wants %in% rownames(installed.packages())
if(any(!has)) install.packages(wants[!has])
[1] 26.26667
[1] 5.125102
$cov
[,1]
[1,] 21.88889
$center
[1] 24.33333
$n.obs
[1] 6
[1] 21.88889
[1] 4.678556
[1] 17.2
[1] 4.147288
0% 25% 50% 75% 100%
17.00 21.50 24.00 28.75 30.00
[1] 7.25
[1] 4
[1] 6.6717
[1] 6.792788
[1] 4.865323
fac <- factor(c("C", "D", "A", "D", "E", "D", "C", "E", "E", "B", "E"),
levels=c(LETTERS[1:5], "Q"))
P <- nlevels(fac)
(Fj <- prop.table(table(fac)))
fac
A B C D E Q
0.09090909 0.09090909 0.18181818 0.27272727 0.36363636 0.00000000
First, calculate Shannon index, then diversity measure.
[1] 0.8193845
[1] -0.08611387
[1] -1.772568
R markdown - markdown - R code - all posts