R Notebook

a<-c(5,4,8,7,5,6,5,4,1,5,3,8,9)
b<-c(5,8,7,4,8,9,5,8,6,5,4,8,1)

summary(a)
summary(b)

sd(a)
median(a)

boxplot(a)

library(car)
data("mtcars")
scatterplot.matrix(mtcars)
head(mtcars)
mtcars[,1:4]->mtcars2
head(mtcars2)
scatterplotMatrix(mtcars2)

Biodiversidad

library(vegan)

data("dune")

head(dune)

dim(dune)

rarecurve (dune, col="blue")

dune[1:6, 1:4]

#Abundancia total de cada muestra
apply(dune, 1, sum)

#Abundancia total por especie
apply(dune, 2, sum)

specnumber(dune)

dune.shannon<-diversity(dune)
dune.shannon

data("dune.env")
head(dune.env)
attach(dune.env)

reg1<-lm(dune.shannon~ Management)
summary(reg1)

library(visreg)
visreg(reg1)