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@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================