1
Emmilly Mel
18/02/2022
Estudos do R - reprodução da aula de descrição de comunidades biológicas
Utilizando a base de dados do estudo “Bird communities along elevation gradient in Yushan Mountain, Taiwan”
Essa é uma reprodução dos códigos usados na aula de descrição de comunidades biológicas para fins de estudo e prática no R.
Carregando todas as bibliotecas de dados que vou precisar usar.
ybirds.spe <- read.delim ('https://raw.githubusercontent.com/zdealveindy/anadat-r/master/data/ybirds_spe.txt', row.names = 1)
ybirds.env <- read.delim ('https://raw.githubusercontent.com/zdealveindy/anadat-r/master/data/ybirds_env.txt', row.names = 1)
ybirds.biomass <- read.delim ('https://raw.githubusercontent.com/zdealveindy/anadat-r/master/data/ybirds_biomass.txt', row.names = 1)Dando dados aos objetos
dim(ybirds.spe)## [1] 50 591. Quantas espécies em cada amostra?
specnumber(ybirds.spe)## 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26
## 24 25 28 25 24 20 16 16 19 28 29 29 16 16 16 19 26 18 26 24 27 30 25 28 20 22
## 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50
## 21 20 24 22 27 18 14 13 13 12 14 12 10 11 11 8 8 7 7 6 5 5 5 5“specnumber” é uma função para saber a quantidade de espécies no dado que estou pedindo.
2. Quandos indivíduos em cada amostra?
rowSums(ybirds.spe)## 1 2 3 4 5 6 7 8 9 10 11 12 13
## 58.23 67.12 56.92 57.93 58.97 61.65 57.00 61.61 49.18 54.18 51.29 54.09 36.22
## 14 15 16 17 18 19 20 21 22 23 24 25 26
## 36.51 37.40 37.91 57.17 40.03 62.21 47.98 43.11 50.54 45.05 51.24 63.55 60.07
## 27 28 29 30 31 32 33 34 35 36 37 38 39
## 65.47 54.08 59.73 59.18 60.81 49.05 46.26 39.21 33.08 36.07 43.28 49.86 37.32
## 40 41 42 43 44 45 46 47 48 49 50
## 34.88 29.12 33.21 29.11 25.36 15.76 9.51 9.89 10.86 10.74 7.31a função “rowSums” é pra saber quantos indivídos tem em cada amostra no dado.
3. Qual a abundância de cada espécie?
head(colSums(ybirds.spe))## ALA AWP BBM BBP BBW BGW
## 11.35 3.38 0.58 0.44 11.92 0.33Mostra somente algumas linhas, mas se removere o head(), mostra a tabela inteira:
colSums(ybirds.spe)## ALA AWP BBM BBP BBW BGW BLB BLS BNB BRB BRD
## 11.35 3.38 0.58 0.44 11.92 0.33 0.37 22.44 117.45 1.91 19.43
## BVB CLT CPO CSE FBB FFP FFW FHP FLT FMB FMF
## 20.95 79.06 0.92 0.11 0.29 68.82 113.34 0.14 67.07 30.22 390.08
## FMY FRF GBT GNB GTM ILT JAY JBR LBC LHC LSB
## 376.75 10.60 35.36 11.34 42.77 3.62 2.24 50.98 1.27 0.93 0.38
## MBW MKP MLB NTC NTH ORC PGW RBF RHB RHT RLT
## 0.22 5.12 4.02 1.74 14.85 2.11 0.33 40.96 32.53 36.03 7.37
## RSB STB SWB SWP THF VRF WBB WBW WBY WES WNB
## 0.35 93.58 30.39 0.50 17.83 50.20 13.97 3.03 20.48 65.33 124.15
## WRN WTB YBW YLT
## 71.24 31.52 34.36 7.26Funções espécificas do pacote “vegan” úteis para a análise de comunidade biológica.
1 - Curvas de acumulação de espécies
sp1 <- specaccum(ybirds.spe, "random")
sp1## Species Accumulation Curve
## Accumulation method: random, with 100 permutations
## Call: specaccum(comm = ybirds.spe, method = "random")
##
##
## Sites 1.0000 2.00000 3.00000 4.00000 5.00000 6.00000 7.00000 8.00000
## Richness 17.3900 28.20000 34.10000 38.44000 40.77000 43.16000 45.35000 46.82000
## sd 7.9986 8.18844 7.86952 6.49058 6.42887 5.80094 5.08588 4.63055
##
## Sites 9.0000 10.00000 11.00000 12.00000 13.00000 14.0000 15.00000 16.00000
## Richness 47.8800 48.65000 49.40000 50.07000 50.60000 51.2600 51.66000 52.16000
## sd 4.6521 4.37942 4.36237 3.99053 3.76856 3.4337 3.34911 3.19634
##
## Sites 17.00000 18.00000 19.00000 20.00000 21.00000 22.00000 23.00000
## Richness 52.45000 53.02000 53.39000 53.70000 53.99000 54.21000 54.44000
## sd 3.27332 3.06818 2.79211 2.55248 2.40578 2.43416 2.37546
##
## Sites 24.00000 25.00000 26.00000 27.0000 28.00000 29.00000 30.00000 31.0000
## Richness 54.71000 54.94000 55.21000 55.4800 55.63000 55.84000 56.17000 56.3100
## sd 2.27101 2.03911 2.04146 2.0225 2.00331 1.98337 1.64566 1.6249
##
## Sites 32.00000 33.00000 34.00000 35.00000 36.00000 37.00000 38.00000
## Richness 56.54000 56.75000 57.00000 57.14000 57.35000 57.46000 57.59000
## sd 1.47997 1.35121 1.27128 1.23108 1.14922 1.11392 1.07398
##
## Sites 39.00000 40.00000 41.00000 42.00000 43.00000 44.00000 45.00000
## Richness 57.73000 57.86000 57.98000 58.18000 58.28000 58.41000 58.53000
## sd 1.00358 0.93225 0.87594 0.74373 0.68283 0.63715 0.57656
##
## Sites 46.00000 47.00000 48.00000 49.00000 50
## Richness 58.62000 58.73000 58.82000 58.89000 59
## sd 0.56461 0.48938 0.38612 0.31447 0Plot
plot(sp1, ci.type = "poly", col = "red", lwd = 2, ci.lty = 0, ci.col = "lightblue")
boxplot(sp1, col = "yellow", add = TRUE, pch="+")
Curvas de rank-abundância
mod <- rad.lognormal(ybirds.spe[5,])## Warning in dpois(y, mu, log = TRUE): non-integer x = 17.910000## Warning in dpois(y, mu, log = TRUE): non-integer x = 7.640000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.140000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.980000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.970000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.880000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.230000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.910000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.490000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.380000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.180000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.830000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.640000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.090000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.070000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.030000mod##
## RAD model: Log-Normal
## Family: poisson
## No. of species: 24
## Total abundance: 58.97
##
## log.mu log.sigma Deviance AIC BIC
## 0.04091811 1.38924830 1.23576270 Inf Infplot(mod)
Esse é um plot lognormal
Para códigos com vários modelos:
mod2 <- radfit(ybirds.spe[1,])## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in nlm(canfun, p = p, x = x, rnk = rnk, logJ = logJ, wt = wt, hessian =
## TRUE, : NA/Inf substituido pelo máximo valor positivo## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in nlm(canfun, p = p, x = x, rnk = rnk, logJ = logJ, wt = wt, hessian =
## TRUE, : NA/Inf substituido pelo máximo valor positivo## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in nlm(canfun, p = p, x = x, rnk = rnk, logJ = logJ, wt = wt, hessian =
## TRUE, : NA/Inf substituido pelo máximo valor positivo## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in nlm(canfun, p = p, x = x, rnk = rnk, logJ = logJ, wt = wt, hessian =
## TRUE, : NA/Inf substituido pelo máximo valor positivo## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in nlm(mandelfun, p = p, x = x, rnk = rnk, off = off, family = fam, :
## NA/Inf substituido pelo máximo valor positivo## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in nlm(mandelfun, p = p, x = x, rnk = rnk, off = off, family = fam, :
## NA/Inf substituido pelo máximo valor positivo## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in nlm(mandelfun, p = p, x = x, rnk = rnk, off = off, family = fam, :
## NA/Inf substituido pelo máximo valor positivo## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in nlm(mandelfun, p = p, x = x, rnk = rnk, off = off, family = fam, :
## NA/Inf substituido pelo máximo valor positivo## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000## Warning in dpois(y, mu, log = TRUE): non-integer x = 13.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 9.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 6.300000## Warning in dpois(y, mu, log = TRUE): non-integer x = 4.460000## Warning in dpois(y, mu, log = TRUE): non-integer x = 3.020000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.950000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.550000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.340000## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 2.320000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.990000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.820000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.270000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 1.160000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.290000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.220000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000
## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.210000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.170000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.130000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.110000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.060000## Warning in dpois(y, mu, log = TRUE): non-integer x = 0.050000mod2##
## RAD models, family poisson
## No. of species 24, total abundance 58.23
##
## par1 par2 par3 Deviance AIC BIC
## Null 7.6732 Inf Inf
## Preemption 0.1 18.4224 Inf Inf
## Lognormal 0.17479 1.258 3.4343 Inf Inf
## Zipf 0.2892 -1.0678 8.9708 Inf Inf
## Mandelbrot 0.71017 -1.4058 1 6.5589 Inf Infplot(mod2)
### Funções com o “BiodiversityR”
library(BiodiversityR)## Carregando pacotes exigidos: tcltk## BiodiversityR 2.14-1: Use command BiodiversityRGUI() to launch the Graphical User Interface;
## to see changes use BiodiversityRGUI(changeLog=TRUE, backward.compatibility.messages=TRUE)carregando a biblioteca
data("dune")
data("dune.env")
bio <- rankabuncomp(dune, dune.env, factor = 'Management', return.data = TRUE, specnames = c(1:2), legend = FALSE)
bio## Grouping species labelit rank abundance proportion plower pupper accumfreq
## 1 BF Lolipere TRUE 1 18 15.4 3.6 27.2 15.4
## 2 BF Trifrepe TRUE 2 14 12.0 6.6 17.3 27.4
## 3 BF Scorautu FALSE 3 13 11.1 0.7 21.5 38.5
## 4 BF Poaprat FALSE 4 12 10.3 6.3 14.2 48.7
## 5 BF Poatriv FALSE 5 11 9.4 -9.9 28.7 58.1
## 6 BF Bromhord FALSE 6 8 6.8 -5.2 18.9 65.0
## 7 BF Achimill FALSE 7 7 6.0 -5.0 17.0 70.9
## 8 BF Planlanc FALSE 8 6 5.1 -6.9 17.1 76.1
## 9 BF Bracruta FALSE 9 6 5.1 -9.3 19.5 81.2
## 10 BF Bellpere FALSE 10 5 4.3 -4.0 12.5 85.5
## 11 BF Anthodor FALSE 11 4 3.4 -10.6 17.4 88.9
## 12 BF Elymrepe FALSE 12 4 3.4 -10.8 17.6 92.3
## 13 BF Vicilath FALSE 13 3 2.6 -4.6 9.8 94.9
## 14 BF Alopgeni FALSE 14 2 1.7 -5.4 8.8 96.6
## 15 BF Hyporadi FALSE 15 2 1.7 -6.3 9.7 98.3
## 16 BF Sagiproc FALSE 16 2 1.7 -6.3 9.7 100.0
## 17 HF Poatriv TRUE 1 24 11.3 8.4 14.1 11.3
## 18 HF Lolipere TRUE 2 20 9.4 3.7 15.1 20.7
## 19 HF Poaprat FALSE 3 17 8.0 4.9 11.1 28.6
## 20 HF Rumeacet FALSE 4 16 7.5 1.1 13.9 36.2
## 21 HF Planlanc FALSE 5 15 7.0 -0.7 14.7 43.2
## 22 HF Scorautu FALSE 6 14 6.6 5.2 8.0 49.8
## 23 HF Trifrepe FALSE 7 14 6.6 3.3 9.8 56.3
## 24 HF Bracruta FALSE 8 14 6.6 1.9 11.2 62.9
## 25 HF Elymrepe FALSE 9 10 4.7 -3.6 13.0 67.6
## 26 HF Anthodor FALSE 10 9 4.2 -0.8 9.2 71.8
## 27 HF Trifprat FALSE 11 9 4.2 -1.4 9.9 76.1
## 28 HF Alopgeni FALSE 12 8 3.8 -3.1 10.6 79.8
## 29 HF Juncarti FALSE 13 8 3.8 -2.8 10.3 83.6
## 30 HF Agrostol FALSE 14 7 3.3 -2.6 9.1 86.9
## 31 HF Achimill FALSE 15 6 2.8 -0.3 5.9 89.7
## 32 HF Juncbufo FALSE 16 6 2.8 -2.5 8.1 92.5
## 33 HF Eleopalu FALSE 17 4 1.9 -3.4 7.2 94.4
## 34 HF Sagiproc FALSE 18 4 1.9 -1.4 5.2 96.2
## 35 HF Bromhord FALSE 19 4 1.9 -1.4 5.1 98.1
## 36 HF Bellpere FALSE 20 2 0.9 -1.7 3.5 99.1
## 37 HF Ranuflam FALSE 21 2 0.9 -1.7 3.6 100.0
## 38 NM Scorautu TRUE 1 19 12.6 6.0 19.2 12.6
## 39 NM Bracruta TRUE 2 17 11.3 2.6 19.9 23.8
## 40 NM Agrostol FALSE 3 13 8.6 -1.2 18.4 32.5
## 41 NM Eleopalu FALSE 4 13 8.6 -1.4 18.6 41.1
## 42 NM Salirepe FALSE 5 11 7.3 -0.2 14.8 48.3
## 43 NM Trifrepe FALSE 6 11 7.3 -2.0 16.6 55.6
## 44 NM Anthodor FALSE 7 8 5.3 -3.8 14.4 60.9
## 45 NM Ranuflam FALSE 8 8 5.3 -1.2 11.8 66.2
## 46 NM Hyporadi FALSE 9 7 4.6 -3.8 13.1 70.9
## 47 NM Juncarti FALSE 10 7 4.6 -2.7 12.0 75.5
## 48 NM Callcusp FALSE 11 7 4.6 -2.8 12.1 80.1
## 49 NM Airaprae FALSE 12 5 3.3 -2.4 9.0 83.4
## 50 NM Planlanc FALSE 13 5 3.3 -2.6 9.2 86.8
## 51 NM Poaprat FALSE 14 4 2.6 -2.5 7.8 89.4
## 52 NM Comapalu FALSE 15 4 2.6 -1.8 7.1 92.1
## 53 NM Sagiproc FALSE 16 3 2.0 -2.9 6.9 94.0
## 54 NM Achimill FALSE 17 2 1.3 -2.4 5.0 95.4
## 55 NM Bellpere FALSE 18 2 1.3 -2.0 4.7 96.7
## 56 NM Empenigr FALSE 19 2 1.3 -1.9 4.6 98.0
## 57 NM Lolipere FALSE 20 2 1.3 -2.0 4.7 99.3
## 58 NM Vicilath FALSE 21 1 0.7 -1.0 2.3 100.0
## 59 SF Agrostol TRUE 1 28 13.7 7.8 19.7 13.7
## 60 SF Poatriv TRUE 2 28 13.7 6.3 21.1 27.5
## 61 SF Alopgeni FALSE 3 26 12.7 4.5 21.0 40.2
## 62 SF Lolipere FALSE 4 18 8.8 -2.1 19.8 49.0
## 63 SF Poaprat FALSE 5 15 7.4 0.5 14.2 56.4
## 64 SF Elymrepe FALSE 6 12 5.9 -1.0 12.8 62.3
## 65 SF Bracruta FALSE 7 12 5.9 0.8 10.9 68.1
## 66 SF Sagiproc FALSE 8 11 5.4 -0.8 11.5 73.5
## 67 SF Eleopalu FALSE 9 8 3.9 -6.3 14.1 77.5
## 68 SF Scorautu FALSE 10 8 3.9 1.4 6.4 81.4
## 69 SF Trifrepe FALSE 11 8 3.9 0.5 7.3 85.3
## 70 SF Juncbufo FALSE 12 7 3.4 -2.3 9.2 88.7
## 71 SF Bellpere FALSE 13 4 2.0 -0.9 4.8 90.7
## 72 SF Ranuflam FALSE 14 4 2.0 -1.3 5.2 92.6
## 73 SF Juncarti FALSE 15 3 1.5 -2.4 5.3 94.1
## 74 SF Callcusp FALSE 16 3 1.5 -2.4 5.3 95.6
## 75 SF Bromhord FALSE 17 3 1.5 -2.1 5.0 97.1
## 76 SF Rumeacet FALSE 18 2 1.0 -1.5 3.5 98.0
## 77 SF Cirsarve FALSE 19 2 1.0 -1.4 3.3 99.0
## 78 SF Achimill FALSE 20 1 0.5 -0.9 1.9 99.5
## 79 SF Chenalbu FALSE 21 1 0.5 -0.8 1.8 100.0
## logabun rankfreq
## 1 1.3 6.2
## 2 1.1 12.5
## 3 1.1 18.8
## 4 1.1 25.0
## 5 1.0 31.2
## 6 0.9 37.5
## 7 0.8 43.8
## 8 0.8 50.0
## 9 0.8 56.2
## 10 0.7 62.5
## 11 0.6 68.8
## 12 0.6 75.0
## 13 0.5 81.2
## 14 0.3 87.5
## 15 0.3 93.8
## 16 0.3 100.0
## 17 1.4 4.8
## 18 1.3 9.5
## 19 1.2 14.3
## 20 1.2 19.0
## 21 1.2 23.8
## 22 1.1 28.6
## 23 1.1 33.3
## 24 1.1 38.1
## 25 1.0 42.9
## 26 1.0 47.6
## 27 1.0 52.4
## 28 0.9 57.1
## 29 0.9 61.9
## 30 0.8 66.7
## 31 0.8 71.4
## 32 0.8 76.2
## 33 0.6 81.0
## 34 0.6 85.7
## 35 0.6 90.5
## 36 0.3 95.2
## 37 0.3 100.0
## 38 1.3 4.8
## 39 1.2 9.5
## 40 1.1 14.3
## 41 1.1 19.0
## 42 1.0 23.8
## 43 1.0 28.6
## 44 0.9 33.3
## 45 0.9 38.1
## 46 0.8 42.9
## 47 0.8 47.6
## 48 0.8 52.4
## 49 0.7 57.1
## 50 0.7 61.9
## 51 0.6 66.7
## 52 0.6 71.4
## 53 0.5 76.2
## 54 0.3 81.0
## 55 0.3 85.7
## 56 0.3 90.5
## 57 0.3 95.2
## 58 0.0 100.0
## 59 1.4 4.8
## 60 1.4 9.5
## 61 1.4 14.3
## 62 1.3 19.0
## 63 1.2 23.8
## 64 1.1 28.6
## 65 1.1 33.3
## 66 1.0 38.1
## 67 0.9 42.9
## 68 0.9 47.6
## 69 0.9 52.4
## 70 0.8 57.1
## 71 0.6 61.9
## 72 0.6 66.7
## 73 0.5 71.4
## 74 0.5 76.2
## 75 0.5 81.0
## 76 0.3 85.7
## 77 0.3 90.5
## 78 0.0 95.2
## 79 0.0 100.0Plor com ggplot
library(ggplot2)
library(ggrepel)plotgg1 <- ggplot(data=bio, aes(x = rank, y = abundance)) +
scale_x_continuous(expand=c(0, 1), sec.axis = dup_axis(labels=NULL, name=NULL)) +
scale_y_continuous(expand=c(0, 1), sec.axis = dup_axis(labels=NULL, name=NULL)) +
geom_line(aes(colour=Grouping), size=1) +
geom_point(aes(colour=Grouping, shape=Grouping), size=5, alpha=0.7) +
geom_text_repel(data=subset(bio, labelit == TRUE),
aes(colour=Grouping, label=species),
angle=45, nudge_x=1, nudge_y=1, show.legend=FALSE) +
scale_color_brewer(palette = "Set1") +
labs(x = "rank", y = "abundance", colour = "Management", shape = "Management")plotgg1
RA.data <- rankabuncomp(dune, y=dune.env, factor='Management', return.data = TRUE, specnames = c(1:10),legend=FALSE)
plotgg2 <- ggplot(data=RA.data, aes(x = rank, y = abundance)) +
scale_x_continuous(expand=c(0, 1), sec.axis = dup_axis(labels=NULL, name=NULL)) +
scale_y_continuous(expand=c(0, 1), sec.axis = dup_axis(labels=NULL, name=NULL)) +
geom_line(aes(colour=Grouping), size=1) +
geom_point(aes(colour=Grouping), size=5, alpha=0.7) +
geom_text_repel(data=subset(RA.data, labelit == TRUE),
aes(label=species),
angle=45, nudge_x=1, nudge_y=1, show.legend=FALSE) +
scale_color_brewer(palette = "Set1") +
facet_wrap(~ Grouping) +
labs(x = "rank", y = "abundance", colour = "Management")plotgg2