Use

## sum(resid(modelo)^2): use deviance().
deviance(modelo)

## estimativa de sigma, use:
summary(modelo)$sigma

## R²
summary(modelo)$r.squared
summary(modelo)$adj.r.squared

Os demais você calcula a partir dos resíduos que estraí com residuals(modelo). Pode ser que o pacote car (companium of applied regression) tenha funções diretas para essas medidas.

À disposição.

Walmes.