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 --
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@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@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.
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@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@gmail.com> *Sent:* Saturday, June 16, 2012 7:39 PM *To:* R-br@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@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@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
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@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@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@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@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@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@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.
Muito Obrigada! :) Em 17 de junho de 2012 02:08, Sandro Sperandei <ssperandei@hotmail.com>escreveu:
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! [image: Smiley de boca aberta]
S
*From:* Natalia Martins <nsmbarreto@gmail.com> *Sent:* Saturday, June 16, 2012 8:26 PM *To:* r-br@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@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@gmail.com> *Sent:* Saturday, June 16, 2012 7:39 PM *To:* R-br@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@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@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 <%2819%29%208306-4743>
------------------------------ _______________________________________________ R-br mailing list R-br@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@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
Olha só quem posta na lista... que bom que resolveu seu problema! 2012/6/17 Natalia Martins <nsmbarreto@gmail.com>
Muito Obrigada!
:)
Em 17 de junho de 2012 02:08, Sandro Sperandei <ssperandei@hotmail.com>escreveu:
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! [image: Smiley de boca aberta]
S
*From:* Natalia Martins <nsmbarreto@gmail.com> *Sent:* Saturday, June 16, 2012 8:26 PM *To:* r-br@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@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@gmail.com> *Sent:* Saturday, June 16, 2012 7:39 PM *To:* R-br@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@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@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 <%2819%29%208306-4743>
------------------------------ _______________________________________________ R-br mailing list R-br@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@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@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.
From: Natalia Martins Sent: Saturday, June 16, 2012 8:26 PM To: r-br@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@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@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@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@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@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.
Teu R² da regressão de Y1 em Y2 não absolutamente nenhuma relação com estar simulando certo ou errado. Se a correlação entre elas for baixa, menor R², só isso. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
Se vc aumentar os valores de 28 para 198, por exemplo, vai conseguir correlações próximas a 1. Cuidado também na hora da amostragem S From: Natalia Martins Sent: Saturday, June 16, 2012 7:39 PM To: R-br@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@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.
participantes (4)
-
FHRB Toledo -
Natalia Martins -
Sandro Sperandei -
Walmes Zeviani