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 <[email protected]> > *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.
