的R中一個定義模型循環

Question

我檢查this問題有關定義在R車型，我想什麼做的是在一個型號爲循環交換變量，這樣每一個變量是一個時間目標變量而所有其他變量都是該迭代的預測變量。

df <- data.frame(customer = c(1:5), product1 = c(1,0,1,1,0), product2 = c(0,1,0,0,1), product3 = c(0,1,1,1,0)) 

customer product1 product2 product3 
1  1  1  0  0 
2  2  0  1  1 
3  3  1  0  1 
4  4  1  0  1 
5  5  0  1  0 

於是我想創建一個for循環使用3次迭代在這種情況下：

mdl <- product1 ~ product2 + product3 
mdl <- product2 ~ product1 + product3 
mdl <- product3 ~ product1 + product2 

要在這裏澄清，我的問題我試圖創建for循環的：

for(j in 1:ncol(df)){ 
    mdl <- df$[j] ~ df[-j] # include all variables except target variable 
    print(mdl) 
    } 

這裏我得到的輸出：

df[j] ~ df[-j] 
df[j] ~ df[-j] 
df[j] ~ df[-j] 
df[j] ~ df[-j] 

雖然我預計一些諸如所需的輸出：

product1 ~ product2 + product3 
product2 ~ product1 + product3 
product3 ~ product1 + product2 

如果你想知道爲什麼我想知道這一點。我想用它在一個for循環，運行的預測模型，如下例所示：

naiveBayes(mdl, df, type = "raw") 

我希望我的問題是清楚的，希望有人可以幫助我。

Answer 1

隨着setdiff和lapply你可以實現公式的組合。

varNames = colnames(DF)[-1] 


lapply(varNames, function(x) paste0(x ,"~", paste0(setdiff(varNames,x),collapse="+"))) 
#[[1]] 
#[1] "product1~product2+product3" 
# 
#[[2]] 
#[1] "product2~product1+product3" 
# 
#[[3]] 
#[1] "product3~product1+product2" 

把這些變成你的模型，你可以這樣做：

modelList = lapply(varNames,function(x) { 

depVar = x 
indepVar = setdiff(varNames,x) 

formulaVar = as.formula(paste0(depVar ,"~", paste0(indepVar,collapse="+"))) 

nbModel = naiveBayes(formulaVar, df, type = "raw") 

outputList = list(indepVar = paste0(indepVar,collapse=","),depVar = depVar,nbModel = nbModel) 

return(outputList) 

}) 

這將返回一個包含相關的增值經銷商，獨立var和最後樸素貝葉斯模型的列表對象。

要訪問其中的任何一個，您可以modelList[[1]],length(modelList)給出此列表中的模型數量。

對於替代方法來生成組合，見?combn和?expand.grid

combn(varNames,2) 
#  [,1]  [,2]  [,3]  
#[1,] "product1" "product1" "product2" 
#[2,] "product2" "product3" "product3" 

head(expand.grid(varNames,varNames,varNames)) 
#  Var1  Var2  Var3 
#1 product1 product1 product1 
#2 product2 product1 product1 
#3 product3 product1 product1 
#4 product1 product2 product1 
#5 product2 product2 product1 
#6 product3 product2 product1