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) | ||
