11 Painotusten laskeminen
Tässä vaiheessa päästään vihdoinkin laskemaan painotukset.
Painotukset lasketaan komennolla ipcw:
toy.rep <- ipcw(toy.rep, id = "id", tstart = tstart, tstop = tstop,
cens = cens, arm = "arm", bas.cov = c("age"),
conf = c("TDconf"), trunc = 0.01, type = 'kaplan-meier')Tämä analyysi heittää warningseja runsaasti, mutta ajoin oleellisesti ihan samat analyysit R:n helppidokumentissa käytetyllä esimerkkidatasetillä (data(“SHIdat”)) ja samassa dokumentissa listatulla koodilla läpi, ja aivan samantyyppiset warningit näyttivät silläkin tulevan.
Saadaan aikaiseksi seuraavalta näyttävä datasetti:
| id | age | arm | swtrtdt | TDconf | cens | weights | tstart | tstop | event | weights.trunc |
|---|---|---|---|---|---|---|---|---|---|---|
| 1 | 20 | A | NA | 1 | 0 | 1 | 0 | 41 | 0 | 1 |
| 1 | 20 | A | NA | 1 | 0 | 1 | 41 | 49 | 1 | 1 |
| 3 | 40 | A | 226 | 1 | 0 | 1 | 0 | 41 | 0 | 1 |
| 3 | 40 | A | 226 | 1 | 0 | 1 | 41 | 49 | 0 | 1 |
| 3 | 40 | A | 226 | 1 | 0 | 1 | 49 | 154 | 0 | 1 |
| 3 | 40 | A | 226 | 1 | 0 | 1 | 154 | 155 | 0 | 1 |