Как использовать определенный факторный уровень mrds: функция обнаружения ddf

Question

Мне нужно рассчитать плотность нескольких видов птиц в районе исследования по данным дистанционной выборки (назовем это "AllObservations"). Я стремлюсь: 1. использовать общую функцию обнаружения для всех видов птиц, используя mrds:ddf с видами птиц в качестве фактора

MyDDF=ddf(method='ds',dsmodel=~mcds(key='hr',formula=~1+factor(species)),data=AllObservations, meta.data=list())

2. подмножество данных, относящихся только к одному виду (назовем это "OneSpObservations"), а потом

3. используйте dsm::dsm, чтобы соответствовать модели поверхности плотности

   dsm(ddf.obj=chosenDdf, formula = D~lon+lat, family=gaussian(link='identity'), group=FALSE,engine='glm',convert.units=1,segment.data=OneSpSegments,observation.data=OneSpObservations)

Тем не менее, как на самом деле использовать фактор (виды) ДДФ в следующей части? Если я использую то же самое OneSpObservations и только меняя название вида, я получаю идентичные результаты.