Nothing
`pp.rho` <-
function(mu,sigma=.553,
rho= c(0, .25, .5, .75, 1),
COLORS=c("black","blue","green","red","black"),nsim=10^5){
nmu<-length(mu)
nr<-length(rho)
out1<-rep(NA,nmu)
for (i in 1:nmu){
out1[i]<-hbpp(mu[i],sigma^2)
}
out2<-out3<-matrix(rep(NA,nmu*nr),nmu,nr)
for (i in 1:nmu){
for (j in 1:nr){
V<-make.v(2,rho[j],sigma^2)
out2[i,j]<-hbpp(c(mu[i],mu[i]),V,simulate=TRUE,nsim=nsim)
}
}
for (i in 1:nmu){
for (j in 1:nr){
V<-make.v(3,rho[j],sigma^2)
out3[i,j]<-hbpp(c(mu[i],mu[i],mu[i]),V,simulate=TRUE,nsim=nsim)
}
}
out<-list(mu=mu,
out1=out1,out2=out2,out3=out3,col1=c("black"),col2=COLORS,col3=COLORS,cparms=rho)
return(out)
}
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.