Значение прогноза ГЛМ в р

Question

Мои данные включают: Год (2000-2010), Квартал (Q1-Q4), Площадь (A1-A5), Усилие (количество крючков), Улов (количество рыбы). Я использовал GLM, чтобы предсказать значение улова. Вот мой код:

 `model = glm( log(catch) ~ Year + Quarter + Area + log(Effort), data = 
  data, family = gussian())`

И я хочу получить значение прогноза при любых условиях в 3 квартале и на 1000 крючков! так что я написал:

`X0 = model.martrix(model)
 effort.pred=1000
 X0[,”log(Effort)”] <- log(effort.pred)
 X0[,”QuarterQ2”] <- 0
 X0[,”QuarterQ3”] <- 1
 X0[,”QuarterQ4”] <- 0
 X0[,”AreaA2”] <- 1
 X0[,”AreaA3”] <- 1
 X0[,”AreaA4”] <- 1
 X0[,”AreaA5”] <- 1
 pred = X0%*%martrix(coef(model))`

Но я не могу получить значение ... Я просто хочу получить столик

 Year    Quarter    Area     pred
 2000      Q3        A1        ?
 2000      Q3        A2        ?
 2000      Q3        A3        ?
 2000      Q3        A4        ?
 2000      Q3        A5        ?
 2001      Q3        A1        ?
                .
                .
                .
                .
 2010      Q3        A3        ?
 2010      Q3        A4        ?
 2010      Q3        A5        ?

Кто-нибудь может мне помочь, пожалуйста !!!!