I can't confirm without a reproducible example but I'm fairly sure the problem is that status
is considered a factor in your dataframe.
You can put
dfit <- survfit2(Surv(time, as.numeric(status)) ~ intervention, data = Drugs_discon)
summary(dfit)
ggsurvfit(dfit)
This should produce the required plot