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.