Плутанина між рівнем фактора та міткою факторів

Здається, існує різниця між рівнями та мітками фактора в Р. Дотепер я завжди вважав, що рівні - це "справжня" назва рівнів факторів, а мітки - це імена, що використовуються для виведення (наприклад, таблиці та графіки) . Очевидно, що це не так, як показує наступний приклад:

df <- data.frame(v=c(1,2,3),f=c('a','b','c'))
str(df)
'data.frame':   3 obs. of  2 variables:
 $ v: num  1 2 3
 $ f: Factor w/ 3 levels "a","b","c": 1 2 3

df$f <- factor(df$f, levels=c('a','b','c'),
  labels=c('Treatment A: XYZ','Treatment B: YZX','Treatment C: ZYX'))
levels(df$f)
[1] "Treatment A: XYZ" "Treatment B: YZX" "Treatment C: ZYX"

Я думав, що рівні ('a', 'b', 'c') якимось чином можна отримати доступ до сценаріїв, але це не працює:

> df$f=='a'
[1] FALSE FALSE FALSE

Але це робить:

> df$f=='Treatment A: XYZ' 
[1]  TRUE FALSE FALSE

Отже, моє запитання складається з двох частин:

• Яка різниця між рівнями та мітками?

• Чи можливо мати різні назви рівнів факторів для сценаріїв та результатів?

Передумови: Для більш довгих сценаріїв сценарій з короткими рівнями фактора здається набагато простішим. Однак для звітів та сюжетів цей короткий коефіцієнт може бути недостатнім, і його слід замінити точнішими назвами.

Дуже коротко: рівні - це вхід, мітки - вихід у factor()функції. Фактор має лише levelатрибут, який задається labelsаргументом у factor()функції. Це відрізняється від концепції міток у статистичних пакетах, таких як SPSS, і на початку може бути заплутаним.

Що ви робите в цьому рядку коду

df$f <- factor(df$f, levels=c('a','b','c'),
  labels=c('Treatment A: XYZ','Treatment B: YZX','Treatment C: ZYX'))

говорить R, що існує вектор df$f

• який ви хочете перетворити на фактор,
• в яких різні рівні кодуються як a, b і c
• і для яких ви хочете, щоб рівні були позначені як Лікування А тощо.

Функція фактора шукатиме значення a, b і c, перетворює їх у числові класи фактора та додає значення мітки до levelатрибута фактора. Цей атрибут використовується для перетворення внутрішніх числових значень у правильні мітки. Але, як бачите, labelатрибута немає .

> df <- data.frame(v=c(1,2,3),f=c('a','b','c'))    
> attributes(df$f)
$levels
[1] "a" "b" "c"

$class
[1] "factor"

> df$f <- factor(df$f, levels=c('a','b','c'),
+   labels=c('Treatment A: XYZ','Treatment B: YZX','Treatment C: ZYX'))    
> attributes(df$f)
$levels
[1] "Treatment A: XYZ" "Treatment B: YZX" "Treatment C: ZYX"

$class
[1] "factor"

Я написав пакет "lfactors", який дозволяє посилатися на рівні або мітки.

# packages
install.packages("lfactors")
require(lfactors)

flips <- lfactor(c(0,1,1,0,0,1), levels=0:1, labels=c("Tails", "Heads"))
# Tails can now be referred to as, "Tails" or 0
# These two lines return the same result
flips == "Tails"
#[1]  TRUE FALSE FALSE  TRUE  TRUE FALSE
flips == 0 
#[1]  TRUE FALSE FALSE  TRUE  TRUE FALSE

Зауважте, що lfactor вимагає, щоб рівні були числовими, щоб їх не можна переплутати з мітками.