update rng

snaketron · snaketron · commit fa486a46da0e · 2024-02-14T15:16:29.000+01:00
diff --git a/inst/stan/dgu_rep.stan b/inst/stan/dgu_rep.stan
@@ -11,14 +11,11 @@ functions {
   }
   
   int zibb_rng(int y, int n, real mu, real phi, real kappa) {
-    if(y==0) {
-      if (bernoulli_rng(kappa) == 1) {
-        return (0);
-      } else {
-        return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));
-      }
+    if (bernoulli_rng(kappa) == 1) {
+      return (0);
+    } else {
+      return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));
     }
-    return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));
   }
   
   real z_rng(real a, real b, real zi) {
diff --git a/inst/stan/gu.stan b/inst/stan/gu.stan
@@ -11,14 +11,11 @@ functions {
   }
   
   int zibb_rng(int y, int n, real mu, real phi, real kappa) {
-    if(y==0) {
-      if (bernoulli_rng(kappa) == 1) {
-        return (0);
-      } else {
-        return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));
-      }
+    if (bernoulli_rng(kappa) == 1) {
+      return (0);
+    } else {
+      return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));
     }
-    return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));
   }
   
   real z_rng(real a, real b, real zi) {
diff --git a/inst/stan/gu_rep.stan b/inst/stan/gu_rep.stan
@@ -11,14 +11,11 @@ functions {
   }
   
   int zibb_rng(int y, int n, real mu, real phi, real kappa) {
-    if(y==0) {
-      if (bernoulli_rng(kappa) == 1) {
-        return (0);
-      } else {
-        return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));
-      }
+    if (bernoulli_rng(kappa) == 1) {
+      return (0);
+    } else {
+      return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));
     }
-    return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));
   }
   
   real z_rng(real a, real b, real zi) {

Original file line number	Diff line number	Diff line change
`@@ -11,14 +11,11 @@ functions {`
`11`	`11`	`}`
`12`	`12`
`13`	`13`	`int zibb_rng(int y, int n, real mu, real phi, real kappa) {`
`14`		`- if(y==0) {`
`15`		`- if (bernoulli_rng(kappa) == 1) {`
`16`		`- return (0);`
`17`		`- } else {`
`18`		`- return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));`
`19`		`- }`
	`14`	`+ if (bernoulli_rng(kappa) == 1) {`
	`15`	`+ return (0);`
	`16`	`+ } else {`
	`17`	`+ return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));`
`20`	`18`	`}`
`21`		`- return (beta_binomial_rng(n, mu * phi, (1 - mu) * phi));`
`22`	`19`	`}`
`23`	`20`
`24`	`21`	`real z_rng(real a, real b, real zi) {`