Übergeben Sie einen Vektor, der die Variablen in die lm () - Formel

War ich versucht, zu automatisieren, ein Stück von meinem code, so dass die Programmierung weniger aufwendig.

Im Grunde war ich versucht zu tun, eine schrittweise Selektion der Variablen mit fastbw() im rms Paket. Ich würde gerne weitergeben der Liste der Variablen ausgewählt, von fastbw() in eine Formel als y ~ x1+x2+x3, "x1" "x2" "x3" wird die Liste der Variablen ausgewählt, von fastbw()

Hier ist der code, den ich versucht und hat nicht funktioniert

olsOAW0.r060 <- ols(roll_pct~byoy+trans_YoY+change18m, 
                    subset= helper=="POPNOAW0_r060", 
                    na.action = na.exclude, 
                    data = modelready)

OAW0 <- fastbw(olsOAW0.r060, rule="p", type="residual", sls= 0.05)

vec <- as.vector(OAW0$names.kept, mode="any")

b <- paste(vec, sep ="+") ##I even tried b <- paste(OAW0$names.kept, sep="+")

bestp.OAW0.r060 <- lm(roll_pct ~ b , 
                      data = modelready, 
                      subset = helper =="POPNOAW0_r060",    
                      na.action = na.exclude)

Ich bin neu in R und noch nicht gezogenen steilen Lernkurve, so entschuldigen für offensichtliche Programmier-Fehler.

InformationsquelleAutor Anand | 2012-02-11
  1. 56

    Du bist fast da. Sie müssen nur paste die gesamte Formel zusammen, so etwas wie dieses:

    paste("roll_pct ~ ",b,sep = "")

    zwingen, es zu einer tatsächlichen Formel mit as.formula und dann pass , dass zu lm. Technisch, denke ich lm möglicherweise zwingen eine Zeichenkette selbst, aber die Nötigung selbst ist in der Regel sicherer. (Einige Funktionen, die erwarten, dass die Formeln nicht aus dem Zwang für Sie, werden andere.)

    • JJoran, danke, ich habe immer noch Fehler .bitte werfen Sie einen Blick hpi < - "einfügen" ("byoy", "change18m" , "change24m" , "change18m0", "change24m0", "lag3byoy" , "lag3change18m" , "lag3change24m", "lag6byoy", "lag6change18m", "lag6change24m", "trans_YoY" , sep ="+") hpi.form <- wie.formula(paste("roll_pct~", "hpi", sep = "")) lmNGC0.r060 <- lm(hpi.form, subset= Helfer== "POPNOANGC0_r060",na.action = na.ausschließen, data = modelready) > Fehler im Modell.frame.default(formula = als.formula(paste("roll_pct~", "hpi",:variable lengths differ (found for 'hpi')
    • Sie setzen hpi in Anführungszeichen in der zweiten Zeile. Geben Sie hpi.form an die Konsole und du wirst sehen, warum diese falsch ist.
    • Joran, noch einmal vielen Dank, ich denke, dass lm() ist nicht nötigt das Objekt "hpi" in der Formel als Zeichenkette. Als ich versuchte, die Eingabe der Variablen in die lm () - Formel, es arbeitete
    • Nein, Sie klebte es zusammen falsch. Sie schrieb: paste("roll_pct~", "hpi", sep = "") eher als paste("roll_pct~", hpi, sep = "") wie ich bereits erwähnt habe in meiner Antwort.
    • Joran, vielen Dank, Sie es bemerkte . Habe es diesmal funktionieren..... danke für die Geduld mit mir
    • funktioniert nicht mit Formel in Werte::Aggregat-Funktion
    • Ja, aggregate ist einer der Fälle, wo Sie brauchen, um zu zwingen es, sich über as.formula.

    InformationsquelleAutor joran
  2. 18

    Würden Sie wirklich brauchen, zu verwenden Zusammenbruch statt der seb bei der Definition von b ist. 

    b <- paste(OAW0$names.kept, collapse="+")

    Dann können Sie es in joran Antwort

    paste("roll_pct ~ ",b,sep = "")

    oder nutzen Sie einfach:

    paste("roll_pct ~ ",paste(OAW0$names.kept, collapse="+"),sep = "")
    InformationsquelleAutor cconnell
  3. 8

    Ich lief in ähnliche Frage heute, wenn Sie wollen, um es noch allgemeiner, wo Sie noch nicht behoben haben Klasse Namen, die Sie verwenden können, 

    frmla <- as.formula(paste(colnames(modelready)[1], paste(colnames(modelready)[2:ncol(modelready)], sep = "", 
                              collapse = " + "), sep = " ~ "))

    Davon ausgegangen, dass Sie Klasse haben, Variablen oder der abhängigen Variablen in der ersten Spalte aber die Indizierung kann einfach gewechselt werden, um Letzte Spalte:

    frmla <- as.formula(paste(colnames(modelready)[ncol(modelready)], paste(colnames(modelready)[1:(ncol(modelready)-1)], sep = "", 
                              collapse = " + "), sep = " ~ "))

    Dann weiter mit lm mit:

    bestp.OAW0.r060 <- lm(frmla , data = modelready, ... )
    • Diese Antwort ist fast 3 Jahre alt, aber sehr einfach und elegant ist und mir viel Zeit erspart - upvote.
    InformationsquelleAutor discipulus
  4. 1

    Wenn Sie auf der Suche nach etwas weniger ausführlich:

    fm <- as.formula( paste( colnames(df)[i], ".", sep=" ~ ")) 
                                      # i is the index of the outcome column

    Hier ist es in einer Funktion:

    getFormula<-function(target, df) {

  i <- grep(target,colnames(df))
  as.formula(paste(colnames(df)[i], 
                   ".", 
                   sep = " ~ "))
}
fm <- getFormula("myOutcomeColumnName", myDataFrame)
rp <- rpart(fm, data = myDataFrame) # Use the formula to build a model
    InformationsquelleAutor Travis Heeter
  5. 0

    nur zu vereinfachen und zu sammeln über Antworten, basierend auf einer Funktion

    my_formula<- function(colPosition, trainSet){
    dep_part<- paste(colnames(trainSet)[colPosition],"~",sep=" ")
    ind_part<- paste(colnames(trainSet)[-colPosition],collapse=" + ")
    dt_formula<- as.formula(paste(dep_part,ind_part,sep=" "))
    return(dt_formula)
}

    Verwendet werden:

    my_formula( dependent_var_position, myTrainSet)
    InformationsquelleAutor user3169330
Schreibe einen Kommentar