DS.prior {BayesGOF} | R Documentation |
Prior Diagnostics and Estimation
Description
A function that generates the uncertainty diagnostic function (U-function
) and estimates DS(G,m)
prior model.
Usage
DS.prior(input, max.m = 8, g.par,
family = c("Normal","Binomial", "Poisson"),
LP.type = c("L2", "MaxEnt"),
smooth.crit = "BIC", iters = 200, B = 1000,
max.theta = NULL)
Arguments
input |
For |
max.m |
The truncation point |
g.par |
Vector with estimated parameters for specified conjugate prior distribution |
family |
The distribution of |
LP.type |
User selects either |
smooth.crit |
User selects either |
iters |
Integer value that gives the maximum number of iterations allowed for convergence; default is 200. |
B |
Integer value for number of grid points used for distribution output; default is 1000. |
max.theta |
For |
Details
Function can take m=0
and will return the Bayes estimate with given starting parameters. Returns an object of class DS.GF.obj
; this object can be used with plot command to plot the U-function (Ufunc
), Deviance Plots (mDev
), and DS-G comparison (DS_G
).
Value
LP.par |
|
g.par |
Parameters for |
LP.max.uns |
Vector of all LP-Fourier coefficients prior to smoothing, where the length is the same as |
LP.max.smt |
Vector of all smoothed LP-Fourier coefficients, where the length is the same as |
prior.fit |
Fitted values for the estimated prior. |
UF.data |
Dataframe that contains values required for plotting the U-function. |
dev.df |
Dataframe that contains deviance values for values of |
m.val |
The value of |
sm.crit |
Smoothing criteria; either |
fam |
The user-selected family. |
LP.type |
User-selected representation of |
obs.data |
Observed data provided by user for |
Author(s)
Doug Fletcher, Subhadeep Mukhopadhyay
References
Mukhopadhyay, S. and Fletcher, D., 2018. "Generalized Empirical Bayes via Frequentist Goodness of Fit," Nature Scientific Reports, 8(1), p.9983, https://www.nature.com/articles/s41598-018-28130-5.
Mukhopadhyay, S., 2017. "Large-Scale Mode Identification and Data-Driven Sciences," Electronic Journal of Statistics, 11(1), pp.215-240.
Examples
data(rat)
rat.start <- gMLE.bb(rat$y, rat$n)$estimate
rat.ds <- DS.prior(rat, max.m = 4, rat.start, family = "Binomial")
rat.ds
plot(rat.ds, plot.type = "Ufunc")
plot(rat.ds, plot.type = "DSg")
plot(rat.ds, plot.type = "mDev")