Logarithmische Regression

Modelle zur Korrelations- und Regressionsanalyse

Moderator: EDi

Werekorden
Beiträge: 83
Registriert: So Feb 04, 2018 7:52 pm

Re: Logarithmische Regression

Beitrag von Werekorden »

Hallo,

ARRRGGHHHHH Einfach alles auf die neusten Versionen updaten hat geholfen. Es gibt keine Fehlermeldungen mehr.

leibt trotzdem noch das Formelproblem.

Danke bis hierhin.

VG
Andreas
jogo
Beiträge: 2085
Registriert: Fr Okt 07, 2016 8:25 am

Re: Logarithmische Regression

Beitrag von jogo »

Hallo Andreas,

meinst Du mit alles auch R und RStudio?
Ja, da gab es mal bei einem Versionssprung Änderungen in der Nutzung von Grafikbibliotheken.

Gruß, Jörg
Werekorden
Beiträge: 83
Registriert: So Feb 04, 2018 7:52 pm

Re: Logarithmische Regression

Beitrag von Werekorden »

Ja bei mir hat es am Windows-Rechner nur geholfen auf R 3.5.2 und auf R-Studio 1.1.463 zu updaten.

Bei meinem Mac klappen gerade die Updates der packages nicht, deshalb kann ich nicht sagen ob es eher die Packages sind oder die Software.
Werekorden
Beiträge: 83
Registriert: So Feb 04, 2018 7:52 pm

Re: Logarithmische Regression

Beitrag von Werekorden »

Hallo,


Ich habe noch immer das Problem, dass die Formel nicht so dargestellt wird wie sie es laut dem Code aus dem R Graphic CCookbook dargestellt werden sollte.

Bei der Ploterstellung sieht die Formel wie folgt aus: y=c(54.53)+c(-1.504)x, r^2 =0.99
Die beiden "c" sollten so in der Formel aber nicht auftauchen.

Hier nochmal mein kompletter Code inklusive des Codes aus dem Buch:

Code: Alles auswählen

library(Rcpp)
  library(rlang)
  library(gcookbook)
  library(ggplot2)
  library(hrbrthemes)
  library(magick)
  library(plyr)
  library(dplyr)
  library(MASS)
  library(scales)
  library(Cairo)
  library(extrafont)
  extrafont::loadfonts()
  font_import()
  loadfonts()

  #Dateneingabe
  Konz_x <- c(47143480, 4714348, 4714, 2357, 1179, 471)
  CT_Wert_y <- c(18.68, 22.595, 33.004, 34.290, 35.369, 35.271)
  data.frame(Konz_x, CT_Wert_y)
  p <- data.frame(Konz_x, CT_Wert_y)
  p
  
  #Model-Erstellung für Formeldarstellung in plot 
  model <- lm(CT_Wert_y ~ log(Konz_x))
  summary(model)
  model
  eqn <- as.character(as.expression(substitute(italic(y)==a+b*italic(x)*"," ~~ italic(r)^2~ "=" ~ r2,
         list(a = format(coef(model) [1], digits = 4), b = format(coef(model)[2], digits = 4), r2 = format(summary(model)$r.squared, digits = 2)
          ))))
  eqn
  parse(text = eqn)
  
 #Plot-Erstellung
  point <- ggplot(p, aes(x=Konz_x, y=CT_Wert_y)) 
  point + geom_point() + annotation_logticks(sides = "b")  + scale_x_log10() + stat_smooth(method=lm, se=0.99, colour="purple") +
    labs(title = 'Linearity Plot of FTD-Kit Data',
         subtitle = 'Linearity Plot to find Range of Linearity with Coefficient of Determination',
         x = 'Copies/ml patient sample',
         y = 'Ct-value',
         caption = 'NAME / NAME@ADRESSE.com') + 
    annotate("text", label = eqn, parse = TRUE, x=2600 , y=20.5) +
    theme_ipsum() + 
    theme(plot.title = element_text(color = "#3a4fa7ff"), plot.caption = element_text(color = "#3a4fa7ff", face = 'bold'))
  logo <-image_read("C:/User/Bilddatei") 
  grid::grid.raster(logo, x = 0.87, y = 0.88, just = c('left', 'bottom'), width = unit(0.8, 'inches'))

Vielleicht hat noch jemand eine Idee.
Athomas
Beiträge: 768
Registriert: Mo Feb 26, 2018 8:19 pm

Re: Logarithmische Regression

Beitrag von Athomas »

Code: Alles auswählen

eqn <- as.character(as.expression(substitute(italic(y)==a+b*italic(x)*"," ~~ italic(r)^2~ "=" ~ r2,
                                             list(a = format(coef(model)[[1]], digits = 4), b = format(coef(model)[[2]], digits = 4), r2 = format(summary(model)$r.squared, digits = 2)
                                             ))))
oder besser

Code: Alles auswählen

eqn <- as.character(as.expression(substitute(italic(y)==a+b*italic(x)*"," ~~ italic(r)^2~ "=" ~ r2,
                                             list(a  = format(summary(model)$coefficients[1,1], digits = 4), 
                                                  b  = format(summary(model)$coefficients[2,1], digits = 4), 
                                                  r2 = format(summary(model)$r.squared, digits = 2)
                                             ))))
Werekorden
Beiträge: 83
Registriert: So Feb 04, 2018 7:52 pm

Re: Logarithmische Regression

Beitrag von Werekorden »

Hi Athomas,

Danke für den Code und besonderen Dank für die geordnete Version ;)

Ich bin gerade überrascht, dass man sowohl eine Zeile und Spalte angeben muss um das "c" wegzubekommen und die entsprechende Variable allein zu erhalten, obwohl es nur eine Zeile gibt.

Ist das immer so? Sorry aber ich will gerne wissen warum es so funktioniert und so wie im Buch nicht funktioniert.


Dank dir auf jeden Fall schonmal.

VG
Andreas
Athomas
Beiträge: 768
Registriert: Mo Feb 26, 2018 8:19 pm

Re: Logarithmische Regression

Beitrag von Athomas »

Sorry aber ich will gerne wissen warum es so funktioniert und so wie im Buch nicht funktioniert.
Ich glaube, die Erklärung wird nicht einfach zu finden sein!
Wenn es Dich so brennend interessiert, hilft nur eins: forschen, forschen, forschen :lol: !
Antworten