R/seqgen.R

In TraMineR: Trajectory Miner: a Sequence Analysis Toolkit

## ============================
## Random sequences generation
## ============================

seqgen <- function(n, length, alphabet, p=NULL) {

	## FROM http://www.cs.chalmers.se/~dubhashi/Courses/BioAlg/Markov%20chains%20and%20DNA%20sequence%20analysis%20in%20R.htm

	m <- matrix(nrow=n, ncol=length)

	for (i in 1:n) 
		m[i,] <- sample(alphabet, length, replace=TRUE, p)

	colnames(m) <- paste("[",1:length,"]", sep="")
	rownames(m) <- 1:n

	return(m)
}