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