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

[Natalia Martins]

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