Added some descriptive statistics verification tests.
git-svn-id: https://svn.apache.org/repos/asf/commons/proper/math/trunk@600358 13f79535-47bb-0310-9956-ffa450edef68
This commit is contained in:
parent
8140c9da58
commit
ef6b260fe9
|
@ -76,7 +76,7 @@ a browser instance and then entering
|
|||
|
||||
at an R prompt. Poking about in the test case files and the online docs should
|
||||
bring you up to speed fairly quickly. Here are some basic things to get
|
||||
you started. I should note at this point that I by no means an expert R
|
||||
you started. I should note at this point that I am by no means an expert R
|
||||
programmer, so some things may not be implemented in the the nicest way.
|
||||
Comments / suggestions for improvement are welcome!
|
||||
|
||||
|
|
|
@ -0,0 +1,83 @@
|
|||
# Licensed to the Apache Software Foundation (ASF) under one or more
|
||||
# contributor license agreements. See the NOTICE file distributed with
|
||||
# this work for additional information regarding copyright ownership.
|
||||
# The ASF licenses this file to You under the Apache License, Version 2.0
|
||||
# (the "License"); you may not use this file except in compliance with
|
||||
# the License. You may obtain a copy of the License at
|
||||
#
|
||||
# http://www.apache.org/licenses/LICENSE-2.0
|
||||
#
|
||||
# Unless required by applicable law or agreed to in writing, software
|
||||
# distributed under the License is distributed on an "AS IS" BASIS,
|
||||
# WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
||||
# See the License for the specific language governing permissions and
|
||||
# limitations under the License.
|
||||
#
|
||||
#------------------------------------------------------------------------------
|
||||
# R source file to validate summary statistics tests in
|
||||
# org.apache.commons.math.stat.CertifiedDataTest
|
||||
#
|
||||
# To run the test, install R, put this file and testFunctions
|
||||
# into the same directory, launch R from this directory and then enter
|
||||
# source("<name-of-this-file>")
|
||||
#
|
||||
# R functions used
|
||||
# mean(x)
|
||||
# sd(x)
|
||||
#------------------------------------------------------------------------------
|
||||
tol <- 1E-14 # error tolerance for tests
|
||||
#------------------------------------------------------------------------------
|
||||
# Function definitions
|
||||
|
||||
source("testFunctions") # utility test functions
|
||||
options(digits=15) # bump displayed digits to 15
|
||||
|
||||
verifyMean <- function(values, expectedMean, tol, desc) {
|
||||
results <- mean(values)
|
||||
if (assertEquals(expectedMean, results, tol, "mean")) {
|
||||
displayPadded(c(desc," mean test"), SUCCEEDED, WIDTH)
|
||||
} else {
|
||||
displayPadded(c(desc, " mean test"), FAILED, WIDTH)
|
||||
}
|
||||
}
|
||||
|
||||
verifySigma <- function(values, expectedSigma, tol, desc) {
|
||||
results <- sd(values)
|
||||
if (assertEquals(expectedSigma, results, tol, "std")) {
|
||||
displayPadded(c(desc," std test"), SUCCEEDED, WIDTH)
|
||||
} else {
|
||||
displayPadded(c(desc, " std test"), FAILED, WIDTH)
|
||||
}
|
||||
}
|
||||
|
||||
cat("Descriptive test cases\n")
|
||||
|
||||
# Michelson data
|
||||
values <- c(299.85,299.74,299.90,300.07,299.93,299.85,299.95,299.98,299.98,
|
||||
299.88,300.00,299.98,299.93,299.65,299.76,299.81,300.00,300.00,299.96,299.96,
|
||||
299.96,299.94,299.96,299.94,299.88,299.80,299.85,299.88,299.90,299.84,299.83,
|
||||
299.79,299.81,299.88,299.88,299.83,299.80,299.79,299.76,299.80,299.88,299.88,
|
||||
299.88,299.86,299.72,299.72,299.62,299.86,299.97,299.95,299.88,299.91,299.85,
|
||||
299.87,299.84,299.84,299.85,299.84,299.84,299.84,299.89,299.81,299.81,299.82,
|
||||
299.80,299.77,299.76,299.74,299.75,299.76,299.91,299.92,299.89,299.86,299.88,
|
||||
299.72,299.84,299.85,299.85,299.78,299.89,299.84,299.78,299.81,299.76,299.81,
|
||||
299.79,299.81,299.82,299.85,299.87,299.87,299.81,299.74,299.81,299.94,299.95,
|
||||
299.80,299.81,299.87)
|
||||
expectedMean <- 299.852400000000
|
||||
expectedSigma <- 0.0790105478190518
|
||||
verifyMean(values, expectedMean, tol, "Michelson")
|
||||
verifySigma(values, expectedSigma, tol, "Michelson")
|
||||
|
||||
# Mavro data
|
||||
values <- c(2.00180,2.00170,2.00180,2.00190,2.00180,2.00170,2.00150,2.00140,
|
||||
2.00150,2.00150,2.00170,2.00180,2.00180,2.00190,2.00190,2.00210,2.00200,
|
||||
2.00160,2.00140,2.00130,2.00130,2.00150,2.00150,2.00160,2.00150,2.00140,
|
||||
2.00130,2.00140,2.00150,2.00140,2.00150,2.00160,2.00150,2.00160,2.00190,
|
||||
2.00200,2.00200,2.00210,2.00220,2.00230,2.00240,2.00250,2.00270,2.00260,
|
||||
2.00260,2.00260,2.00270,2.00260,2.00250,2.00240)
|
||||
expectedMean <- 2.00185600000000
|
||||
expectedSigma <- 0.000429123454003053
|
||||
verifyMean(values, expectedMean, tol, "Mavro")
|
||||
verifySigma(values, expectedSigma, tol, "Mavro")
|
||||
|
||||
displayDashes(WIDTH)
|
|
@ -39,6 +39,10 @@ source("regressionTestCases")
|
|||
|
||||
# inference
|
||||
source("chiSquareTestCases")
|
||||
|
||||
# descriptive
|
||||
source("descriptiveTestCases")
|
||||
|
||||
#------------------------------------------------------------------------------
|
||||
# if output has been diverted, change it back
|
||||
if (sink.number()) {
|
||||
|
|
Loading…
Reference in New Issue