R/r.FGSTN.R

In a-mahdavi/SSMFGSN.EM: EM algorithms for the SSMFGSN distributions

r.FGSTN <- function(n , xi, s, la, nu){
		m <- length(la) ; a=seq(1,2*m-1,by=2)
		y <- 0
		for(i in 1:n){
		x0 <- rt(1,nu)
		a.x0 <- t(outer(x0,a,'^'))
		x1<- rnorm(1)
		if( x1 < as.numeric(t(la)%*%a.x0) )
		y[i] <- x0
		else
		y[i] <- -x0
		}
		return(xi + s*y)
		}