Pass-Formel, um eine Funktion in R?

Jede Hilfe bei diesem würde sehr geschätzt. Ich bin mit dem Lumley-Umfrage Paket-und versuche zu vereinfachen mein code, aber getroffen haben eine leichte Haken.

Den svymean Funktion aus dem Paket wie folgt aufgerufen wird, in meinem code, wobei das erste argument ist eine Formel, die angibt, welche Variablen ich will, und das zweite argument ist, dass dataset:

svymean(~hq_ehla, FraSvy, na.rm=TRUE)

Ich versuche, eine Funktion erstellen, die das herausziehen der Mittelwert (Anteile) und Standardfehler für kategoriale Variablen, also habe ich die folgende Funktion:

stats <- function(repstat, num) {
    estmean <- as.numeric(round(100 * repstat[num], digits=0))
    estse <- round(100 * sqrt(attributes(repstat)$var[num,num]), digits=1)
    return(list(mean=estmean, se=estse))
}

Das funktioniert, also wenn ich das herausziehen der Mittelwert und se von meinem ersten Kategorie, verwende ich zum Beispiel:

stats(svymean(~hq_ehla, FraSvy, na.rm=TRUE), 1)$mean
stats(svymean(~hq_ehla, FraSvy, na.rm=TRUE), 1)$se

Was ich möchte in der Lage sein zu tun ist, vereinfacht dies etwas viel kürzer, wo würde ich vielleicht nur noch schreiben:

stats(FraSvy, "hq_ehla", 1)$mean

Oder so ähnlich. Problem ist, dass ich kann nicht herausfinden, wie übergeben Sie eine Formel auf eine Funktion mit einer Variablen-Namen.

InformationsquelleAutor RickyB | 2013-02-03
  1. 12

    Können Sie reformulate zu konstruieren, die Formel und rufen svymean innerhalb Ihrer Funktion. Verwenden ... übergeben na.rm oder anderen Argumente zu svymean

    stats <- function(terms, data,  num, ...) {
  .formula <- reformulate(terms)
  repstat <- svymean(.formula, data, ...)
  estmean <- as.numeric(round(100 * repstat[num], digits=0))
  estse <- round(100 * sqrt(attributes(repstat)$var[num,num]), digits=1)
  return(list(mean=estmean, se=estse))
}

stats(data = FraSvy, terms = "hq_ehla", 1, na.rm = TRUE)$mean

    Haben Sie einen Blick auf diese Antwort für mehr details auf programmitically erstellen Formel-Objekte

    Oder Sie passieren konnten, ein Formel-Objekt innerhalb der Funktion.

    stats2 <- function(formula, data,  num, ...) {

  repstat <- svymean(formula, data, ...)
  estmean <- as.numeric(round(100 * repstat[num], digits=0))
  estse <- round(100 * sqrt(attributes(repstat)$var[num,num]), digits=1)
  return(list(mean=estmean, se=estse))
}


stats2(data = FraSvy, formula = ~hq_ehla, 1, na.rm = TRUE)$mean
    • Perfekt...keine Ahnung warum ich nicht finden konnte, dass, Dank!
    InformationsquelleAutor mnel
  2. 0

    den coef und SE Funktionen machen Ihr Leben leichter..

    # construct a function that takes the equation part of svymean as a string
# instead of as a formula.  everything else gets passed in the same
# as seen by the `...`
fun <- function( var , ... ) svymean( reformulate( var ) , ... )

# test it out.
result <- fun( "hq_ehla" , FraSvy , na.rm = TRUE )

# print the results to the screen
result

# also your components
coef( result )
SE( result )

# and round it
round( 100 * coef( result ) )
round( 100 * SE( result ) )
    InformationsquelleAutor Anthony Damico
Schreibe einen Kommentar