Displaying 1 result from an estimated 1 matches for "lrm_reduc".
Did you mean:
lrm_reduced
2005 Jul 12
1
Design: predict.lrm does not recognise lrm.fit object
...he results with an AIC/MMI approach. The script contains:
# Determine full logistic regression
lrm_logist = lrm(PresAbs ~ Size + X2ndpc + soil + AAR + tjan.jun,
data=training)
# Backward selection of variables in model
lrm_stp = fastbw(lrm_logist, rule="p", sls=0.05)
# Fit reduced model
lrm_reduced = lrm.fit(training[,lrm_stp$parms.kept[-1]], training$PresAbs)
# Predict using parameters from reduced model against a new dataset
predict(lrm_reduced, testing, type="fitted.ind")
It is the last command that fails, reporting the error:
Error in getOldDesign(fit): fit was not created wi...