[Rd] Lemon drops

Achim Zeileis Achim.Zeileis at wu-wien.ac.at
Tue Mar 15 20:05:40 CET 2005


On 15 Mar 2005 18:26:46 +0100 Peter Dalgaard wrote:

> Marc Schwartz <MSchwartz at 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.

Maybe:
  coef[,1,drop=FALSE][1:nrow(coef)]
?
Z

> 
> -- 
>    O__  ---- Peter Dalgaard             Blegdamsvej 3  
>   c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
>  (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45)
>  35327918
> ~~~~~~~~~~ - (p.dalgaard at biostat.ku.dk)             FAX: (+45)
> 35327907
> 
> ______________________________________________
> R-devel at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-devel
>



More information about the R-devel mailing list