X <- data[, 15:855]
    Y <- data[, 1:14]
    Y.matrix <- data.matrix(Y)
    X.matrix <- data.matrix(X)
    summary(manova(Y.matrix ~ X.matrix))

No error degrees of freedom
           Df
  X.matrix 94