Olá Gilenio, seria algo assim:

### simula dados 

n <- 1000; beta <- c(-3,1)
y <- integer(n)
y[1] <- rbinom(1,1,.5)
for (i in 2:n)
  y[i] <- rbinom(1, 1, 1/(1+exp(-beta[1]-beta[2]*y[i-1])))

### estima modelo 

m <- glm(y ~ y2, family=binomial, data=data.frame(y=y[2:n], y2=y[2:n-1]))
summary(m)

 
Elias T. Krainski

De: Gilenio Borges Fernandes <gilenio@ufba.br>
Para: r-br@listas.c3sl.ufpr.br
Enviadas: Sábado, 14 de Julho de 2012 20:52
Assunto: [R-br] Ajuda para Simulação de variáveis binárias correlacionadas no R

Simulação de variáveis binárias correlacionadas no R é bem popular, usando os pacotes bindata, bild, binarySimCLF, BinNor e binGroup. Mas eu estou precisando simular variáveis binárias correlacionadas num modelo auto regressivo de primeira ordem de regressão logística, a partir dos parâmetros do preditor linear, incluindo o referente à variável dependente defasada (Yi-1). Alguém poderia me indicar algum pacote ou função no R para este fim?
Desde já, agradeço
--
Gilenio Borges Fernandes
Universidade Federal da Bahia
Instituto de Matemática
Departamento de Estatística
Av. Adhemar de Barros, s/n – Ondina.
40.170-110 - Salvador - BA, Brasil
Tel.: (071)3283-6280/6336  Fax:  (071)3283-6276
URL: http://lattes.cnpq.br/6764860618464860


_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.