R'yi bir regresyonda referans olarak belirli bir faktör seviyesini kullanmaya nasıl zorlayabilirim?

Bir regresyonda ikili açıklayıcı değişkenler kullanırsam, R'ye belirli bir seviyeyi referans olarak kullanmasını nasıl söyleyebilirim?

Varsayılan olarak sadece bir seviye kullanıyor.

lm(x ~ y + as.factor(b)) 

ile b {0, 1, 2, 3, 4}. Diyelim ki R tarafından kullanılan sıfır yerine 3 kullanmak istiyorum.

İşleve bakın relevel(). İşte bir örnek:

set.seed(123)
x <- rnorm(100)
DF <- data.frame(x = x,
                 y = 4 + (1.5*x) + rnorm(100, sd = 2),
                 b = gl(5, 20))
head(DF)
str(DF)

m1 <- lm(y ~ x + b, data = DF)
summary(m1)

Şimdi faktörünü değiştirmek biçinde DFkullanımıyla relevel()fonksiyonu:

DF <- within(DF, b <- relevel(b, ref = 3))
m2 <- lm(y ~ x + b, data = DF)
summary(m2)

Modeller farklı referans seviyelerini tahmin etmiştir.

> coef(m1)
(Intercept)           x          b2          b3          b4          b5 
  3.2903239   1.4358520   0.6296896   0.3698343   1.0357633   0.4666219 
> coef(m2)
(Intercept)           x          b1          b2          b4          b5 
 3.66015826  1.43585196 -0.36983433  0.25985529  0.66592898  0.09678759

Diğerleri relevel, verilerinizdeki tüm analizler için temel düzeyi değiştirmek istiyorsanız (veya verileri değiştirerek yaşamaya istekliyseniz) en iyi çözüm olan komuttan bahsetmişlerdir .

Verileri değiştirmek istemiyorsanız (bu tek seferlik bir değişikliktir, ancak gelecekte varsayılan davranışı tekrar istiyorsanız), Ckontrastları ayarlamak için (büyük harf notu) işlevinin bir kombinasyonunu ve ile contr.treatmentsişlevi kullanabilirsiniz. temel olmasını istediğiniz seviyeyi seçmek için temel bağımsız değişken.

Örneğin:

lm( Sepal.Width ~ C(Species,contr.treatment(3, base=2)), data=iris )

relevel()Komut sorunuza kısayol yöntemidir. Yaptığı şey faktörü yeniden sıralamaktır, böylece ref seviyesi ne olursa olsun önce olur. Bu nedenle, faktör seviyelerinizi yeniden düzenlemek de aynı etkiye sahip olacak, ancak size daha fazla kontrol sağlayacaktır. Belki 3,4,0,1,2 seviyelerine sahip olmak istediniz. Bu durumda...

bFactor <- factor(b, levels = c(3,4,0,1,2))

Bu yöntemi tercih ediyorum çünkü kodumda yalnızca referansın ne olduğunu değil, diğer değerlerin konumunu da (bunun sonuçlarına bakmak yerine) görmek benim için daha kolay.

NOT: Bunu sıralı bir faktör YAPMAYIN. Belirli bir sıraya sahip bir faktör ve sıralı bir faktör aynı şey değildir. lm()Bunu yaparsanız polinom zıtlıkları istediğinizi düşünmeye başlayabilirsiniz.

Bunun eski bir soru olduğunu biliyorum, ancak benzer bir sorunla karşılaştım ve şunu buldum:

lm(x ~ y + relevel(b, ref = "3")) 

tam olarak istediğini yapıyor.

Ayrıca contrasts, regresyon işlevleri tarafından saygı duyulduğu görülen bir öznitelikle sütunu manuel olarak da etiketleyebilirsiniz :

contrasts(df$factorcol) <- contr.treatment(levels(df$factorcol),
   base=which(levels(df$factorcol) == 'RefLevel'))

Dplyr / tidyverse versiyon arayanlar için. Gavin Simpson çözümüne dayalı olarak:

# Create DF
set.seed(123)
x <- rnorm(100)
DF <- data.frame(x = x,
                 y = 4 + (1.5*x) + rnorm(100, sd = 2),
                 b = gl(5, 20))

# Change reference level
DF = DF %>% mutate(b = relevel(b, 3))

m2 <- lm(y ~ x + b, data = DF)
summary(m2)