#Simulated data
simDat <- skipTrack.simulate(n = 100, skipProb = c(.7, .2, .1))
#Run model fit (should typically run with much more than 50 reps)
modFit <- skipTrack.fit(Y = simDat$Y, cluster = simDat$cluster, chains = 2, reps = 50)
summary(modFit, burnIn = 25) #Recommended burnIn with real data is at least 750
Run the code above in your browser using DataLab