<HTML><HEAD></HEAD>
<BODY dir=ltr>
<DIV dir=ltr>
<DIV style="FONT-FAMILY: 'Calibri'; COLOR: #000000; FONT-SIZE: 12pt">
<DIV>Se vc aumentar os valores de 28 para 198, por exemplo, vai conseguir
correlações próximas a 1. Cuidado também na hora da amostragem</DIV>
<DIV><BR>S</DIV>
<DIV> </DIV>
<DIV
style="FONT-STYLE: normal; DISPLAY: inline; FONT-FAMILY: 'Calibri'; COLOR: #000000; FONT-SIZE: small; FONT-WEIGHT: normal; TEXT-DECORATION: none">
<DIV style="FONT: 10pt tahoma">
<DIV> </DIV>
<DIV style="BACKGROUND: #f5f5f5">
<DIV style="font-color: black"><B>From:</B> <A title=nsmbarreto@gmail.com
href="mailto:nsmbarreto@gmail.com">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">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-STYLE: normal; DISPLAY: inline; FONT-FAMILY: 'Calibri'; COLOR: #000000; FONT-SIZE: small; FONT-WEIGHT: normal; TEXT-DECORATION: none">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" target=_blank
value="+551926071398">(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>
<P>
<HR>
_______________________________________________<BR>R-br mailing
list<BR>R-br@listas.c3sl.ufpr.br<BR>https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br<BR>Leia
o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo
reproduzível.</DIV></DIV></DIV></BODY></HTML>