Displaying 2 results from an estimated 2 matches for "hakulinen".
2011 Feb 15
0
Delta method using numerical derivatives
...in SAS).
I would like something akin to the following predictnl() function:
## get some data and fit a model
require(Epi)
data(lungDK)
fit <- glm(D/Y ~ as.factor(A5)+as.factor(A5):P5-1,
data=lungDK, subset=(P5>=1970), weight=Y,
family=poisson("identity")) # (Hakulinen and Dyba)
## calculate some summary statistics (e.g. cumulative incidence to 85
years)
require(plyr)
cuminc <- function(object,year=P5,newdata=NULL)
ddply(newdata,deparse(substitute(year)),function(data)
sum(5*predict(object,newdata=data)))[,-1,drop=TRUE]
## Now: get summary statistics...
2008 Jan 29
2
Direct adjusted survival?
Hello,
I am trying to find an R function to compute 'direct adjusted survival'
with standard errors. A SAS-macro to do this is presented in Zhang X,
Loberiza FR, Klein JP, Zhang MJ. A SAS macro for estimation of direct
adjusted survival curves based on a stratified Cox regression model.
Comput Methods Programs Biomed 2007;88:95-101. It appears that this
method is not implemented in R.