Generates the suite of functions related to the one sample normal experiment with a onesided alternative hypothesis of interest when the variance is known.
1  norm1KV.1sided(sigma, theta0, mu, tau)

sigma 
Scalar. The known standard deviation of the population of interest. 
theta0 
Scalar. The critical value of the mean under the null hypothesis: theta<=theta0. 
mu 
Scalar. The mean of the normal prior density on theta. See documentation
for 
tau 
Scalar. The standard deviation for the normal prior density on theta. See
documentation for 
norm1KV.1sided
is used to generate a suite of functions for a
onesample normal experiment with a onesided alternative hypothesis when the
variance is known. That is, when
X ~ Normal(theta,sigma2)
H0: theta <= theta0 vs. H1: theta > theta0
using the following prior on theta
pi(theta) = Normal(mu,tau2),
where Normal(mu,tau2) is Normal density with mean mu
and variance
tau2
.
The functions that are generated are useful in examining the prior and
posterior densities of the parameter theta
, as well as constructing the
Bayes Factor and determining the sample size via an average error based
approach.
The arguments of norm1KV.1sided
are passed to each of the additional
functions upon their creation as default values. That is, if mu
is
set to 1 in the call to norm1KV.1sided
, each of the functions returned
will have the defaualt value of 1 for mu
. If an argument is not
specified in the call to norm1KV.1sided
, then it remains a required
parameter in all functions created.
norm1KV.1sided
returns a list of 4 functions:
logm 
Returns a list of three vectors: the log marginal density under
the null hypothesis ( logm(xbar, n, sigma, theta0, mu, tau)

logbf 
Returns a vector. The value of the log Bayes Factor given the observed data provided and the prior parameters specified. The function has the following usage: logbf(xbar, n, sigma, theta0, mu, tau) See 
prior 
Returns a vector. The value of the prior density. The function has the following usage: prior(theta, mu, tau)

post 
Returns a vector. The value of the posterior density. The function has the following usage: post(theta, xbar, n, sigma, mu, tau)

binom1.1sided
,binom1.2sided
,
binom2.1sided
,binom2.2sided
,
norm1KV.2sided
,norm2KV.2sided
norm1UV.2sided
,ssd
,BAEssd
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30  ############################################################
# Generate the suite of functions for a onesample normal
# with a onesided test. Consider the hypothesis
# H0: theta<=0 vs. H1: theta>0
#
# with a normal prior for theta with prior mean 2 and
# prior standard deviation 1. The known standard
# deviation for the population is 5.
# generate suite
f5 < norm1KV.1sided(sigma=5,theta0=0,mu=2,tau=1)
# attach suite
attach(f5)
# plot the prior and posterior density when the observed data is
# n = 30, xbar = 1.
mus < seq(5,5,0.01)
mu1 < prior(mus)
mu2 < post(mus,xbar=1,n=30)
plot(c(mu1,mu2)~rep(mus,2),type="n",ylab="Density",xlab="mu",main="")
lines(mu1~mus,lty=1,lwd=2)
lines(mu2~mus,lty=2,lwd=2)
# perform sample size calculation with TE bound of 0.25 and weight 0.5
ssd.norm1KV(alpha=0.25,w=0.5,logm=logm)
# detain suite
detach(f5)

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.
Please suggest features or report bugs with the GitHub issue tracker.
All documentation is copyright its authors; we didn't write any of that.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.