2レベル未満の場合の式からのドロップファクタ

Question

ファクタが2レベル未満かどうかに基づいて式からタームを削除したいと考えています。

私の設定では、データのサブセットごとに同じロジットモデルを多数実行したいと考えています。問題は、一部のサブセットでは、一部の要因では1つの値しか取れないことです。ここ

は一例であり：

depvar <- sample(0:1, 100, replace=T) 
bin1 <- rep(1, 100) 
cont1 <- rnorm(100) 
df <- data.frame(depvar=depvar, bin1=bin1, cont1=cont1) 
formla <- as.formula("depvar ~ bin1 + cont1") 
logit.reg <- glm(formla, data=df, family=binomial(link=logit)) 

私はそれが唯一の値をとるときに何とか自動的 式からbin1をドロップできるようにしたいです。

Answer 1

十分です/他の場合には、単純なように私には思える：BIN1はあなたの例では要因ではないことを

if(length(unique(df$bin)) > 1) { 
    formla <- as.formula("depvar ~ bin1 + cont1") 
} else { 
    formla <- as.formula("depvar ~ cont1") 
} 

注（1かかわらずにそれを強制するためにas.factor()を使用することができます）