[chronojump] Fixed time to pmax in photocells



commit 71c2c4248f06fb13e1defb35fe2092fb0297d113
Author: Xavier Padullés <testing chronojump org>
Date:   Sat Jan 16 01:50:21 2021 +0100

    Fixed time to pmax in photocells

 r-scripts/sprintPhotocells.R | 4 ++--
 1 file changed, 2 insertions(+), 2 deletions(-)
---
diff --git a/r-scripts/sprintPhotocells.R b/r-scripts/sprintPhotocells.R
index 3b061858..c78be9c4 100644
--- a/r-scripts/sprintPhotocells.R
+++ b/r-scripts/sprintPhotocells.R
@@ -238,9 +238,9 @@ drawSprintFromPhotocells <- function(sprintDynamics, splitTimes, positions, titl
                 plot((sprintDynamics$t.fitted - sprintDynamics$T0), sprintDynamics$p.fitted, type="l", axes 
= FALSE, xlab="", ylab="", col="red"
                      , ylim=c(0,sprintDynamics$pmax.fitted + .1 * sprintDynamics$pmax.fitted), xlim = xlims
                      , yaxs= "i", xaxs = "i")
-                abline(v = sprintDynamics$tpmax.fitted, col="red", lty = 2)
+                abline(v = (sprintDynamics$tpmax.fitted- sprintDynamics$T0), col="red", lty = 2)
                 axis(line = 5, side = 4, col ="red", at = seq(0, sprintDynamics$pmax.fitted, by = 200))
-                axis(3, at = sprintDynamics$tpmax.fitted, labels = round(sprintDynamics$tpmax.fitted, 3))
+                axis(3, at = sprintDynamics$tpmax.fitted- sprintDynamics$T0, labels = 
round(sprintDynamics$tpmax.fitted, 3))
                 #text(sprintDynamics$tpmax.fitted, sprintDynamics$pmax.fitted, paste("Pmax fitted =", 
round(sprintDynamics$pmax.fitted, digits = 2)),  pos = 3)
                 # mtext(side = 1, line = 5, at = splitTimes[length(splitTimes)]*0.75, cex = 1.5,
                 #       substitute(P(t) == A*e^(-K*t)*(1-e^(-K*t)) + B*(1-e^(-K*t))^3,


[Date Prev][Date Next]   [Thread Prev][Thread Next]   [Thread Index] [Date Index] [Author Index]