[R-br] Modelos Mistos

Prezados, Estou utilizando da função lme (pacote nlme) para ajustar um modelo misto. Gostaria de saber porém se é possível flexibilizar a estrutura de correlação dos efeitos aleatórios para modelos fora da função "pdMat" e suas respectivas classes ( "pdClass"). Posso (vide código a seguir) propor uma estrutura bloco diagonal, com blocos na classe "corClasses" ou outras formas e/ou funções que eu mesmo descreveria? lme(Y~X1+X2-1, data=dados.frame, random=list(subject=pdBlocked(list(corARMA (~Z1+Z2-1), corARMA(~ Z3+Z4-1))))) Grato David