Skip to contents

Creates an object representing the hyper-g/n mixture of g-priors on coefficients for BAS. This is a special case of the tCCH prior

Usage

hyper.g.n(alpha = 3, n = NULL)

Arguments

alpha

a scalar > 0, recommended 2 < alpha <= 3

n

The sample size; if NULL, the value derived from the data in the call to `bas.glm` will be used.

Value

returns an object of class "prior", with the family and hyerparameters.

Details

Creates a structure used for bas.glm. This is a special case of the tCCH, where hyper.g.n(alpha=3, n) is equivalent to tCCH(alpha=1, beta=2, s=0, r=1.5, v = 1, theta=1/n)

See also

Author

Merlise Clyde

Examples

n <- 500
hyper.g.n(alpha = 3, n = n)
#> $family
#> [1] "hyper-g/n"
#> 
#> $class
#> [1] "TCCH"
#> 
#> $hyper.parameters
#> $hyper.parameters$alpha
#> [1] 1
#> 
#> $hyper.parameters$beta
#> [1] 2
#> 
#> $hyper.parameters$s
#> [1] 0
#> 
#> $hyper.parameters$r
#> [1] 1.5
#> 
#> $hyper.parameters$v
#> [1] 1
#> 
#> $hyper.parameters$theta
#> [1] 0.002
#> 
#> 
#> $n
#> [1] 500
#> 
#> attr(,"class")
#> [1] "prior"