File:20linjaa.csv

cat("Loading necessary packages and functions.\n")

library(OpasnetBaseUtils)
library(reshape)
library(xtable)

cat("Loading data and transforming it into the right format.\n")

data <- opasnet.csv("1/12/20linjaa.csv")
data <- melt(data, id.var = 1:4)
names(data)[1] <- "transect"

cat("Showing the first rows of data.\n")

head(data)

q.wiki <- c(0,1,2,3,4)

out <- diversity(amount = data$value, species = data$variable, transect = data$transect)

cat("Gamma diversities for transects.\n")
print(xtable(out[[3]]), type = 'html', html.table.attributes = "class='sortable'")

out <- diversity.table(amount = data$value, species = data$variable, transect = data$transect)

cat("Summary of diversity indices.\n")
print(xtable(out), type = 'html', html.table.attributes = "class='sortable'")