# NOT RUN {
# create a GeneralSubstitution object
# provide an Alphabet object and the rates
p<-GeneralSubstitution(alphabet=BinaryAlphabet(), rate.list=list("1->0"=1,"0->1"=1))
# get equlibrium distribution
getEquDist(p)
# get equilibrium distribution via virtual field
p$equDist
# sample from equilibrium distribution
sampleState(p)
# }
Run the code above in your browser using DataLab