Создание участка с ошибкой/кумулятивным инцидентом/событием с svykm

Question

Я использовал пакет twang для оценки веса с учетом моих данных и использую функции пакета опроса.

Я могу создать стандартную кривую выживания без проблем, используя следующий код:

km <- svykm(Surv(rec_time,rec)~procedure, design=design.mnps) 
plot(km) 

Я посмотрел в Интернете, и если он не был объектом svykmlist, следующий будет работать:

km <- svykm(Surv(rec_time,rec)~procedure, design=design.mnps, function(x) 1-x) 
plot(km) 

я получаю следующее сообщение об ошибке:

1: In plot.window(...) : "FUN" is not a graphical parameter 

и т.д.

Я также пробовал FUN = «событие», которое отлично работает с традиционным объектом съемки, но не с объектом svykm. Я также попытался манипулировать объектом съемки, который не увенчался успехом.

Есть ли что-нибудь, что я еще не пробовал, что может работать с объектом svykmlist?

Спасибо, Matt

Answer 1

, когда вы размещаете на SO, пожалуйста, используйте минимальный воспроизводимый пример. объект design.mnps не пригодится. отредактируйте свой вопрос и уточните, почему вы не просто инвертируете свою функцию в формуле, как это? благодаря

пример из? svykm

library(survey) 
data(pbc, package="survival") 
pbc$randomized <- with(pbc, !is.na(trt) & trt>0) 
biasmodel<-glm(randomized~age*edema,data=pbc) 
pbc$randprob<-fitted(biasmodel) 
dpbc<-svydesign(id=~1, prob=~randprob, strata=~edema, data=subset(pbc,randomized)) 

s2<-svykm(Surv(time,status>0)~I(bili>6), design=dpbc) 
plot(s2) 

sz<-svykm(Surv(time,status>0)~I(1 - I(bili>6)), design=dpbc) 
plot(sz)