各位大侠新年好：
      我有如下数据，分析过程中发现了些问题，详述如下：

病人编号pid	是否使用药物d_or	使用药物次数d_time	第几次随访 sf_time
01	否	0	1
01	否	0	3
01	是	12	13
01	是	89	23
02	是	67	6
02	否	0	21
02	是	.	22
03	是	602	4
03	是	11	21
04	是	34	3
04	是	.	4
04	否	0	5
04	否	0	6
05	否	0	6
05	否	.	21

目的：了解病人使用药物次数 或者 是否使用药物随时间变化的趋势图（或者叫做重复测量资料概貌分析）
提示：原始数据pid有6000多例； 是否使用药物（d_or)变量有部分缺失； 变量d_time最小值0，最大值902，且有大量缺失；
          变量 sf_time最小值1，最大值26，因随访时间不等，重新整理的结果。
例出程序之一如下：proc gplot data=a1;
plot d_time*sf_time=pid/vaxis=axis1 haxis=axis2;
axis1 length=80.0 label=(a=90 h=0.8 'value of inject_times');
axis2 length=80.0 offset=(2) order=(0 to 26 by 1)
label=(h=0.8 'sf_times');
run;
出现问题：图乱码，请问程序错在哪？或者为达到此次目的，更好的建议？