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

[FHRB Toledo]

Olha só quem posta na lista... que bom que resolveu seu problema!

2012/6/17 Natalia Martins <nsmbarreto em gmail.com>

> 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
>
>
>
> _______________________________________________
> 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/b8ad3000/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/b8ad3000/attachment.obj>