I want to show the change of a numerical score in 68 indivials at two time points, before and after a treatment.
I' struggling with the user written command profileplot. It generates a plot with different colors per observation.

profileplot var1 var2 , by(id) msymbol(i) xlabel(1 "Before" 2 "After")

But, I want three colors only: Green for those who decrease, red for those who increase and blue for those with no change. I still want each individual observation presented as a line (i.e. not a mean per group). I've generated a variable for this purpose ("Change") which tells which group the individual belongs to (1,2 or 3)
Does anyone know if this is possible, or, if I should use another command?