A random effect formulation for R package nlme is read in as a string of characters from an input file: ranEff <- "pdCompSymm(~1+Age)" I need to convert 'ranEff' to a formula class. However, as shown below: > as.formula(ranEff) ~1 + Age the "pdCompSymm" is lost in the conversion. Any solutions? Thanks! Gang