Dear Rxperts!
Thanks for the suggestions...(esp. Dr. Mathemedician)! :)

I add one more factor level (a level higher to "a") in the example..
a <- data.frame(p=rep(c("A","B"),each=10,len=30),
a=rep(c(1,2,3),each=10),id=seq(30),
b=round(runif(30,10,20)),
c=round(runif(30,40,70)))

tabular(((p=factor(p))*(a=factor(a))+1) ~ (N = 1) + (b + c)*
(mean+sd),data=a)

             b               c
 p a   N  mean  sd    mean  sd
 A 1   10 16.30 2.497 52.30  9.358
   2    0   NaN    NA   NaN     NA
   3   10 15.60 2.716 60.30  8.001
 B 1    0   NaN    NA   NaN     NA
   2   10 15.40 2.366 57.70 10.414
   3    0   NaN    NA   NaN     NA
   All 30 15.77 2.473 56.77  9.601

How do I remove the lines which N=0 or the columns with NaN or NA? I could
not make the subsetting....
I can't generate the above resulting table look like..

             b               c
 p a   N  mean  sd    mean  sd
 A 1   10 16.30 2.497 52.30  9.358
     3   10 15.60 2.716 60.30  8.001
 B  2   10 15.40 2.366 57.70 10.414
   All 30 15.77 2.473 56.77  9.601


Thanks so much,
Santosh


On Fri, Apr 19, 2013 at 3:15 PM, David Winsemius <dwinsemius@comcast.net>wrote:

>
> On Apr 19, 2013, at 2:03 PM, Santosh wrote:
>
> > Rounding was done to replicate the problem I faced in the original data
> set...
>
> I got an error every time I tried to use digits=0. My basic rule is that
> if a parameter choice consistently errors out ... stop using it. I'm
> fundamentally a pragmatist.
>
> > 2 independent identifiers in a group (shown only for 1 of 7 groups) may
> have same value of b or c.
>
> I don't understand what that means  ....  or why it should affect the
> ability of the program to deliver 'a'-level-specific estimates of the
> mean/sd values of 'b' & 'c' in the environment of data.frame-'a'. I showed
> that the problem was in not using factor() on the LHS of the formula when
> the grouping variable is numeric.
>
> --
> David.
>
> >
> >
> > On Fri, Apr 19, 2013 at 1:53 PM, David Winsemius <dwinsemius@comcast.net>
> wrote:
> >
> > On Apr 19, 2013, at 1:29 PM, Santosh wrote:
> >
> > > Dear expeRts!
> > > when I try summarize data using "tabular" of tables package..
> > >
> > > here is an example..
> > >> a <-
> > >
> data.frame(a=rep(10,10),b=round(runif(10,10,20)),c=round(runif(10,40,50)))
> > >> a
> > >    a  b  c
> > > 1  10 14 44
> > > 2  10 17 40
> > > 3  10 13 40
> > > 4  10 18 40
> > > 5  10 11 49
> > > 6  10 16 46
> > > 7  10 15 43
> > > 8  10 19 49
> > > 9  10 19 49
> > > 10 10 13 45
> > >> tabular((a+1)~(n=1)+Format(dig=0)*(b + c) * (mean + sd),data=a)
> > > Error in term2table(rows[[i]], cols[[j]], data, n) :
> > >  Duplicate values: a and b
> > >
> > > Is there a way to let the function  accept the "apparently" duplicate
> > > values?
> >
> > I don't consider myself an expert, just one who reads the manuals and
> makes modest adjustments to the help-page examples.  Perhaps if you
> explained what you thought the correct answer would be ( and can you please
> explain why you are setting digits=0?
> >
> > tabular( 1~Format(digits=2) *(b+c)* (mean + sd),data=a)
> > #-------------------------
> >      b         c
> >      mean sd   mean sd
> >  All 14.9  3.0 44.3  2.9
> > #-----------------
> > # Putting a numeric variable on the LHS appears to fail but wrapping it
> in factor() allows success:
> >
> > >
> a<-data.frame(a=rep(c(1,2),10),b=round(runif(10,10,20)),c=round(runif(10,40,50)))
> > >
> > > tabular( (factor(a) +1)~Format(digits=2) *(b+c)* (mean + sd),data=a)
> > #----------
> >            b         c
> >  factor(a) mean sd   mean sd
> >  1         15.2  3.7 45.0  2.6
> >  2         13.0  3.5 48.2  1.2
> >  All       14.1  3.7 46.6  2.6
> >
> >
> > --
> >
> > David Winsemius
> > Alameda, CA, USA
> >
> >
>
> David Winsemius
> Alameda, CA, USA
>
>

	[[alternative HTML version deleted]]

