Факторы в aov () - PullRequest
       5

Факторы в aov ()

1 голос
/ 27 апреля 2011

У меня возникла проблема с подключением по примеру 1 в R Guide . Вот пример

> datafilename="http://personality-project.org/r/datasets/R.appendix1.data"
> data.ex1 = read.table(datafilename,header=T)   #read the data into a table 
> aov.ex1 = aov(Alertness~Dosage,data=data.ex1)  #do the analysis of variance
> summary(aov.ex1)                                    #show the summary table

Но когда я применил aov к своим собственным данным, все изменилось.

> test.data <- data.frame(fac=letters[c(1:3,1:3)], x=1:6)
> test.result <- aov(fac~x, data=test.data)
Error in storage.mode(y) <- "double" :
  invalid to change the storage mode of a factor
In addition: Warning message:
In model.response(mf, "numeric") :
  using type="numeric" with a factor response will be ignored

Я полностью сбит с толку. В чем разница между test.data и data.ex1 в примере руководства по R?

> str(test.data)
'data.frame':   6 obs. of  2 variables:
 $ fac: Factor w/ 3 levels "a","b","c": 1 2 3 1 2 3
 $ x  : int  1 2 3 4 5 6
> str(data.ex1)
'data.frame':   18 obs. of  2 variables:
 $ Dosage   : Factor w/ 3 levels "a","b","c": 1 1 1 1 1 1 2 2 2 2 ...
 $ Alertness: int  30 38 35 41 27 24 32 26 31 29 ...

1 Ответ

4 голосов
/ 27 апреля 2011

должно быть aov(x ~ fac, data = test.data), что работает. Формула должна быть фактором отклика, а не фактором отклика.

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...