Skip to content
/ Rwtd Public

Rough R translation of functions from the Stata wtdttt package

License

Notifications You must be signed in to change notification settings

mbg-unsw/Rwtd

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 

History

9 Commits
 
 
 
 
 
 
 
 

Repository files navigation

Rwtd

n.b. this code has been superseded by the Rwtdttt package

Rough R translation of functions from the Stata wtdttt package

(c) 2023 Malcolm Gillies malcolm.gillies@unsw.edu.au

This code is a direct translation of portions from the Stata module wtdttt, available at: https://econpapers.repec.org/software/bocbocode/s458265.htm by Katrine Bødkergaard Nielsen (kani@ph.au.dk) and Henrik Støvring (stovring@ph.au.dk)

These functions estimate parameters of the ordinary and reverse Waiting Time Distribution (WTD) by maximum likelihood, with application to modelling the recurrence of drug dispensing. Please see the Stata package documentation for references.

As with wtdttt, this code is made available under the terms of the GPL v3 (or later).

Currently, only the lognormal distribution is implemented.

Example data in Stata format, taken directly from the wtdttt package, is included for demonstration purposes

library(haven)

x <- read_dta("ref/wtddat_dates.dta")

x.w <- wtdttt(x$rx1time, as.Date('2014-01-01'), as.Date('2014-12-31'),
      reverse=F)

summary(x.w)

x <- read_dta("ref/ranwtddat_discdates.dta")

x.w <- ranwtdttt(x$rxdate, x$pid, as.Date('2014-01-01'),
      as.Date('2014-12-31'), reverse=F)

summary(x.w)

x.w <- ranwtdttt(x$rxdate, x$pid, as.Date('2014-01-01'), as.Date('2014-12-31'),
    reverse=T)

summary(x.w)

About

Rough R translation of functions from the Stata wtdttt package

Topics

Resources

License

Stars

Watchers

Forks

Languages