parcial1.R

setwd("~/Dropbox/Clases_UPJ/Estudios de Vegetación/Ciclo 1910/Examen_1")
library(BiodiversityR)

## Loading required package: tcltk

## Warning in system2("/usr/bin/otool", c("-L", shQuote(DSO)), stdout = TRUE):
## running command ''/usr/bin/otool' -L '/Library/Frameworks/R.framework/
## Resources/library/tcltk/libs//tcltk.so'' had status 1

## Loading required package: vegan

## Loading required package: permute

## Loading required package: lattice

## This is vegan 2.5-4

## BiodiversityR 2.11-1: Use command BiodiversityRGUI() to launch the Graphical User Interface; 
## to see changes use BiodiversityRGUI(changeLog=TRUE, backward.compatibility.messages=TRUE)

dapab <- read.csv("dap_ab.csv", header=T, sep=",")
attach(dapab)
##Los cambios principales de la sucesión involucran un incremento progresivo en la riqueza y la estructura de la vegetación
## Histogramas para el DAP (dap) 
hist(dap[suc=="Temprana"], breaks=10, prob=T,ylim = c(0,12), xlim = c(0,1))
lines(density(dap[suc=="Temprana"]))

hist(dap[suc=="Intermedia"], breaks=10, prob=T,ylim = c(0,12), xlim = c(0,1))
lines(density(dap[suc=="Intermedia"]))

hist(dap[suc=="Tardia"], breaks=seq(0,1, by=0.1), prob=T,ylim = c(0,12), xlim = c(0,1))
lines(density(dap[suc=="Tardia"]))

## Histogramas para el área basal (ab)
hist(ab[suc=="Temprana"], breaks=10, prob=T,ylim = c(0,20), xlim = c(0,0.6))

hist(ab[suc=="Intermedia"], breaks=10, prob=T,ylim = c(0,20), xlim = c(0,0.6))

hist(ab[suc=="Tardia"], breaks=seq(0,0.6, by=0.05), prob=T,ylim = c(0,20), xlim = c(0,0.6))

detach(dapab)
ivisp <- read.csv("ivi.csv", header=T, sep=",")
attach(ivisp)
ivi.temp=importancevalue(ivisp, site='trans', species='sp', count='abu', basal='ab',   factor='suc', level='Temprana')
write.csv (ivi.temp, file = "ivi_temp.csv")
ivi.inter=importancevalue(ivisp, site='trans', species='sp', count='abu', basal='ab',   factor='suc', level='intermedia')
write.csv (ivi.inter, file = "ivi_inter.csv")
ivi.tard=importancevalue(ivisp, site='trans', species='sp', count='abu', basal='ab',   factor='suc', level='Tardia')
write.csv (ivi.tard, file = "ivi_tard.csv")
detach(ivisp)
ivifam <- read.csv("ivif_orig.csv", header=T, sep=",")
attach(ivifam)
ivif.temp=importancevalue(ivifam, site='trans', species='fam', count='abun', basal='ab',   factor='suc', level='temprana')
write.csv (ivi.temp, file = "ivif_temp.csv")
ivif.inter=importancevalue(ivifam, site='trans', species='fam', count='abun', basal='ab',   factor='suc', level='intermedia')
write.csv (ivi.inter, file = "ivif_inter.csv")
ivif.tard=importancevalue(ivifam, site='trans', species='fam', count='abun', basal='ab',   factor='suc', level='tardia')
write.csv (ivi.tard, file = "ivif_tard.csv")