Não foi possível enviar o arquivo. Será algum problema com as permissões?
Diferenças
Aqui você vê as diferenças entre duas revisões dessa página.
Ambos lados da revisão anterior Revisão anterior Próxima revisão | Revisão anterior Próxima revisão Ambos lados da revisão seguinte | ||
pessoais:eder [2011/09/25 20:50] eder [Códigos] |
pessoais:eder [2011/09/26 08:01] eder [section 5] |
||
---|---|---|---|
Linha 14: | Linha 14: | ||
* {{:pessoais:inlarrblup.r|GWS}} Seleção Genômica Ampla Via ML REML INLA | * {{:pessoais:inlarrblup.r|GWS}} Seleção Genômica Ampla Via ML REML INLA | ||
* {{:pessoais:reml_inla.r|Script}} Modelo seleção Genótipo ambiente via REML ML INLA | * {{:pessoais:reml_inla.r|Script}} Modelo seleção Genótipo ambiente via REML ML INLA | ||
+ | * {{:pessoais:linearregression.rnw|Script}} Regressão Linear - inferência via Mínimos quadrados, ML, REML, Gibbs, Metropolis, INLA, dclone ... (Em construção) | ||
+ | |||
===== Disciplinas 2011/1 ===== | ===== Disciplinas 2011/1 ===== | ||
* [[http://www.leg.ufpr.br/doku.php/disciplinas:ce210-2010-02|CE-210: Inferência estatística II]] | * [[http://www.leg.ufpr.br/doku.php/disciplinas:ce210-2010-02|CE-210: Inferência estatística II]] | ||
Linha 61: | Linha 63: | ||
} | } | ||
num <- (sum(dnorm(y,b0_n+b1_n*x,sigma_n,log=TRUE))#+ | num <- (sum(dnorm(y,b0_n+b1_n*x,sigma_n,log=TRUE))#+ | ||
- | # sum(dnorm(y,mu0[1],V0[1,1],log=TRUE))+ | + | #sum(dnorm(b0_n,mu0[1],V0[1,1],log=TRUE))+ |
- | #sum(dnorm(y,mu0[1],V0[2,2],log=TRUE))+ | + | #sum(dnorm(b1_n,mu0[1],V0[2,2],log=TRUE))+ |
- | # sum(log(dinvgamma(y,v0,tau0))) | + | #sum(log(dinvgamma(sigma_n,v0,tau0))) |
) * u^4 | ) * u^4 | ||
den <- (sum(dnorm(y,b01_n+b11_n*x+b21_n*x^2,sigma1_n,log=TRUE))#+ | den <- (sum(dnorm(y,b01_n+b11_n*x+b21_n*x^2,sigma1_n,log=TRUE))#+ | ||
- | # sum(dnorm(y,mu02[1],V02[1,1],log=TRUE))+ | + | # sum(dnorm(b01_n,mu02[1],V02[1,1],log=TRUE))+ |
- | # sum(dnorm(y,mu02[2],V02[2,2],log=TRUE))+ | + | # sum(dnorm(b11_n,mu02[2],V02[2,2],log=TRUE))+ |
- | # sum(dnorm(y,mu02[3],V02[3,3],log=TRUE))+ | + | # sum(dnorm(b21_n,mu02[3],V02[3,3],log=TRUE))+ |
- | # sum(log(dinvgamma(y,v02,tau02))) | + | # sum(log(dinvgamma(sigma1_n,v02,tau02))) |
) * dnorm(u,0,2) | ) * dnorm(u,0,2) | ||
u = runif(1, 0, 1) | u = runif(1, 0, 1) |