[Mattice-commits] r253 - pkg/R

[noreply at r-forge.r-project.org]

Author: andrew_hipp
Date: 2012-05-03 04:15:18 +0200 (Thu, 03 May 2012)
New Revision: 253

Modified:
   pkg/R/ouSim.hansenBatch.R
Log:
Fixing error in which alpha and sigma.sq are being passed as mean +/- CI to ouSim rather than just the mean

Modified: pkg/R/ouSim.hansenBatch.R
===================================================================
--- pkg/R/ouSim.hansenBatch.R	2012-05-02 21:34:37 UTC (rev 252)
+++ pkg/R/ouSim.hansenBatch.R	2012-05-03 02:15:18 UTC (rev 253)
@@ -3,7 +3,7 @@
   analysis <- object
   # if(class(analysis) == "hansenBatch") analysis <- summary(analysis)
   if(identical(rootState, NULL)) rootState <- analysis$thetaMatrix[treeNum, ][tree at root] # rootstate taken to be the optimum at the root
-  outdata <- ouSim(tree, rootState, sqrt.alpha = analysis$modelAvgAlpha, variance = analysis$modelAvgSigmaSq, theta = analysis$thetaMatrix[treeNum, ], ...)
+  outdata <- ouSim(tree, rootState, sqrt.alpha = analysis$modelAvgAlpha['mean'], variance = analysis$modelAvgSigmaSq['mean'], theta = analysis$thetaMatrix[treeNum, ], ...)
   class(outdata) <- "ouSim"
   return(outdata)
 }