Re: [R] Counting observations split by a factor when there are NAs in the data

From: Peter Dalgaard <>
Date: Tue 11 Jul 2006 - 05:47:06 EST

"Jenifer Larson-Hall" <> writes:

> I am a very novice R user, a social scientist (linguist) who is trying
> to learn to use R after being very familiar with SPSS. Please be kind!
> My concern:
> I cannot figure out a way to get an accurate count of observations of
> one column of data split by a factor when there are NAs in the data.
> I know how to use commands like tapply and summaryBy to obtain other
> summary statistics I am interested in, such as the following:
> tapply(RLWTEST, list(STATUS), mean, na.rm=T)
> summaryBy(RLWTEST~STATUS, data=lh.forgotten, FUN=c(mean, sd, min, max),
> na.rm=T)
> However, with tapply I know I cannot use length to get a count where
> there are NAs. summaryBy appears to work the same way. I do know how to
> get a count of the entire column using sum:
> sum(!$RLWTEST))
> However, this does not give me a count split up by my factor (STATUS). I
> have looked through Daalgard (2002) and Verzani (2005), and have


> searched the help files, but with no luck.

How about

  tapply(!, STATUS, sum)

or maybe just


> Thank you in advance for your help. I love R and am interested in making
> it more accessible to social scientist types like me. I know it can do
> everything SPSS can and more, but sometimes the very simplest things
> seem to be a lot harder in R.

   O__  ---- Peter Dalgaard             ุster Farimagsgade 5, Entr.B
  c/ /'_ --- Dept. of Biostatistics     PO Box 2099, 1014 Cph. K
 (*) \(*) -- University of Copenhagen   Denmark          Ph:  (+45) 35327918
~~~~~~~~~~ - (                  FAX: (+45) 35327907

______________________________________________ mailing list
PLEASE do read the posting guide!
Received on Tue Jul 11 05:52:46 2006

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Tue 11 Jul 2006 - 06:16:30 EST.

Mailing list information is available at Please read the posting guide before posting to the list.