s-news
[Top] [All Lists]

Hmisc multiple imputation with owner glm family

To: <s-news@lists.biostat.wustl.edu>
Subject: Hmisc multiple imputation with owner glm family
From: "GIORGI Roch" <roch.giorgi@mail.ap-hm.fr>
Date: Thu, 31 May 2007 19:12:30 +0200
Thread-index: Acejpt7q7VMlHkiFQ6ixuj6GegMeXg==
Thread-topic: Hmisc multiple imputation with owner glm family
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/


<Prev in Thread] Current Thread [Next in Thread>
  • Hmisc multiple imputation with owner glm family, GIORGI Roch <=