standard deviation of the absolute value of the log

git-svn-id: svn+ssh://hemlock.ucr.edu/srv/svn/misc/trunk/origins_of_life@526 25fa0111-c432-4dab-af88-9f31a2f6ac42

diff --git a/kinetic_formalism.Rnw b/kinetic_formalism.Rnw
index cd5577541a98cb7541793ba3a7cf8801b51dbb15..a1fa109a5e3f9c7652ef5d1cbbac263af879eb20 100644 (file)
--- a/kinetic_formalism.Rnw
+++ b/kinetic_formalism.Rnw
@@ -250,7 +250,7 @@ is $10$, yielding:
 
 \begin{equation}
  % cu_f = 10^{\mathrm{stdev}\left|\log cu_\mathrm{vesicle}\right|}
-  cu_f = 10^{\left|\left<\log cu_\mathrm{vesicle} \right>\right|\mathrm{stdev} \log cu_\mathrm{vesicle}}
+  cu_f = 10^{\left|\left<\log cu_\mathrm{vesicle} \right>\right|\mathrm{stdev} \left|\log cu_\mathrm{vesicle}\right|}
   \label{eq:curvature_forward}
 \end{equation}
 
@@ -262,8 +262,8 @@ $\Sexpr{format(digits=3,to.kcal(10^(0.13*0.213)))}
 
 % 1.5 to 0.75 3 to 0.33
 <<fig=TRUE,echo=FALSE,results=hide,width=7,height=7>>=
-grid <- expand.grid(x=seq(0,max(c(sd(log(c(1,3))),
-                      sd(log(c(1,0.33))),sd(log(c(0.33,3))))),length.out=20),
+grid <- expand.grid(x=seq(0,max(c(sd(abs(log(c(1,3)))),
+                      sd(abs(log(c(1,0.33)))),sd(abs(log(c(0.33,3)))))),length.out=20),
                     y=seq(0,max(c(mean(log(c(1,3)),
                       mean(log(c(1,0.33))),
                       mean(log(c(0.33,3)))))),length.out=20))
@@ -277,8 +277,8 @@ print(wireframe(z~x*y,grid,cuts=50,
 rm(grid)
 @ 
 <<fig=TRUE,echo=FALSE,results=hide,width=7,height=7>>=
-grid <- expand.grid(x=seq(0,max(c(sd(log(c(1,3))),
-                      sd(log(c(1,0.33))),sd(log(c(0.33,3))))),length.out=20),
+grid <- expand.grid(x=seq(0,max(c(sd(abs(log(c(1,3)))),
+                      sd(abs(log(c(1,0.33)))),sd(abs(log(c(0.33,3)))))),length.out=20),
                     y=seq(0,max(c(mean(log(c(1,3)),
                       mean(log(c(1,0.33))),
                       mean(log(c(0.33,3)))))),length.out=20))