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

[Natalia Martins]

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 <ssperandei em hotmail.com>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 <nsmbarreto em gmail.com>
> *Sent:* Saturday, June 16, 2012 7:39 PM
> *To:* R-br em 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 em 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 em 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.
>



-- 

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