Тест Левена в R для приведенной ниже таблицы данных

Question

dput(head(Clean, 100))

structure(list(`Reporting unit` = c("Albury Wodonga Health [Albury Campus]", 
"Albury Wodonga Health [Albury Campus]", "Albury Wodonga Health [Albury Campus]", 
"Albury Wodonga Health [Albury Campus]", "Albury Wodonga Health [Albury Campus]", 
"Albury Wodonga Health [Albury Campus]", "Albury Wodonga Health [Albury Campus]", 
"Albury Wodonga Health [Albury Campus]", "Albury Wodonga Health [Albury Campus]", 
"Albury Wodonga Health [Albury Campus]"), `Reporting unit type` = c("Hospital", 
"Hospital", "Hospital", "Hospital", "Hospital", "Hospital", "Hospital", 
"Hospital", "Hospital", "Hospital"), State = c("NSW", "NSW", 
"NSW", "NSW", "NSW", "NSW", "NSW", "NSW", "NSW", "NSW"), `Local Hospital Network (LHN)` = c("Albury Wodonga Health", 
"Albury Wodonga Health", "Albury Wodonga Health", "Albury Wodonga Health", 
"Albury Wodonga Health", "Albury Wodonga Health", "Albury Wodonga Health", 
"Albury Wodonga Health", "Albury Wodonga Health", "Albury Wodonga Health"
), `Peer group` = c("Large hospitals", "Large hospitals", "Large hospitals", 
"Large hospitals", "Large hospitals", "Large hospitals", "Large hospitals", 
"Large hospitals", "Large hospitals", "Large hospitals"), `Time period` = c("2011–12", 
"2012–13", "2013–14", "2014–15", "2015–16", "2016–17", "2011–12", 
"2012–13", "2013–14", "2014–15"), Category = c("Cellulitis", 
"Cellulitis", "Cellulitis", "Cellulitis", "Cellulitis", "Cellulitis", 
"Chronic Obstructive Pulmonary Disease (without complications)", 
"Chronic Obstructive Pulmonary Disease (without complications)", 
"Chronic Obstructive Pulmonary Disease (without complications)", 
"Chronic Obstructive Pulmonary Disease (without complications)"
), `Total number of stays` = c(111, 116, 141, 155, 210, 196, 
109, 116, 75, 132), `Number of overnight stays` = c(92, 98, 115, 
123, 166, 155, 108, 113, 71, 122), `Percentage of overnight stays` = c(0.83, 
0.84, 0.82, 0.79, 0.79, 0.79, 0.99, 0.97, 0.95, 0.92), `Average length of stay (days)` = c(3.9, 
3.3, 3.1, 2.5, 2.6, 2.7, 5.8, 4.6, 5.7, 4.4), `Peer group average (days)` = c(3.7, 
3.5, 3.3, 3.2, 3, 3, 4.8, 4.4, 4.2, 3.9), `Total overnight patient bed days` = c(356, 
326, 351, 306, 431, 418, 622, 518, 405, 538)), row.names = c(NA, 
-10L), class = c("tbl_df", "tbl", "data.frame"))

Здравствуйте, я пытаюсь выполнить тест levene для указанного выше фрейма данных между равноправной группой (крупные и средние больницы) и средней длиной, но получаю следующую ошибку:

Ошибка в leveneTest.formula (Чистая $ group ~ Чистая $ Average length of stay (days),: тест Левена не подходит для количественных объяснительных переменных.

Может кто-нибудь помочь, пожалуйста .

Answer 1

Решение может решить проблему.

1. Если используется интерфейс формулы, в нем не нужно указывать имя data.frame, только в аргументе data.
2. Порядок терминов формулы должен быть обратным, объясняющая переменная должна быть в RHS.

И поскольку пакет, в который leveneTest может не загружаться, я добавляю префикс вызова функции.

car::leveneTest(`Average length of stay (days)` ~ `Peer group`, data = Clean)

Кроме того, имена столбцов синтаксически недопустимы, пробелы в именах нуждаются в обратном тикании.

names(Clean) <- make.names(names(Clean))

Answer 2

Обратный порядок в вашей формуле:

leveneTest(Clean$`Average length of stay (days)`~Clean$`Peer group`)