Я получаю объект ошибки типа "встроенный" не является поднабором при репликации кода

Question

Он написан Оортом и Джаком по вопросу об освобождении от МЛ с использованием OpenMx в r и опубликован в Research Synthesis Methods. Часть кода с ошибкой приведена здесь.

> fit <- mxFitFunctionML()
>  for (i in 1:length(cormat)){
+   data[[i]] <- mxData(observed = cormat[[i]], type = 'cov', numObs = samplesizes[[i]])   
+ matrixR[[i]] <-    mxMatrix(   type = 'Symm',   nrow = nrow(cormat[[i]]),   ncol = ncol(cormat[[i]]),   free = TRUE,   values = vech(startvalues[[i]]),
+   labels = vech(cormatlabels[[i]]),   name = 'R')
+   matrixD[[i]] <- mxMatrix(   type = 'Diag',   nrow = nrow(cormat[[i]]),   ncol = ncol(cormat[[i]]),   free = FALSE,   values = 1,   byrow = TRUE,   name = 'D')
+    matrixI <- mxMatrix(          type = "Iden",   nrow = nrow(cormat[[i]]),   ncol = ncol(cormat[[i]]),          name = "I")
+   exp[[i]] <- mxExpectationNormal(covariance="S",dimnames = rownames(cormat[[i]]))
+    model[[i]] <- mxModel(paste('sample',i,sep="") , data[[i]], matrixR[[i]],           matrixI,matrixD[[i]], algebraS, exp[[i]], fit)  }
Error in exp[[i]] <- mxExpectationNormal(covariance = "S", dimnames = rownames(cormat[[i]])) : 
  object of type 'builtin' is not subsettable

Я провел его в прошлом году без проблем! GunnarG