de.bund.bfr.knime.testflows/workflows/InitializeParentsFlocks/FSK Creator (#1)/drop/model.R
##########################################################################################
# initialize_parents_animals = function(n_iter, Npos, Ntotal)
#
# DESCRIPTION:
# Initialize prevalence in parent flocks
#
# INPUT:
# n. iter integer number of iterations
# Npos integer number of positive animals in a flock
# Ntotal integer total number of animals in a flock
#
# then there is "elements"; it is not a part of the model it is a help varibale in order to write results in a file on disk, thus it does not appear in the list of independent variables below
# elements string can be "flocks" or "animals" and is used in the model script to create file names which store the end results to be used in subsequent model steps
#
# OUTPUT:
# Prevalence vector prevalence of ESBL among or within flock (depending on the values of Npos and Ntotal)
#
# @author: Carolina Plaza-Rodriguez, Federal Institute for Risk Assessment, Germany, 2015
# Alteration and modularisation by Guido correia Carreira, Federal Institute for Risk Assessment, Germany, 2016
#########################################################################################
initialize_parents_animals = function(n_iter, Npos, Ntotal){
Prevalence = 100*(rbeta(n_iter, shape1 = Npos + 1, shape2 = Ntotal - Npos + 1 ))
return(Prevalence)
}
result <- initialize_parents_animals(n_iter, Npos, Ntotal)