commit da836be5f37fbc90baee74fd260b650f77932a5a
parent 5f82063967a4d453cfb43477a586b3ecde0c806e
Author: eamoncaddigan <eamon.caddigan@gmail.com>
Date: Fri, 21 Aug 2015 16:46:35 -0400
Wasn't using smryMCMC.
Diffstat:
1 file changed, 0 insertions(+), 27 deletions(-)
diff --git a/Jags-Yord-Xnom1grp-Mnormal.R b/Jags-Yord-Xnom1grp-Mnormal.R
@@ -112,31 +112,6 @@ genMCMC = function( datFrm, yName , qName,
#===============================================================================
-smryMCMC = function( codaSamples , compVal , #RopeEff=NULL ,
- saveName=NULL ) {
- summaryInfo = NULL
- mcmcMat = as.matrix(codaSamples,chains=TRUE)
- summaryInfo = rbind( summaryInfo ,
- "mu" = summarizePost( mcmcMat[,"mu"] ) )
- summaryInfo = rbind( summaryInfo ,
- "sigma" = summarizePost( mcmcMat[,"sigma"] ) )
- summaryInfo = rbind( summaryInfo ,
- "effSz" = summarizePost(
- ( mcmcMat[,"mu"] - compVal ) / mcmcMat[,"sigma"] ,
- compVal=compVal , ROPE=NULL ) )
- for ( colName in grep( "thresh" , colnames(mcmcMat) , value=TRUE ) ) {
- summaryInfo = rbind( summaryInfo ,
- summarizePost( mcmcMat[,colName] ) )
- rownames(summaryInfo)[nrow(summaryInfo)] = colName
- }
- if ( !is.null(saveName) ) {
- write.csv( summaryInfo , file=paste(saveName,"SummaryInfo.csv",sep="") )
- }
- return( summaryInfo )
-}
-
-#===============================================================================
-
plotMCMC = function( codaSamples , datFrm , yName , qName, compVal , #RopeEff=NULL ,
showCurve=FALSE ,
saveName=NULL , saveType="jpg" ) {
@@ -146,8 +121,6 @@ plotMCMC = function( codaSamples , datFrm , yName , qName, compVal , #RopeEff=NU
q = as.numeric(as.factor(datFrm[[qName]]))
nQlevels = max(q)
-
-
#-----------------------------------------------------------------------------
# Plots for each question
for (i in 1:nQlevels) {