Changing colour in scatter plot

[Martin Imelda Borg]

Hi I'm plotting deviance residuals following a survival analysis

I would like to change the colour for those that have suffered a failure vs those who haven't, similar to this

Up to now I have managed to do the following

stcox var1 var2 var3
predict mg, mgale
predict xb 
predict dev, deviance
scatter dev xb || lfit dev _t

However I don't know how to insert a colour option for those who have failed / those who have been censored...

[Maarten Buis]

sysuse cancer, clear
stcox i.drug age
predict dev, deviance
separate dev, by(_d)
label var dev0 "censored"
label var dev1 "died"
scatter dev0 dev1 _t || lfit dev _t,  ///
  legend(order(1 2)) ytitle(deviance residual)

[Martin Imelda Borg]

Thanks for this
However with the code you supplied in post 2 I get this graph



However when I change the code to:

 
 stcox var1 var2 var3 predict dev, deviance separate dev, by(_d)  label var dev0 "censored"  label var dev1 "died" predict xb  scatter dev0 dev1 xb _t || lfit dev _t,  ///   legend(order(1 2)) ytitle(deviance residual)

I get the following graph:


Where blue = censore, red = died

Why don't I get the scatter dots with your code? The line of best fit with your code seems steeper compared to my code I edited her...

Also why are mine seem to be clustered all in the region of <5?
Not sure if this would be in your remit but,
I also seem to get a downward trend line, but shouldn't the line be at 0 indicating a good fit as per Colletts book?

[Martin Imelda Borg]

Maarten Buis appreciate your thoughts on this …