ououbmprior: Draw prior samples for OUOUBM model

In ouxy: Model of Adaptive Trait Evolution

Description

Simulate sample for parameters in OUOUBM model given a set of hyper parameters

Usage

ououbmprior(
  prior.model.params = prior.model.params,
  prior.reg.params = prior.reg.params
)

Arguments

prior.model.params	A vectors of hyper parameters for model parameter containing (alpha.y.min, alpha.y.max) for alpha.y, (alpha.x.min, alpha.x.max) for alpha.x,(theta.y.min, theta.y.max) for theta.y, (sigmasq.x.min, sigmasq.x.max) for sigmasq.x, (tau.y.min, tau.max) for rate parameter of tau
prior.reg.params	A vector of hyper paramter for regression parameters. (b0.min,b0.max) for b0, (b1.min,b1.max) for b1, (b2.min,b2.max) for b2

Details

The function requires user to input hyper parameters for α_y,α_x,θ_x, σ^2_x,τ and hyper parameters for regression parameters b_0,b_1,b_2 from uniform distribution with its minimum and maximum values.

Value

Returns the samples of model parameters and regression parameter

Examples

prior.model.params<-c(0,3,0,3,-5,5,0,1,0,2)
names(prior.model.params)<-c(

"alpha.y.min","alpha.y.max","alpha.x.min","alpha.x.max",

"theta.x.min","theta.x.max","sigmasq.x.min","sigmasq.x.max",

"tau.min","tau.max")
prior.reg.params<-c(-3, 3, -3, 3, -3, 3)
names(prior.reg.params)<-c("b0.min", "b0.max", "b1.min", "b1.max", "b2.min", "b2.max")
ououbmprior(prior.model.params=prior.model.params,prior.reg.params=prior.reg.params)

ouxy documentation built on July 2, 2020, 4:05 a.m.