# NOT RUN {
#Load data
data(presence)
#Set seed
set.seed(9842)
#Hyperparameters for each prior distribution
gamma <-0.01
alpha0<-0.01
alpha1<-0.01
#Execute the LDA for the Binomial entry
res<-rlda.bernoulli(presence, 10, alpha0, alpha1, gamma,
5000, TRUE, FALSE)
#Plot the results
plot(res)
# }
Run the code above in your browser using DataLab