search for: rdirichlet

Hi all, I have a problem using rdirichlet{gtools}. For Dir( a1, a2, ..., a_n), its mode can be found at $( a_i -1)/ ( \sum_{i}a_i - n)$; The means are $a_i / (\sum_{i} a_i ) $; I tried to study the above properties using rdirichlet from gtools. The code are: ############## library(gtools) alpha = c(1,3,9) #totoal=13 mean.expect = c(1...

...mpiled windows package has also been sent to Brian Ripley for inclusion in the windows package directory. This version of gregmisc fixes a number of bugs and introduces several new functions. Highlights of this version: =========================== New Features: - Introduction of pdirichlet() and rdirichlet() functions. - Addition of fast.prcomp() and fast.svd(), are more efficient than the built-in versions on very wide matrixes. (Short term solution until patches are integrated into the R source.) - Add permute() function as a convenience wrapper for sample() Bug Fixes: - Fixed bug in barplo...

Hola comunidad, Les consulto mi duda. Quisiera generar N números aleatorios tal que su suma esté predeterminada a un valor k. Sé que con la distribución Dirichlet, por ejemplo, se podrían obtener números aleatorios tal que su suma es la unidad. Por ejemplo MCMCpack::rdirichlet(1,rep(10,7)) genera 7 números tal que su suma es 1 (el parámetro rep(10,7) es cualquier cosa). Pero me pregunto si se podrá realizar algo parecido variando la suma (no necesariamente 1) y que los valores no solo sean positivos (como en la Dirichlet). ¡Gracias de antemano! -- «Pídeles sus títul...

...k interface ... done Warning in download.packages(unique(pkgs), destdir = tmpd, available = available, : no package 'bayesm' at the repositories The package is certainly there. It could be my version of linux, I know. All I want to do is sample from a dirichlet prior, so I want rdirichlet and ddirichlet commands basically. If anyone can help or suggest an alternate package that would be great. E -- View this message in context: http://www.nabble.com/bayesm-package-not-downloading-via-any-mirror-or-repository-tp20906495p20906495.html Sent from the R help mailing list archive at Nab...

...gtools::invalid gtools::keywords gtools::logit gtools::logratio2foldchange gtools::mixedorder gtools::mixedsort gtools::odd gtools::permutations gtools::permute gtools::quantcut gtools::rdirichlet gtools::running gtools::scat gtools::setTCPNoDelay gtools::smartbind gtools::sprint gtools::strmacro I am new to R, and I googled everything about the error message. I don't really understand what is "requires the gtools...

...nts of a Vector combine Combine R Objects With a Column Labeling the Source contrast.lm Compute and test arbitrary contrasts for regression objects (Depreciated, use 'fit.contrasts' instead.) rdirichlet Functions for the Dirichlet Distribution estimable Compute and test estimable linear functions of the fitted coefficients (including contrasts) of regression objects factorial Compute factorial function fast.pr...

...trivial example was a little too trivial and is causing some confusion. What's happening in the real function I'm writing is... 1. In R: Draw tens-of-thousands of times from a handful to Gamma RVs with different parameters to initialize some variables. (Technically, I'm calling gtools::rdirichlet which calls stats::rgamma) 2. Transfer the initialized variables to a function in C++ 3. In C++: Draw millions of times from a Categorical(p) distribution, where "p" is recalculated after each draw based on the current state of the RVs in my system. (The heart of this is actually a Unifor...

Thank you for this. I'd like to be sure I understand the intuition correctly. Is the following true from what you said? I can just fix the seed at the C++ level and the results will still be (pseudo) random because the initialization at the R level is (pseudo) random. On Thu, Jul 30, 2020 at 3:36 PM Duncan Murdoch <murdoch.duncan at gmail.com> wrote: > I wouldn't trust the C++