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

[Sandro Sperandei]

Assim sendo, Natalia, parece que seus dados não são muito correlacionados mesmo... Como disse, com base na sua matriz de covariância, o que vc pode esperar é uma correlação de 0,19 entre os dados. Se essa matriz foi gerada a partir de dados reais, é isso mesmo! É o que vc vai usar... C’est la vie! 

S



From: Natalia Martins 
Sent: Saturday, June 16, 2012 8:26 PM
To: r-br em listas.c3sl.ufpr.br 
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 <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 
  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 






--------------------------------------------------------------------------------
_______________________________________________
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.
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120617/406714d4/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: wlEmoticon-openmouthedsmile[1].png
Tipo: application/octet-stream
Tamanho: 1106 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20120617/406714d4/attachment.obj>