# simple estimation on iris data, clustering by "Species"
res_poisson = femlm(Sepal.Length ~ Sepal.Width + Petal.Length +
Petal.Width | Species, iris)
# we plot the residuals
plot(resid(res_poisson))
Run the code above in your browser using DataLab