Deutsches R-Forum

Code: Alles auswählen

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)

fm <- lm(CT_Wert_y ~ log(Konz_x))
summary (fm)
fm
plot(Konz_x, CT_Wert_y, xlim = c(300, 50000000), ylim = c(0, 40), log="x", xlab = "Copies/ml patient sample", 
     ylab = "Ct-Values", main ="Linear regression CAP Kit", pch = 19)
# abline(fm, col = "red",lwd=2, untf = FALSE)

b <- coef(fm); b0 <- b[1]; b1 <- b[2]
mycurve <- function(x) b0 + b1*log(x)

curve(mycurve, 300, 50000000, add=TRUE, col="green")

Code: Alles auswählen

b <- coef(fm); b0 <- b[1]; b1 <- b[2]
curve(b0 + b1*log(x), 300, 50000000, add=TRUE, col="blue")
# oder
curve(cbind(1, log(x)) %*% coef(fm), 300, 50000000, add=TRUE, col="yellow")
# oder
curve(tcrossprod(coef(fm), cbind(1, log(x))), 300, 50000000, add=TRUE, col="brown")

Code: Alles auswählen

lines(Konz_x, fitted(fm), col="orange")

Code: Alles auswählen

axis(1, at=c(1e3, 1e4, 1e5, 1e6, 1e7,1e8), 
             labels=FALSE, lwd.ticks=2)
axis(1, at=c(seq(1e2, 1e3, 1e2), 
             seq(1e3, 1e4, 1e3),
             seq(1e4, 1e5, 1e4),
             seq(1e5, 1e6, 1e5),
             seq(1e6, 1e7, 1e6),
             seq(1e7, 1e8, 1e7)
             ), labels = FALSE)

Code: Alles auswählen

library(Rcpp)
library(rlang)
library(gcookbook)
library(ggplot2)
library(hrbrthemes)
library(magick)
library(plyr)
library(dplyr)
library(MASS)
library(scales)

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
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="red") +
  labs(title = 'Linearity Plot of FTD Kit Data',
       subtitle = 'Linearity Plot with Coefficient of Determination',
       x = 'Copies/ml patient sample',
       y = 'Ct-value',
       caption = 'NAME HIER / ??@??.com') + 
  theme_ipsum() + 
  theme(plot.title = element_text(color = "#EC0108"), plot.caption = element_text(color = "#EC0108", face = 'bold'))
logo <-image_read("/Users/NAME/Pictures/Dateiname.png")
grid::grid.raster(logo, x = 0.07, y = 0.01, just = c('left', 'bottom'), width = unit(0.5, 'inches'))

Code: Alles auswählen

  theme_ipsum

Code: Alles auswählen

extrafont::loadfonts()

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'))