[R-br] simulação de mutinomiais

[Walmes Zeviani]

Então é só escrever a função logística, conforme se queira, e determinar
essas probabilidades e simular. Um ponto de partida pode ser esse

logis <- function(x, A, B) 1/(1+exp(-(A+B*x)))
curve(logis(x, 1, 1), -9, 4)
curve(logis(x, 2.3, 1), add=TRUE)
curve(logis(x, 3.6, 1), add=TRUE)
x <- -2
px <- logis(-2, c(1,2.3,3.6), 1)
abline(v=x, h=px, lty=3)
probs <- diff(c(0,px,1))
probs
sum(probs)
amostr <- sample(1:4, size=1000, prob=probs, replace=TRUE)
prop.table(table(amostr))

À disposição.
Walmes.

==========================================================================
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W)
Departamento de Estatística - Universidade Federal do Paraná
fone: (+55) 41 3361 3573
VoIP: (3361 3600) 1053 1173
e-mail: walmes em ufpr.br
twitter: @walmeszeviani
homepage: http://www.leg.ufpr.br/~walmes
linux user number: 531218
==========================================================================
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120517/e41390d6/attachment.html>