如何将prodlim函数与公式中的非二进制变量一起使用？

我试图（最终）使用prodlim函数按组绘制数据。 我正在调整和改编其他人（没有问题）编写的代码，我对prodlim库/函数不是很熟悉。当然还有其他方法可以做我想做的事情，但我正努力让它与前一个人所做的保持一致。
当我将数据分为两组时，我的代码是有效的，但当我尝试针对4组情况进行调整时，我得到了一个错误。
值得注意的是，数据是使用StatTransfer从SAS传来的，StatTransfer一直运行良好

我不熟悉编码，但我已经比较了我尝试使用的数据帧。第二个变量只是第一个变量的子集（代码在其中工作），具有所有相同的变量，并且我试图根据这两个变量分组的变量都是整数值。

Hist（medpop$dz_time，medpop$dz_status）

工作正常，因此问题一定出在prodlim函数上，遗憾的是，我还不太了解我查找到的关于它的内容：/但是文档似乎表明它支持连续或分类变量，而且似乎也不限于二进制。据我所知，这些选项似乎都不适用

这项工作：

M <- prodlim(Hist(dz_time, dz_status)~med, data=pop)

M
N <- prodlim(Hist(dz_time, dz_status)~strength, data=medpop)
N <- prodlim(Hist(dz_time, dz_status)~strength, data=pop, subset=pop$med==1)