Therneau, Terry M., Ph.D. mayo.edu> writes:
> This issue has arisen within my anova.coxph routine, but is as
> easily illustrated with glm.
> testdata <- data.frame(y= 1:5,
> n= c(8,10,6,20,14),
> sex = c(0,1,0,1,1),
> age
On Thu, 23 Apr 2015, Therneau, Terry M., Ph.D. wrote:
This issue has arisen within my anova.coxph routine, but is as easily
illustrated with glm.
testdata <- data.frame(y= 1:5,
n= c(8,10,6,20,14),
sex = c(0,1,0,1,1),
age = c(30,
> "Save the model frame in case you need to refit something next month"
> does not appear to be a safe approach to reproducible research.
Is this a standard recommendation? It will not work in many cases. E.g.,
if
you use lm() to model the sum of some variables the model.frame contains
only the
This issue has arisen within my anova.coxph routine, but is as easily
illustrated with glm.
testdata <- data.frame(y= 1:5,
n= c(8,10,6,20,14),
sex = c(0,1,0,1,1),
age = c(30,20,35,25,40))
fit <- glm(cbind(y,n) ~ age + sex, bin
