Objekt nicht gefunden-Fehler in der Funktion

Ich habe eine Funktion, die bestimmt ist, um die Ausgabe der Zusammenfassung eines ANCOVA und Plotten Sie die Ergebnisse:

statAncova <- function (dataframe, response, covariate, Factor) {

  library(ggplot2)
  mod <- aov(response ~ covariate + Factor, data=dataframe)
  pred <- predict(mod)
  plotMod <- ggplot(data = cbind(mod$model, pred), aes(covariate, response, color=Factor)) + 
    geom_point() +
    facet_grid(. ~ Factor) + 
    geom_line(aes(y=pred))

  return(list(mod, plotMod))

}

Wenn ich versuche zu verwenden, Funktion wie diese:

statAncova(mtcars, drat, hp, cyl)

Bekomme ich diesen Fehler:

Error in eval(expr, envir, enclos) : object 'drat' not found

Was mache ich falsch?

  • Es gibt mehr fundamentale Fragen hier, vor allem, dass R nicht erkennen response, covariate oder Factor überhaupt. Sie müssen übergeben Sie als Zeichenfolgen, und bauen alle Formeln von hand.
InformationsquelleAutor luciano | 2013-06-25
  1. 3

    R erwartet für die es ein Objekt namens 'verflixt' in der Umwelt, aber 'verflixt' ist ein Mitglied der Datenrahmen mtcars.

    Ich weiß, es gibt wahrscheinlich eine elegantere Lösung, aber eine Möglichkeit, es zu beheben wäre:

    statAncova <- function (dataframe, response, covariate, Factor) {

  library(ggplot2)
  eval(parse(text=paste0("mod <- aov(",response," ~ ",covariate," + ",Factor,", data=dataframe)")))
  pred <- predict(mod)
  eval(parse(text=paste0("plotMod <- ggplot(data = cbind(mod$model, pred), aes(",covariate,", ",response,", color=",Factor,")) + 
  geom_point() +
  facet_grid(. ~ ",Factor,") + 
  geom_line(aes(y=pred))")))

  return(list(mod, plotMod))

}

statAncova(mtcars, "drat", "hp", "cyl")

    Alternativ könnte man einfach die einzelnen Variablen, die Sie interessiert sind in:

    statAncova <- function (response, covariate, Factor) {

  dataframe <- data.frame(response, covariate, Factor)

  library(ggplot2)
  mod <- aov(response ~ covariate + Factor, data=dataframe)
  pred <- predict(mod)
  plotMod <- ggplot(data = cbind(mod$model, pred), aes(covariate, response, color=Factor)) + 
    geom_point() +
    facet_grid(. ~ Factor) + 
    geom_line(aes(y=pred))

  return(list(mod, plotMod))

}

statAncova(mtcars$drat, mtcars$hp, mtcars$cyl)
    InformationsquelleAutor Sam Dickson
Schreibe einen Kommentar