data("patent")
patentMix <- stepFlexmix(Patents ~ lgRD, k=3,
model=FLXMRglm(family="poisson"),
nrep=5, data=patent)
plot(Patents ~ lgRD, data=patent,
pch=as.character(cluster(patentMix)))
ordering <- order(patent$lgRD)
apply(fitted(patentMix), 2, function(y)
lines(sort(patent$lgRD), y[ordering]))
Run the code above in your browser using DataLab