[chronojump] MIF. Plotting tau

[Xavier Padullés <xpadulles src gnome org>]

commit 9c381357566048fe166e005a9eab6ac7de09cc08
Author: Xavier Padullés <x padulles gmail com>
Date:   Fri Dec 1 19:26:10 2017 +0100

    MIF. Plotting tau

 r-scripts/maximumIsometricForce.R |   19 +++++++++++++++++--
 1 files changed, 17 insertions(+), 2 deletions(-)
---
diff --git a/r-scripts/maximumIsometricForce.R b/r-scripts/maximumIsometricForce.R
index 83dcbbd..d1d83cb 100644
--- a/r-scripts/maximumIsometricForce.R
+++ b/r-scripts/maximumIsometricForce.R
@@ -314,7 +314,22 @@ drawDynamicsFromLoadCell <- function(
         
         #Plotting RFD
         #lines(dynamics$time, dynamics$rfd/10)
-
+        
+        #Plotting tau
+        abline(v = 0, col = "green", lty = 3)
+        abline(v = dynamics$tau.fitted, col = "green", lty = 3)
+        arrows(x0 = 0, y0 = dynamics$fmax.fitted*0.6321206,
+               x1 = dynamics$tau.fitted, y1 = dynamics$fmax.fitted*0.6321206)
+        text(x = (dynamics$tau.fitted / 2), y = dynamics$fmax.fitted*0.6321206,
+             labels = paste("τ =", round(dynamics$tau.fitted, digits = 2), "s"), pos = 3, cex = 1.5)
+        
+        arrows(x0 = dynamics$tau.fitted, y0 = 0,
+               x1 = dynamics$tau.fitted, y1 = dynamics$fmax.fitted*0.6321206)
+        
+        text(x = (dynamics$tau.fitted), y = dynamics$fmax.fitted*0.6321206 / 2,
+              labels = "63% of fmax", pos = 4, cex = 1.5)
+        
+        #Plotting fmax.raw
         text( x = dynamics$tfmax.raw, y = dynamics$fmax.raw,
               labels = paste("Fmax = ", round(dynamics$fmax.raw, digits=2), " N", sep=""), pos = 3, cex = 
1.5)
         points(x = dynamics$tfmax.raw, y = dynamics$fmax.raw)
@@ -354,7 +369,7 @@ drawDynamicsFromLoadCell <- function(
         legendText = c(
                       paste("Fmax =", round(dynamics$fmax.fitted, digits = 2), "N"),
                       paste("K = ", round(dynamics$k.fitted, digits = 2),"s⁻¹"),
-                      paste("Tau = ", round(dynamics$tau.fitted, digits = 2),"s")
+                      paste("τ = ", round(dynamics$tau.fitted, digits = 2),"s")
                       )
         legendColor = c("blue", "blue")