## Script: 09_1varEstadoParces.R                  /
## Sobre:  calculo de variables a nivel agregado  /
##  por unidad muestral. Calculo de N y G.        /
##
## Profesor: Christian Salas Eljatib
## E-mail: cseljatib AT gmail DOT com
## Web: https://eljatib.com
#====================================================
library(datana)
data(radiatapl)
#?radiatapl

df<-radiatapl
#revisando los datos
str(df)

#cuantas parcelas hay?
table(df$plot)

head(df)
tail(df)
#empleemos una sola

head(df)
#primero calculemos la densidad
df$fe <- 10000/150
#num.arboles por parcela
nha.ppar<-tapply(df$fe,df$plot,sum)
nha.ppar

#luego el area basal de cada arbol
df$g <- (pi/40000)*df$dbh^2
#el area basal de cada arbol expandida a la hectarea
df$garb.ha <- df$g * df$fe

#area basal por parcela
gha.ppar<-tapply(df$garb.ha,df$plot,sum)
gha.ppar

##Variables de estado de rodal por parcela
nha<-nha.ppar;gha<-gha.ppar
varest.ppar<-data.frame(nha,gha)
varest.ppar

#calculo del diam. del arbol de area basal media
varest.ppar$dg<- sqrt((varest.ppar$gha/varest.ppar$nha)*(40000/pi))
varest.ppar

##guardando las varias de estado por parcela
write.csv(varest.ppar, file="varEstaParcelas.csv",row.names = F)




#╔═════════════════╗
#║ Fin del script! ║
#║ Atte.           ║
#║ El profesor     ║
#╚═════════════════╝