MANOVA в R - ошибка dimnames

Question

Я хочу использовать тест MANOVA в R, чтобы проверить статистическую разницу между вектором средних значений ("test") и вектором (равной длины) единиц ("random"). Мои данные как ниже:

    ID      openwater    closedshrubland      barren    cropnatural         crop   decidbroad     mixed
  test      0.8435707           1.037015     0.90612      0.8724474    0.9837284     1.035505   1.05059
random              1                  1           1              1            1            1         1

Когда я запускаю следующий код:

manovares1 <- manova(cbind(openwater,closedshrubland,barren,cropnatural,crop,decidbroad,mixed) ~ as.factor(ID),data=test)

Распечатка дает мне следующее сообщение:

Error in print.aov(list(coefficients = c(1, -0.156429334, 1, 0.0370145719999997,  : 
  length of 'dimnames' [1] not equal to array extent

И summary(manovares1) дает:

              Df
as.factor(ID)  1**

Есть идеи, что я делаю не так? (Подобный вопрос был задан здесь, но ответ, похоже, не помогает в моей ситуации.

Спасибо за любую помощь, которую вы можете предложить!

Answer 1

Я не уверен, что происходит с вашей базой данных, но держу пари, что это проблема.

Каковы ваши удостоверения личности, это открытая вода и т. Д.? В этом случае это то, что не так, потому что каждая группа в факторе имеет ровно один случай. Ни одна программа не даст вам лучших результатов при рассмотрении дисперсии. Так как нет ни одного.