Nothing
rMSTT <-
function(n, mu, Sigma, lambda, nu=1){
if(is.null(n))
stop("sample size must be specified")
if(is.null(mu))
stop("mu must be specified")
if(is.null(Sigma))
stop("Sigma must be specified")
if(is.vector(Sigma)) Sigma=matrix(Sigma)
if(is.null(lambda))
stop("lambda must be specified")
if(round(n)!=n | n<=0)
stop("sample size must be a positive integer")
if(!is.matrix(mu)) mu<-matrix(mu, ncol=length(mu),nrow=1)
if (!isSymmetric(Sigma))
stop("Sigma is not symmetric")
if (min(eigen(Sigma)$values)<=0)
stop("Sigma is not a positive-definite matrix")
if (ncol(Sigma) != ncol(mu))
stop("The dimension of mu does not agree with the dimension of Sigma")
if (ncol(Sigma) != length(lambda))
stop("The dimension of lambda does not agree with the dimension of Sigma")
if (ncol(mu) != length(lambda))
stop("The dimension of mu does not agree with the dimension of lambda")
if(nu<=0)
stop("nu must be positive")
p<-ncol(Sigma)
lambda=as.vector(lambda)
delta <- lambda/c(sqrt(1 + t(lambda)%*%lambda))
y<- matrix(0,n,p)
for (i in 1:n){
u=rgamma(1,nu/2,nu/2)
Z<-rMSN(1, mu=rep(0,p), Sigma=diag(p), lambda=lambda)
y[i,]=c(mu)+u^(-1/2)*matrix.sqrt(Sigma)%*%t(Z)
}
return(y)
}
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.