diff --git a/R/gen_error.R b/R/gen_error.R
index f455dee6b6ade2cafe97a19b195a9832a9b6329a..d0ab68c22cd4fa3da3bbd2aee5139c06ddafb275 100644
--- a/R/gen_error.R
+++ b/R/gen_error.R
@@ -9,6 +9,7 @@
 #' @param FUN_predict function for the prediciton
 #' @param error_type chosen error type 'r','s','sr','or','sor','isor'
 #' @param par parameters of the input, output, structure, and remnant errors
+#' @param seed seed for reproducibility
 #' @keywords aqoslogen
 #' @export
 #' @examples
@@ -17,11 +18,14 @@
 #' }
 #
 gen_error <- function(dataset,target,predictors,FUN_fit,FUN_predict,error_type=c("r","or","sor","isor"),par=list(
-      phi_f1 = list(mean=5.64/1.96,sd=10),
-      phi_g1 = list(mean=5.64/1.96,sd=10),
-      phi_h1 = list(mean=0.1,sd=0.1),
-      beta_r  = list(mean=1,sd=2),
-      sig_r  = list(df=1,scale=0.2))) {
+      phi_f1 = list(4.35/100/1.96),
+      phi_g1 = list(4.35/100/1.96),
+      phi_h1 = list(0.1),
+      beta_r  = list(1),
+      sig_r  = list(1)),
+      seed=NULL) {
+
+  set.seed(seed) # for reproducibility
 
   res <- switch(error_type,