Natalia,

 

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?

 

S

 

 

 

From: Natalia Martins

Sent: Saturday, June 16, 2012 7:39 PM

To: R-br@listas.c3sl.ufpr.br

Subject: [R-br] simulação da normal bivariada

 

Boa noite pessoal,

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?

 

a programação que estou utilizando segue a baixo:

 

# gerando os dados normal bivariada

set.seed(1234) #semente

sigma <- matrix(c(19260.71398,28.92129,28.92129,2.173789), ncol=2)

sim_bnorm <- rmvnorm(n=1000, mean=c(mean(RST$RS),mean(RST$Temp)), sigma=sigma)

bnorm <- data.frame(sim_bnorm)

 

simulado <- sample(bnorm$X2,31)

 

modelo <- lm(simulado~RST$Temp)

summary(modelo)

 

plot( RST$Temp , simulado,xlab="Valores reais",ylab="Valores simulados")

abline(modelo,col="red")

 

Muito obrigada

--

 

 

---

_______________________________________________
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.

_______________________________________________
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.