<html><body><div style="color:#000; background-color:#fff; font-family:times new roman, new york, times, serif;font-size:12pt"><div><span>Olá Gilenio, seria algo assim:</span></div><div><span><br></span></div><div><span>### simula dados </span></div><div><span><br></span></div><div><span>n <- 1000; beta <- c(-3,1)</span></div><div><span>y <- integer(n)</span></div><div><span>y[1] <- rbinom(1,1,.5)</span></div><div><span>for (i in 2:n)</span></div><div><span> y[i] <- rbinom(1, 1, 1/(1+exp(-beta[1]-beta[2]*y[i-1])))</span></div><div><span><br></span></div><div><span>### estima modelo </span></div><div><span><br></span></div><div><span>m <- glm(y ~ y2, family=binomial, data=data.frame(y=y[2:n], y2=y[2:n-1]))</span></div><div><span>summary(m)</span></div><div><span><br></span></div><div> </div><div>Elias T. Krainski<br><blockquote style="border-left: 2px solid rgb(16, 16, 255); margin-left: 5px; margin-top: 5px;
padding-left: 5px;"> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "> <div style="font-family: 'times new roman', 'new york', times, serif; font-size: 12pt; "> <div dir="ltr"> <font size="2" face="Arial"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Gilenio Borges Fernandes <gilenio@ufba.br><br> <b><span style="font-weight: bold;">Para:</span></b> r-br@listas.c3sl.ufpr.br <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Sábado, 14 de Julho de 2012 20:52<br> <b><span style="font-weight: bold;">Assunto:</span></b> [R-br] Ajuda para Simulação de variáveis binárias correlacionadas no R<br> </font> </div> <br><meta http-equiv="x-dns-prefetch-control" content="off"><div id="yiv466610480"><div class="yiv466610480MsoNoSpacing" style="text-align:justify;"><span style="font-size:12.0pt;">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 (Y<sub>i-1</sub>). Alguém poderia me indicar algum pacote ou função no
R para este fim?</span></div><div>Desde já, agradeço</div>-- <br>Gilenio Borges Fernandes<br>Universidade Federal da Bahia<br>Instituto de Matemática<br>Departamento de Estatística<br>Av. Adhemar de Barros, s/n – Ondina.<br>
40.170-110 - Salvador - BA, Brasil<br>Tel.: (071)3283-6280/6336 Fax: (071)3283-6276<br>URL: <a rel="nofollow" target="_blank" href="http://lattes.cnpq.br/6764860618464860">http://lattes.cnpq.br/6764860618464860</a><br><br>
</div><meta http-equiv="x-dns-prefetch-control" content="on"><br>_______________________________________________<br>R-br mailing list<br><a ymailto="mailto:R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br">R-br@listas.c3sl.ufpr.br</a><br><a href="https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br" target="_blank">https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br</a><br>Leia o guia de postagem (<a href="http://www.leg.ufpr.br/r-br-guia" target="_blank">http://www.leg.ufpr.br/r-br-guia</a>) e forneça código mínimo reproduzível.<br><br> </div> </div> </blockquote></div> </div></body></html>