... and presumably the probabilities for class memberships are given by your function of the explanatory variable. So you just plug in your (30?) values of x, no?
-- Bert Bert Gunter "Data is not information. Information is not knowledge. And knowledge is certainly not wisdom." -- Clifford Stoll On Mon, Jun 29, 2015 at 8:50 AM, Christoph Scherber <csche...@gwdg.de> wrote: > Dear Don, > > True, but what I need is a sample *for a given explanatory variable*, where > class memberships are > defined in response to an underlying variable x (and known regression > parameters a, b). > > Have you had a look at the code and do you see a way to modify it for > rmultinom? > > Thanks a lot and best wishes > Christoph > > ## > > ## > set.seed(0) > x = sort(runif(1000,1,4)) > > a = 1.2 > b = 1.2 > y_det= a+b*poly(x,3) > > # for a normal distribution, I would use: > ynorm = rnorm(1000,y_det) > > # but how would I do it for a multinomial? > > ymult=rmultinom(.....) > > > > > > > > Am 29/06/2015 um 17:36 schrieb MacQueen, Don: >> See the combinat package: >> >> combinat::rmultinomial >> Generate random samples from multinomial distributions >> >> >> -Don >> > > > -- > PD Dr. Christoph Scherber > Senior Lecturer > DNPW, Agroecology > University of Goettingen > Grisebachstrasse 6 > 37077 Goettingen > Germany > telephone +49 551 39 8807 > facsimile +49 551 39 8806 > www.gwdg.de/~cscherb1 > > ______________________________________________ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > https://stat.ethz.ch/mailman/listinfo/r-help > PLEASE do read the posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
