| To: | <s-news@lists.biostat.wustl.edu> |
|---|---|
| Subject: | ? Error in stepAIC?? |
| From: | "Hunsicker, Lawrence" <lawrence-hunsicker@uiowa.edu> |
| Date: | Mon, 29 Dec 2008 10:15:29 -0600 |
| Thread-index: | Aclp0KrTdkMmx2cdQJ2J90tVih4GnQ== |
| Thread-topic: | [S] ? Error in stepAIC?? |
|
Good morning, all:
I have tried to run stepAIC on a glmmPQL model. The first function (glmmPQL) works fine. But when I run the stepAIC function, I get the error: ‘Object “tt” not found’ ScrMiss.glmm<-glmmPQL(fixed = Scr.missing ~ Sexe + Age1 + Age.miss + Ethnie + HistoryDiabete + DurDiab + DurDiab.miss + Retinopathie + PiedDiabetique + Hypertension + DurBP + DurBP.miss + SBP1 + SBP.miss + DBP1 + DBP.miss + AntecedentTabac + Hyperlipidemie + LVH + CAD + MI + CHF + Stroke + PVD + MedicamentDiabete + NoBPMeds + Diuretique + ACEInhibitor + Calciumchannelblockers + ARBs + Statins + OtherLipid + NoAnticoag + Aspirin + OtherAnticoagulant + NumberBPClasses , random = ~1|Region/Pays/Center, family = binomial(link=logit),data="" = F, niter=15) minmod.glmm <- stepAIC(ScrMiss.glmm, direction = "both") > minmod.glmm <- stepAIC(ScrMiss.glmm, direction = "both") Start: AIC= 126349.7 Scr.missing ~ Sexe + Age1 + Age.miss + Ethnie + HistoryDiabete + DurDiab + DurDiab.miss + Retinopathie + PiedDiabetique + Hypertension + DurBP + DurBP.miss + SBP1 + SBP.miss + DBP1 + DBP.miss + AntecedentTabac + Hyperlipidemie + LVH + CAD + MI + CHF + Stroke + PVD + MedicamentDiabete + NoBPMeds + Diuretique + ACEInhibitor + Calciumchannelblockers + ARBs + Statins + OtherLipid + NoAnticoag + Aspirin + OtherAnticoagulant + NumberBPClasses Problem in switch(data.class(object),: Object "tt" not found Use traceback() to see the call stack I checked to see whether the glmmPQL object (ScrMiss.glmm) had the “extract.AIC()” method, and it does. What am I doing wrong this time? As always, my thanks to anyone that can get me on the right track again. Larry Hunsicker
|
| <Prev in Thread] | Current Thread | [Next in Thread> |
|---|---|---|
| ||
| Previous by Date: | Re: Writing a function taking a data.frame column name as a paramete, Hunsicker, Lawrence |
|---|---|
| Next by Date: | [no subject], yiwu ye |
| Previous by Thread: | Re: Writing a function taking a data.frame column name as a paramete, Hunsicker, Lawrence |
| Next by Thread: | [no subject], yiwu ye |
| Indexes: | [Date] [Thread] [Top] [All Lists] |