From bc678b2525257a40c615cb7487acb4fcae132a31 Mon Sep 17 00:00:00 2001 From: Franz Innerbichler Date: Sun, 10 May 2026 11:40:46 +0200 Subject: [PATCH] =?UTF-8?q?Bug=20fix=20f=C3=BCr=20linear=20reg=20bei=20XL?= =?UTF-8?q?=20upload?= MIME-Version: 1.0 Content-Type: text/plain; charset=UTF-8 Content-Transfer-Encoding: 8bit --- .gitignore | 4 ++++ app.R | 31 ++++++++++++++++--------------- 2 files changed, 20 insertions(+), 15 deletions(-) create mode 100644 .gitignore diff --git a/.gitignore b/.gitignore new file mode 100644 index 0000000..5b6a065 --- /dev/null +++ b/.gitignore @@ -0,0 +1,4 @@ +.Rproj.user +.Rhistory +.RData +.Ruserdata diff --git a/app.R b/app.R index ea7f1b6..45f14b1 100644 --- a/app.R +++ b/app.R @@ -2164,7 +2164,8 @@ server <- function(input, output, session) { # tab <- sim2() if (is.character(tab)) stop(tab) - browser() + #browser() + log_conc <- tab$log_dose noDilSer = (ncol(tab)-1)/2 noDil <- nrow(tab) Conctab <- perConcTab(tab, noDilSer) @@ -2173,11 +2174,11 @@ server <- function(input, output, session) { slopeSt <- slopeTe <- matrix(NA, nrow=noDil-2,ncol=2) for (i in 1:(noDil-2)) { avs <- Conctab[noDilSer+1,] - threes <- data.frame(lnC=log(Conc[i:(i+2)]), resp=avs[i:(i+2)]) + threes <- data.frame(lnC=log_conc[i:(i+2)], resp=avs[i:(i+2)]) lm3St <- lm(resp ~ lnC, data=threes) slopeSt[i,] <- lm3St$coefficients avt <- Conctab[noDilSer*2+4,] - threet <- data.frame(lnC=log(Conc[i:(i+2)]), resp=avt[i:(i+2)]) + threet <- data.frame(lnC=log_conc[i:(i+2)], resp=avt[i:(i+2)]) lm3Te <- lm(resp ~ lnC, data=threet) slopeTe[i,] <- lm3Te$coefficients } @@ -2185,9 +2186,9 @@ server <- function(input, output, session) { indS <- which(abs(slopeSt[,2]) == max(abs(slopeSt[,2]))) indT <- which(abs(slopeTe[,2]) == max(abs(slopeTe[,2]))) - pl_ <- slopeSt[indS,1]+slopeSt[indS,2]*log(Conc) - pl_T <- slopeTe[indT,1]+slopeTe[indT,2]*log(Conc) - pl_df <- data.frame(lnC=log(Conc), plotS=pl_, plotT=pl_T) + pl_ <- slopeSt[indS,1]+slopeSt[indS,2]*log_conc + pl_T <- slopeTe[indT,1]+slopeTe[indT,2]*log_conc + pl_df <- data.frame(lnC=log_conc, plotS=pl_, plotT=pl_T) all_l <- melt(data.frame(tab), id.vars="log_dose",variable.name="replname",value.name="readout") isRef <- rep(c(1,0), 1,each=nrow(all_l)/2) @@ -2240,9 +2241,9 @@ server <- function(input, output, session) { # fit intercept for test and ref and common slope - pl_restS <- sum_mLin$coefficients[1,1] + sum_mLin$coefficients[2,1]*log(Conc) - pl_restT <- sum_mLin$coefficients[1,1] + sum_mLin$coefficients[3,1] + sum_mLin$coefficients[2,1]*log(Conc) - pl_rest <- data.frame(lnC=log(Conc), plotS=pl_restS, plotT=pl_restT) + pl_restS <- sum_mLin$coefficients[1,1] + sum_mLin$coefficients[2,1]*log_conc + pl_restT <- sum_mLin$coefficients[1,1] + sum_mLin$coefficients[3,1] + sum_mLin$coefficients[2,1]*log_conc + pl_rest <- data.frame(lnC=log_conc, plotS=pl_restS, plotT=pl_restT) pr2 <- p + geom_line(data=pl_rest,aes(x=lnC,y=plotS),color="#4545BA", inherit.aes = F) + @@ -2285,9 +2286,9 @@ server <- function(input, output, session) { indS <- which(abs(slopeSt[,2]) == max(abs(slopeSt[,2]))) indT <- which(abs(slopeTe[,2]) == max(abs(slopeTe[,2]))) - pl_ <- slopeSt[indS,1]+slopeSt[indS,2]*log(Conc) - pl_T <- slopeTe[indT,1]+slopeTe[indT,2]*log(Conc) - pl_df <- data.frame(lnC=log(Conc), plotS=pl_, plotT=pl_T) + pl_ <- slopeSt[indS,1]+slopeSt[indS,2]*log_conc + pl_T <- slopeTe[indT,1]+slopeTe[indT,2]*log_conc + pl_df <- data.frame(lnC=log_conc, plotS=pl_, plotT=pl_T) all_l <- melt(data.frame(tab), id.vars="log_dose",variable.name="replname",value.name="readout") isRef <- rep(c(1,0), 1,each=nrow(all_l)/2) @@ -2338,9 +2339,9 @@ server <- function(input, output, session) { # alternativ: modAB <- lm(readout ~ log_dose+isSample, circle) sum_mLin <- summary(mLin) - pl_restS <- sum_mLin$coefficients[1,1] + sum_mLin$coefficients[2,1]*log(Conc) - pl_restT <- sum_mLin$coefficients[1,1] + sum_mLin$coefficients[3,1] + sum_mLin$coefficients[2,1]*log(Conc) - pl_rest <- data.frame(lnC=log(Conc), plotS=pl_restS, plotT=pl_restT) + pl_restS <- sum_mLin$coefficients[1,1] + sum_mLin$coefficients[2,1]*log_conc + pl_restT <- sum_mLin$coefficients[1,1] + sum_mLin$coefficients[3,1] + sum_mLin$coefficients[2,1]*log_conc + pl_rest <- data.frame(lnC=log_conc, plotS=pl_restS, plotT=pl_restT) pr2 <- p + geom_line(data=pl_rest,aes(x=lnC,y=plotS),color="blue", inherit.aes = F) +