R/brier_score.R

Defines functions brier.score

Documented in brier.score

## goal: to summarize pi^ (predicted P[Y=1]) using different scores

## for 2 classes only
## Brier score = sum_i (y_i - p_i^)^2

brier.score <- function (predictedArr, truthArr) {
	## make sure the 2 arrays are the same length
	if (length(predictedArr) != length(truthArr)) {
	   print ("ERROR: length NOT equal!!")
        }
 	## sum over all the samples
	temp.vec <- (truthArr - predictedArr)^2
	sum (temp.vec)
}

Try the iterativeBMA package in your browser

Any scripts or data that you put into this service are public.

iterativeBMA documentation built on Nov. 17, 2017, 1:26 p.m.