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.
