MyNixOS website logo
Description

Statistical Tests of Familial Hypotheses.

Provides functionality for testing familial hypotheses. Supports testing centers belonging to the Huber family. Testing is carried out using the Bayesian bootstrap. One- and two-sample tests are supported, as are directional tests. Methods for visualizing output are provided.

familial

R-CMD-check codecov

Overview

An R package for familial inference. Briefly, this package provides tests for hypotheses of the form

\mathrm{H}_0:\mu(\lambda)=\mu_0\text{ for some }\lambda\in\Lambda\quad\text{vs.}\quad\mathrm{H}_1:\mu(\lambda)\neq\mu_0\text{ for all }\lambda\in\Lambda,

where \{\mu(\lambda):\lambda\in\Lambda\} is a family of centers, e.g., that induced by the Huber loss function with parameter \lambda. In contrast to classic statistical tests such as the t or sign tests for the mean or median, familial tests do not depend on a single (sometimes arbitrarily chosen) center.

Presently, familial supports tests of the Huber family of centers, which includes the mean and median. Testing is carried out using a Bayesian approach whereby the posterior probabilities of the competing hypotheses \mathrm{H}_0 and \mathrm{H}_1 are from the Bayesian bootstrap. One- and two-sample tests are supported, as are directional tests. Methods for visualizing output are provided.

Installation

To install the latest stable version from CRAN, run the following code:

install.packages('familial')

To install the latest development version from GitHub, run the following code:

devtools::install_github('ryan-thompson/familial')

Usage

The center.test() function performs a test of centers, with the default being the Huber family of centers.

library(familial)
set.seed(1)

# One-sample test with point null
x <- MASS::galaxies
center.test(x, mu = 21000)
## -----------------------------------------------
## familial test of centers with huber family
## -----------------------------------------------
## mu = 21000 
## posterior probabilities: 
##    H0    H1 
## 0.542 0.458 
## optimal decision: indeterminate
# One-sample test with interval null
center.test(x, mu = c(20500, 21500))
## -----------------------------------------------
## familial test of centers with huber family
## -----------------------------------------------
## mu = 20500 21500 
## posterior probabilities: 
##    H0    H1 
## 0.959 0.041 
## optimal decision: H0
# Two-sample test
x <- MASS::cabbages[MASS::cabbages$Cult == 'c39', 'HeadWt']
y <- MASS::cabbages[MASS::cabbages$Cult == 'c52', 'HeadWt']
center.test(x, y)
## -----------------------------------------------
## familial test of centers with huber family
## -----------------------------------------------
## mu = 0 
## posterior probabilities: 
##    H0    H1 
## 0.008 0.992 
## optimal decision: H1
# Two-sample paired directional test
x <- MASS::anorexia[MASS::anorexia$Treat == 'FT', 'Postwt']
y <- MASS::anorexia[MASS::anorexia$Treat == 'FT', 'Prewt']
center.test(x, y, paired = T, alternative = 'greater')
## -----------------------------------------------
## familial test of centers with huber family
## -----------------------------------------------
## mu = 0 
## posterior probabilities: 
##    H0    H1 
## 0.006 0.994 
## optimal decision: H1

Documentation

See the package vignette or reference manual.

Metadata

Version

1.0.5

License

Unknown

Platforms (75)

    Darwin
    FreeBSD
    Genode
    GHCJS
    Linux
    MMIXware
    NetBSD
    none
    OpenBSD
    Redox
    Solaris
    WASI
    Windows
Show all
  • aarch64-darwin
  • aarch64-genode
  • aarch64-linux
  • aarch64-netbsd
  • aarch64-none
  • aarch64_be-none
  • arm-none
  • armv5tel-linux
  • armv6l-linux
  • armv6l-netbsd
  • armv6l-none
  • armv7a-darwin
  • armv7a-linux
  • armv7a-netbsd
  • armv7l-linux
  • armv7l-netbsd
  • avr-none
  • i686-cygwin
  • i686-darwin
  • i686-freebsd
  • i686-genode
  • i686-linux
  • i686-netbsd
  • i686-none
  • i686-openbsd
  • i686-windows
  • javascript-ghcjs
  • loongarch64-linux
  • m68k-linux
  • m68k-netbsd
  • m68k-none
  • microblaze-linux
  • microblaze-none
  • microblazeel-linux
  • microblazeel-none
  • mips-linux
  • mips-none
  • mips64-linux
  • mips64-none
  • mips64el-linux
  • mipsel-linux
  • mipsel-netbsd
  • mmix-mmixware
  • msp430-none
  • or1k-none
  • powerpc-netbsd
  • powerpc-none
  • powerpc64-linux
  • powerpc64le-linux
  • powerpcle-none
  • riscv32-linux
  • riscv32-netbsd
  • riscv32-none
  • riscv64-linux
  • riscv64-netbsd
  • riscv64-none
  • rx-none
  • s390-linux
  • s390-none
  • s390x-linux
  • s390x-none
  • vc4-none
  • wasm32-wasi
  • wasm64-wasi
  • x86_64-cygwin
  • x86_64-darwin
  • x86_64-freebsd
  • x86_64-genode
  • x86_64-linux
  • x86_64-netbsd
  • x86_64-none
  • x86_64-openbsd
  • x86_64-redox
  • x86_64-solaris
  • x86_64-windows