[R-br] Ajuda para Simulação de variáveis binárias correlacionadas no R

[Elias T. Krainski]

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 em ufba.br>
>Para: r-br em 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 em 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.
>
>
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120716/0532d17a/attachment.html>