Para o texto em itálico, considere o CMR abaixo.

library(lattice)

# Nomes das faixas em itálico.
xyplot(Sepal.Length ~ Petal.Length | Species,
       data = iris,
       strip = strip.custom(par.strip.text = list(font = 3)))

Sobre R², se realmente você quer apresentar um, pode considerar o quadrado da correlação entre observado e predito que justamente é o R² em modelos lineares. Esse não é único R² usado para GLM, tem uns como porcentagem de deviance explicada, enfim.

# Correlação entre observado e predito ao quadrado.
cor(iris$Sepal.Length, fitted(m0))^2
summary(m0)$r.squared

À disposição.

Walmes.

​