From p.dalgaard at biostat.ku.dk  Tue Mar 15 18:26:46 2005
From: p.dalgaard at biostat.ku.dk (Peter Dalgaard)
Date: Tue Mar 15 18:33:43 2005
Subject: [Rd] Lemon drops
In-Reply-To: <1110907737.8170.17.camel@horizons.localdomain>
References: <x2r7igon5j.fsf@biostat.ku.dk>
	<1110907737.8170.17.camel@horizons.localdomain>
Message-ID: <x2mzt4om5l.fsf@biostat.ku.dk>

Marc Schwartz <MSchwartz@MedAnalytics.com> writes:

> On Tue, 2005-03-15 at 18:05 +0100, Peter Dalgaard wrote:
> > I bumped into the following situation:
> > 
> > Browse[1]> coef
> >             deg0NA deg4NA deg8NA deg0NP deg4NP deg8NP
> > (Intercept)    462    510    528    492    660    762
> > Browse[1]> coef[,1]
> > [1] 462
> > Browse[1]> coef[,1,drop=F]
> >             deg0NA
> > (Intercept)    462
> > 
> > where I really wanted neither, but
> > 
> > (Intercept)
> >         462
> > 
> > Anyone happen to know a neat way out of the conundrum?
> > 
> > I can think of
> > 
> > rowSums(coef[,1,drop=F])
> > 
> > or of course
> > 
> > val <- coef[,1]
> > names(val) <- rownames(x))
> > 
> > but the first one is sneaky and the second gets a bit tedious...
> 
> Peter,
> 
> How about something like this:
> 
> > x
>             deg0NA deg4NA deg8NA deg0NP deg4NP deg8NP
> (Intercept)    462    510    528    492    660    762
> IV1              1      2      3      4      5      6
> 
> 
> > x[1]
> (Intercept)
>         462
> 
> > str(x[1])
>  Named num 462
>  - attr(*, "names")= chr "(Intercept)"
> 
> 
> Does that work or am I missing something?

You're missing the context. What I really need is something that will
extract a column of a matrix as a vector in the usual way, but will
not get confused if there is only one row.


-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk)             FAX: (+45) 35327907


[Rd] Lemon drops
Peter Dalgaard