预测概率

Question

我需要在R中使用predict（）函数来尝试获取一组值的概率。数据集如下：

Loc: Beach, NonBeach
Sex: Male, Female

Infec: (Response variable) 

带有glm：

mod1 = glm(Infec ~ Loc + Sex, family = binomial)

我想确定所有这些变量的概率值，即使用预测的是男性-海滩，男性-非海滩，女性-海滩，女性-非海滩。我知道我可以使用以下方式的预测函数：

predict(mod1, newdata = list(Loc = "Beach", Sex = "Female"), type = "response") 

获得“海滩和女性”的答案，但是我不确定如何在相同结果中获得所有四个概率，例如：

Female    Beach     xxx
Male      Beach     xxx
Female    NonBeach  xxx
Male      NonBeach  xxx

没有打印出与观察值一样多的值。