Nothing
BMDU <-
function(x,d)
{
k=length(x)
for(m in 1:k)
{
if(x[m]<0)
{cat('Warning: arguments must be non-negative integers')
}
}
if(d>=1 && d<=k)
{
m=0
l=0
u=0
diff=0
s=sum(x)
s1=floor(k/d)
d1=runif(s1,0,1)###First half of the vector
d2=runif(k-s1,1,2)###Second half of the vector
a=c(d1,d2)
p=x+a###Prior for Dirichlet
dr=rdirichlet(10000, p)###Posterior
for(j in 1:k)
{
l[j]=round(quantile(dr[,j],0.025),4)###Lower Limit
u[j]=round(quantile(dr[,j],0.975),4)###Upper Limit
m[j]=round(mean(dr[,j]),4)###Point Estimate
diff[j]=u[j]-l[j]
}
cat('Mean\n')
print(m)
cat('Lower Limit\n')
print(l)
cat('Upper Limit\n')
print(u)
p=prod(diff)
cat('Volume\n')
print(p)
}
else
{cat('warning:size of the division should be less than the size of the input matrix')
}
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.