I note that you have some non-estimable main effect LSmeans, so I am going to suggest the following, assuming the model is fit in PROC MIXED:
proc mixed data=yourdataset;
class treatment time subjectid;
model response_variable=treatment*time;
repeated time/subject=subjectid type=sp(pow)(t);
lsmeans treatment*time;
slice treatment*time / sliceby=time diff;
run;
Some things to note here. I see that your time points are not equally spaced. In yourdataset, define t=time, so that you will have a continous version of time to feed into the spatial power covariance structure. The key here is the SLICE statement. This should enable you to get differences between lsmeans for treatment at each time point where data exists for the treatments. At some points, only one treatment was measured, so no difference can be calculated.
Steve Denham
... View more