analysis.r

In eunices/brve: Binary Regression Varying Effort Modelling

library(brve)

set.seed(2020)

# Modelling data
# y <- rbinom(10, 1, 0.5)
# x1 <- rnorm(10)
# x2 <- rnorm(10)
# effort <- runif(10, 0.5, 1.5)
# frogs <- data.frame(x1=x1, x2=x2, y=y, effort=effort)

# Test data
x1 <- rnorm(20)
x2 <- rnorm(20)
effort_new <- runif(20, 0.5, 1.5)
df_new <- data.frame(x1=x1, x2=x2, effort=effort_new)

# Modelling
data(frogs)
mod <- brve(y~x1+x2, frogs$effort, frogs)
predict(mod, df_new, "abundance", se.fit = FALSE)
plot(mod)
summary(mod)