Add result folder

adultabundance/run_adultabundance.R

@@ -99,5 +99,17 @@ cat("Sample analyzed after ", elapsed.time, ' minutes\n')
 save(fit.mcmc,file=paste('result/result',model_name,'.RData',sep=""))
 
 
+## STATS
+
+load(file=paste('result/result',model_name,'.RData',sep=""))
+fit.mat <- as.matrix(fit.mcmc)
+
+
+mydf <- cbind(summary(fit.mcmc)$statistics[,c("Mean","SD")],summary(fit.mcmc)$quantiles)
+mydf <- cbind(rownames(mydf), mydf)
+rownames(mydf) <- NULL
+colnames(mydf)[1] <- c("Variable")
+write.table(mydf,file=paste('result/stat',model_name,'.csv',sep=""),sep=",", row.names = FALSE)
+