Sorry. I've used: > library(rms) I realise I still have a lot to learn to ask questions well - it took me a long time to compile this one, but I've obviously missed important things. Please see below for the session info.
> sessionInfo() R version 2.15.0 (2012-03-30) Platform: i386-pc-mingw32/i386 (32-bit) locale: [1] LC_COLLATE=English_United Kingdom.1252 LC_CTYPE=English_United Kingdom.1252 [3] LC_MONETARY=English_United Kingdom.1252 LC_NUMERIC=C [5] LC_TIME=English_United Kingdom.1252 attached base packages: [1] splines stats graphics grDevices utils datasets methods base other attached packages: [1] rms_3.5-0 Hmisc_3.9-3 survival_2.36-12 loaded via a namespace (and not attached): [1] cluster_1.14.2 grid_2.15.0 lattice_0.20-6 tools_2.15.0 Many thanks for your help. Cecile On 2 August 2012 01:00, R. Michael Weylandt <michael.weyla...@gmail.com>wrote: > What package(s) are the functions in question from? > > This might also help: > > > http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example > > Michael > > On Wed, Aug 1, 2012 at 2:57 AM, Cecile De Cat <c.de...@leeds.ac.uk> wrote: > > You're right, it's just the 2 columns that are characters that return > > false. But I don't use them in the analysis (it's the experiments' > > names and the participants' names). > > > > So I guess I'm back to my original question (although I can discard > > one possible cause thanks to you): there appears to be only "real > > numbers" in the data used for the lrm analysis, and yet it falls over. > > > > Thanks a lot for your help. > > > > Cecile > > > > > > On 31 July 2012 16:42, R. Michael Weylandt <michael.weyla...@gmail.com> > wrote: > >> What classes are the columns of your data frame? > >> > >> Note that > >> > >> is.finite("a") # False > >> is.finite(factor("a")) # True > >> > >> M > >> > >> On Tue, Jul 31, 2012 at 10:34 AM, Cecile De Cat <c.de...@leeds.ac.uk> > wrote: > >>> Thank you. This is very useful. I do indeed get the following: > >>>> table(sapply(dat, is.finite)) > >>> FALSE TRUE > >>> 28164 253476 > >>> > >>> But the number of observations returned baffles me, as there should > >>> only be 14082 in the data. And when I look at each variable > >>> individually, none appear to violate "is.finite": e.g. > >>> > >>>> table(sapply(dat$Proficiency, is.finite)) > >>> TRUE > >>> 14082 > >>> > >>> Sorry if this is a dumb question, but can you help me understand > >>> what's going on? > >>> > >>> Many thanks. > >>> > >>> Cecile > [[alternative HTML version deleted]] ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.