Bem os valores médios e da matriz de covariancia são provenientes dos dados reais, ou seja, calculei as médias das duas variaveis e matriz de covariancia e as utilizei na simulação.<br><br><div class="gmail_quote">Em 16 de junho de 2012 20:18, Sandro Sperandei <span dir="ltr"><<a href="mailto:ssperandei@hotmail.com" target="_blank">ssperandei@hotmail.com</a>></span> escreveu:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div dir="ltr">
<div dir="ltr">
<div style="font-size:12pt;font-family:'Calibri'">
<div>Natalia,</div>
<div> </div>
<div>Provavelmente, o problema está na sua matriz de correlação. Se vc fizer
cor(bnorm), verá que os dados gerados apresentam sempre correlação próxima a
0,19. De onde vieram os valores usados na matriz?</div>
<div> </div>
<div>S</div>
<div> </div>
<div> </div>
<div style="font-size:small;font-style:normal;text-decoration:none;font-family:'Calibri';display:inline;font-weight:normal">
<div style="FONT:10pt tahoma">
<div> </div>
<div style="BACKGROUND:#f5f5f5">
<div><b>From:</b> <a title="nsmbarreto@gmail.com" href="mailto:nsmbarreto@gmail.com" target="_blank">Natalia Martins</a> </div>
<div><b>Sent:</b> Saturday, June 16, 2012 7:39 PM</div>
<div><b>To:</b> <a title="R-br@listas.c3sl.ufpr.br" href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">R-br@listas.c3sl.ufpr.br</a> </div>
<div><b>Subject:</b> [R-br] simulação da normal bivariada</div></div></div>
<div> </div></div>
<div style="font-size:small;font-style:normal;text-decoration:none;font-family:'Calibri';display:inline;font-weight:normal"><div><div>Boa
noite pessoal,
<div>estou tentando simular dados diários de temperatura e radiação solar usando
uma normal bivariada e após a simulação eu tentei fazer uma reta de regressão
para verificar se minha simulação esta de acordo, ou seja, eu esperaria uma reta
pois os valores deveriam ser iguais. No entanto quando a faço o meu R^2 ajustado
é ruim (sei que tambem devo fazer testes e não apenas olhar o R), alguem pode me
ajudar, pois acredito que estou fazendo algo errado?</div>
<div> </div>
<div>a programação que estou utilizando segue a baixo:</div>
<div> </div>
<div>
<div># gerando os dados normal bivariada</div>
<div>set.seed(1234) #semente</div>
<div>sigma <- matrix(c<a href="tel:%2819260.71398" value="+551926071398" target="_blank">(19260.71398</a>,28.92129,28.92129,2.173789),
ncol=2)</div>
<div>sim_bnorm <- rmvnorm(n=1000, mean=c(mean(RST$RS),mean(RST$Temp)),
sigma=sigma)</div>
<div>bnorm <- data.frame(sim_bnorm)</div>
<div> </div>
<div>simulado <- sample(bnorm$X2,31)</div>
<div> </div>
<div>
<div>modelo <- lm(simulado~RST$Temp)</div>
<div>summary(modelo)</div>
<div> </div>
<div>plot( RST$Temp , simulado,xlab="Valores reais",ylab="Valores
simulados")</div>
<div>abline(modelo,col="red")</div></div>
<div> </div>
<div>Muito obrigada</div>-- <br><br style="FONT-FAMILY:tahoma;FONT-SIZE:13px"><font size="4" face="Times New Roman"><br></font>
<div> </div>
<div> </div></div>
</div></div><p>
</p><hr>
_______________________________________________<br>R-br mailing
list<br><a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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.<p></p></div></div></div></div>
<br>_______________________________________________<br>
R-br mailing list<br>
<a href="mailto:R-br@listas.c3sl.ufpr.br" target="_blank">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></blockquote></div><br><br clear="all"><div><br></div>-- <br>
<br style="font-family:Tahoma;font-size:13px"><font face="'Times New Roman'" size="4">Natália da Silva Martins<br>Bacharel em Estatística - Universidade Estadual de Maringá/ UEM<br>Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP<br>
Contato: <a href="tel:%2819%29%208306-4743" value="+551983064743" target="_blank">(19) 8306-4743</a></font> <br><div><br></div><br>