广义估计方程

请问下将表格宽变长,用R语言GEE评估一组数据重复测量的差别及两两比较。分为4列:id、time(1、2、3、4)、variable(密度)和value(密度的值),就是比较4次密度的变化,用geeglm函数 fit<—geeglm(value~variable+time,id=id,data=…)总报错:contrasts can be applied only to factors with 2 or more levels。请问是不是variable设置有误,我variable那列都写“密度”这两个字