2016-07-31 26 views
10

R pakietu ggpmisc mogą służyć do wskazywania równanie lm modelu i poly modelu na ggplot2 (See here odniesienie). Zastanawiam się, jak wyświetlić wyniki równania modelu nls na ggplot2 przy użyciu ggmisc. Poniżej znajduje się mój MWE. Jakakolwiek pomoc będzie doceniona. Dziękirównanie Pokazano modelu NLS z ggpmisc

library(ggpmisc) 
args <- list(formula = y ~ k * e^x, 
      start = list(k = 1, e = 2)) 
ggplot(mtcars, aes(wt, mpg)) + 
    geom_point() + 
    stat_fit_augment(method = "nls", 
        method.args = args) 
+3

Użyj "adnotacji". –

Odpowiedz

2

Zainspirowany posta powiązane. Użyj geom_text, aby dodać etykietę po wyodrębnieniu parametrów.

nlsFit <- 
    nls(formula = mpg ~ k * e^wt, 
     start = list(k = 1, e = 2), 
     data = mtcars) 

nlsParams <- 
    nlsFit$m$getAllPars() 

nlsEqn <- 
    substitute(italic(y) == k %.% e^italic(x), 
      list(k = format(nlsParams['k'], digits = 4), 
        e = format(nlsParams['e'], digits = 2))) 

nlsTxt <- 
    as.character(as.expression(nlsEqn)) 

ggplot(mtcars, aes(wt, mpg)) + 
    geom_point() + 
    stat_fit_augment(method = "nls", 
        method.args = args) + 
    geom_text(x = 5, y = 30, label = nlsTxt, parse = TRUE)