# How do i predict with standard errors using betareg package in R?

I'm using 'betareg' package in R to perform beta regression. predict() function with se.fit=T is supposed to return standard errors along with the prediction but it doesn't. Is there any other way I can get the standard error outputs?

I'm open to using other packages that can perform beta regression too.

## 2 Answers

I've thought about two alternatives to get the prediction interval of the mean of the response:

I have used the package 'effects' and the function allEffects()

BetaReg <- betareg(Value~x, data = Data)
Effects <- as.data.frame(allEffects(BetaReg, xlevels=list(x=BetaReg$$x)))$$x
Effects


I think it is correct (looks about right on a graph) :)