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 Ambos lados da revisão seguinte | ||
|
pessoais:eder [2011/06/15 22:07] eder [section 5] |
pessoais:eder [2011/06/16 18:11] eder [section 5] |
||
|---|---|---|---|
| Linha 231: | Linha 231: | ||
| qbeta(c(0.025,0.975),alfa+sum(x),beta+sum(n-x)) | qbeta(c(0.025,0.975),alfa+sum(x),beta+sum(n-x)) | ||
| ###------------------------------------------------------------### | ###------------------------------------------------------------### | ||
| - | ###------------------------------------------------------------### | + | ##------------------------------------------------------------### |
| ###-----------------------------------------------------------------### | ###-----------------------------------------------------------------### | ||
| ### Regressão Beta | ### Regressão Beta | ||
| Linha 246: | Linha 246: | ||
| return(ll) | return(ll) | ||
| } | } | ||
| + | |||
| ###-----------------------------------------------------------------### | ###-----------------------------------------------------------------### | ||
| opt <- optim(c(B0=-0.5,B1=-0.51,B2=0.11,phi=35),log.vero,y=FoodExpenditure$food/FoodExpenditure$income, | opt <- optim(c(B0=-0.5,B1=-0.51,B2=0.11,phi=35),log.vero,y=FoodExpenditure$food/FoodExpenditure$income, | ||
| Linha 257: | Linha 257: | ||
| summary(fe_beta) | summary(fe_beta) | ||
| ###-----------------------------------------------------------------### | ###-----------------------------------------------------------------### | ||
| + | log.veroP <- function(par,phi,y,x1,x2){ | ||
| + | mu <- exp((par[1] + par[2] * x1 + par[3] * x2))/(1+exp((par[1] + par[2] * x1 + par[3] * x2)))##logit^-1 | ||
| + | ll <- sum(dbeta(y, mu* phi, (1-mu)*phi,log = TRUE)) | ||
| + | return(ll) | ||
| + | } | ||
| + | |||
| + | opt <- grid.phi <- seq(20,60,l=150) | ||
| + | con <- 1 | ||
| + | for (i in grid.phi){ | ||
| + | opt[con] <- optim(c(B0=-0.5,B1=-0.51,B2=0.11),log.veroP,phi=i,y=FoodExpenditure$food/FoodExpenditure$income, | ||
| + | x1=FoodExpenditure$income, | ||
| + | x2=FoodExpenditure$persons, | ||
| + | hessian = TRUE, control=(list(fnscale=-1)))$value | ||
| + | con <- con+1 | ||
| + | } | ||
| + | |||
| + | plot(grid.phi,2*(max(opt)-opt),type='l') | ||
| + | abline(h=3.84) | ||
| </code> | </code> | ||
