# NOT RUN {
# simple estimation on iris data, using "Species" fixed-effects
res = femlm(Sepal.Length ~ Sepal.Width*Petal.Length +
Petal.Width | Species, iris)
head(model.matrix(res))
# }
Run the code above in your browser using DataLab