> # Visão do início dos dados
> head(galia)
plot species dbh
1 1 Boehmeria caudata 12.0
2 1 Cabralea canjerana 19.1
3 1 Boehmeria caudata 7.1
4 1 Sapium glandulatum 57.2
5 1 Syagrus romanzoffiana 26.7
6 1 Cupania vernalis 8.5
>
> # Cálculo da densidade
> dens = aggregate( galia$dbh, list(sp=galia$species), length )
> names(dens)[2] = "densidade" # acrescenta o nome correto à variável "x"
> head(dens)
sp densidade
1 Acacia polyphylla 1
2 Actinostemon concepcionis 1
3 Actinostemon concolor 16
4 Aspidosperma polyneuron 3
5 Boehmeria caudata 2
6 Cabralea canjerana 2
>
> # Cálculo da Dominância
> dom = aggregate( (pi/4)*galia$dbh^2, list(sp=galia$species), sum )
> names(dom)[2] = "dominancia" # acrescenta o nome correto à variável "x"
> head(dom)
sp dominancia
1 Acacia polyphylla 22.90221
2 Actinostemon concepcionis 51.52997
3 Actinostemon concolor 470.36711
4 Aspidosperma polyneuron 1218.03474
5 Boehmeria caudata 152.68926
6 Cabralea canjerana 307.75827
>
> # Cálculo da Frequência
> freq = aggregate( galia$plot, list(sp=galia$species), function(x) length(unique(x)) )
> names(freq)[2] = "frequencia" # acrescenta o nome correto à variável "x"
> head(freq)
sp frequencia
1 Acacia polyphylla 1
2 Actinostemon concepcionis 1
3 Actinostemon concolor 6
4 Aspidosperma polyneuron 3
5 Boehmeria caudata 1
6 Cabralea canjerana 2
>
> # Junta as variáveis num "data.frame"
> fitosoc = merge(merge(dens,dom),freq)
> head(fitosoc)
sp densidade dominancia frequencia
1 Acacia polyphylla 1 22.90221 1
2 Actinostemon concepcionis 1 51.52997 1
3 Actinostemon concolor 16 470.36711 6
4 Aspidosperma polyneuron 3 1218.03474 3
5 Boehmeria caudata 2 152.68926 1
6 Cabralea canjerana 2 307.75827 2
>
Notas:
Faça as devidas adaptações para os nomes das variáveis.
Não esqueça as transformações de unidades para que os valores de densidade e dominância seja apresentados em unidades apropriadas.
Valores relativos de densidade, dominância e frequência, bem como os famigerados IVI e IVC, seguem facilmente a partir desses cálculos.