Add output of model mean an error for the prediction from nested sampling

Bojan Nikolic · Bojan Nikolic · commit 9dd53c2a3da9 · 2012-11-01T19:22:53.000Z
diff --git a/R/output.R b/R/output.R
@@ -130,4 +130,8 @@ nested.mhist <- function(r, m, x)
                   w,
                   main=paste("Probability distribution of model at ", x),
                   breaks="FD")
+
+    x1 <- sum(pw*w)/sum(w)
+    x2 <- sum((pw-x1)**2*w)/sum(w)
+    print(paste("*** Model prediction and error: ", x1, " +/- ", x2))
   }

Original file line number	Diff line number	Diff line change
`@@ -130,4 +130,8 @@ nested.mhist <- function(r, m, x)`
`130`	`130`	`w,`
`131`	`131`	`main=paste("Probability distribution of model at ", x),`
`132`	`132`	`breaks="FD")`
	`133`	`+`
	`134`	`+ x1 <- sum(pw*w)/sum(w)`
	`135`	`+ x2 <- sum((pw-x1)*2w)/sum(w)`
	`136`	`+ print(paste("*** Model prediction and error: ", x1, " +/- ", x2))`
`133`	`137`	`}`