Dear all,
I'm using Hmisc package to do multiple imputation.
I have to use glm method but with a family that is not yet part of S-PLUS. I
created my own family but I don't know how to call it when I use Hmisc's
fit.mult.impute.
I tried this:
mi <- aregImpute(~Y + sex + age1 + age2 + X1, n.impute=5, data=MCAR)
fit <- fit.mult.impute(Y ~ sex + age1 + age2 + X1, family="MyFamily",
fitter=glm, xtrans=mi, data=MCAR)
but I have the impression that the family used is not MyFamily
fit$family gives:
name link variance
"Gaussian" "Identity: mu" "Constant: 1"
Any help would be much appreciated. I am using Splus6.2 and Hmisc 3.3-2.
Thanks
Roch
_____________________________________________________
Roch GIORGI, MD, PhD
Service de Santé Publique et d'Information Médicale
Hôpital de la Timone
264, rue St Pierre
13385 Marseille cedex 5
Tel: +33 (0)491 384 949
Fax: +33 (0)491 385 749
email: roch.giorgi@ap-hm.fr
Web: http://cybertim.timone.univ-mrs.fr/
|