search for: bl_cholest

...names")= chr [1:886393] "1" "2" "3" "4" ... $ means : num [1:16] 43.472 4.162 0.563 52.815 8.503 ... $ terms :Classes 'terms', 'formula' length 3 srv900 ~ rcs(Age, c(35, 50, 65)) + Sex + rcs(BL_HDL.A, 3) * rcs(BL_CHOLEST.A, 3) + rcs(BL_GGT.A, 3) + BL_TRIGLYC.A + rcs(BL_GGT.A, 3) + ... .. ..- attr(*, "variables")= language list(srv900, rcs(Age, c(35, 50, 65)), Sex, rcs(BL_HDL.A, 3), rcs(BL_CHOLEST.A, 3), rcs(BL_GGT.A, 3), BL_TRIGLYC.A, BP.B, ... .. ..- attr(*, "factors")= in...

Hi, I have a list of fit objects (fit objects from HMISC functions) I create elements in the list in this way lrm.sumtot <- lrm( ae7bepn ~ trarm + sumtot , data=sd.fix) lrm.list[['lrm.sumtot']] <- lrm.sumtot And I can run (anova(lrm.sumtot)) The following also gives the anova I'd expect zz <- lrm.list[['lrm.sumtot']];anova(zz) And similarly for the summary

...took a crack at hacking the plot.Design function, adding a col=col parameter to be passed in the function call to image(), but failed to get the desired effect: # else image(xseqn, y, zmat, xlab = xlab, ylab = laby , col = col) library(Hmisc); library(Design) lr.fit6 <- lrm(death ~ rcs(BL_CHOLEST.A,c(180,220,280))*rcs(BL_HDL.A,c(40,55,70))*Sex, data = pref900) # str(pref900[,c("BL_HDL.A","BL_CHOLEST.A","death")]) $'data.frame': 910659 obs. of 3 variables: $ BL_HDL.A : num 34 35 40 46 39 45 46 34 42 52 ... $ BL_CHOLEST.A: num 181 184 238 134...