Added randdir

R/admix1.R

@@ -693,15 +693,4 @@ admix1 <- function(tietue) {
 #     else
 #         str = ['0.' num2str(prob)];
 #     end;        
-# end;
-
-# %-------------------------------------------------
-
-# function svar=randdir(counts,nc)
-# % Käyttöesim randdir([10;30;60],3)
-
-# svar=zeros(nc,1);
-# for i=1:nc
-#    svar(i,1)=randga(counts(i,1),1);
-# end
-# svar=svar/sum(svar);
+# end;

R/randdir.R

@@ -0,0 +1,15 @@
+#' @title Generates random numbers
+#' @return vector of length `nc` with r.v. realizations from Gamma(rate=1)
+#' @examples randdir(matrix(c(10, 30, 60), 3), 3)
+#' @param counts shape parameter
+#' @param nc number of rows on output
+#' @seealso randga
+#' @export
+randdir <- function (counts, nc) {
+    svar <- zeros(nc, 1)
+    for (i in 1:nc) {
+       svar[i, 1] = randga(counts[i, 1], 1)
+    }
+    svar <- svar / sum(svar)
+    return(svar)
+}