From: Natalia Martins
Sent: Saturday, June 16, 2012 8:26 PM
To: [email protected]
Subject: Re: [R-br] simulação da normal bivariada
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.
Em 16 de junho de 2012 20:18, Sandro Sperandei <[email protected]>
escreveu:
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: [email protected]
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
[email protected]
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
[email protected]
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.
--
Natália da Silva Martins
Bacharel em Estatística - Universidade Estadual de Maringá/ UEM
Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP
Contato: (19) 8306-4743
--------------------------------------------------------------------------------
_______________________________________________
R-br mailing list
[email protected]
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
[email protected]
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.
