求助：glm函数能否进行3分类或者更多分类结局变量的线性回归模型？

例如，iris数据集
取前100行，species是二分类结局，这样即可：
glm(Species~.,data = iris[1:100,],family = 'binomial')
如果是取iris整个数据集，species作为结局有三个分类，glm函数中family这个参数应该写什么呢？
glm(Species~.,data = iris, family = '???')