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

[Natalia Martins]

Muito Obrigada!

:)

Em 17 de junho de 2012 02:08, Sandro Sperandei <ssperandei em 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 em gmail.com>
> *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 <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 <%2819%29%208306-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.
>
>
> _______________________________________________
> 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/20120617/20fd2d39/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/20fd2d39/attachment.obj>