# Re: [R] elegant solution to transform vector into percentages?

Date: Wed 27 Jul 2005 - 08:21:15 EST

Why not write a function ? Here is one.

mytable <- function(x, br){

n <- length(br)

tb <- table( cut( x, breaks=c(-Inf, br, Inf) ) )    tb <- 100 * tb / sum(tb)

```   tb.n      <- paste( c("", br), c(br, ""), sep=" < x <= " )
tb.n   <- paste("x <= ", br, sep="")
tb.n[n+1] <- paste("x > ",  br[n], sep="")

```

names(tb) <- tb.n
return(tb)
}

mytable( store, br=c(0.8, 1.2) )

```       x <= 0.8 0.8 < x <= 1.2        x > 1.2
25             25             50

```

You might want to do a bit more testing especially at the break points. However I do not like the output of the above because the names of the table overlaps. Here is another function that might have nicer output.

mytable2 <- function(x, br){
tb <- as.matrix( mytable( x=x, br=br ) )    colnames(tb) <- "Percentage"
return(tb)
}

mytable2( store, br=c(0.8, 1.2) )

```                Percentage
x <= 0.8               25
0.8 < x <= 1.2         25
x > 1.2                50

```

Hope this helps.

On Tue, 2005-07-26 at 15:48 -0400, jhainm@fas.harvard.edu wrote:
> Hi,
>
> I am looking for an elegant way to transform a vector into percentages of values
> that meet certain criteria.
>
> store<-c(1,1.4,3,1.1,0.3,0.6,4,5)
>
> # now I want to get the precentages of values
> # that fall into the categories <=M , >M & <=N , >N
> # let
> M <-.8
> N <- 1.2
> # In my real example I have many more of these cutoff-points
>
> # What I did is:
>
> out <- matrix(NA,1,3)
>
> out[1,1] <- ( (sum(store<=M )) /length(store) )*100
> out[1,2] <- ( (sum(store> M & store<= N )) /length(store) )*100
> out[1,3] <- ( (sum(store> N )) /length(store) )*100
>
> colnames(out)<-c("percent<=M","percent>M & <=N","percent>N")
> out
>
> But this gets very tedious if I have many cutoff-points. Does anybody know a
> more elegant way to do this task?
>
> Thanks so much.
>
> Cheers,
> Jens
>
> ______________________________________________
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help