[R] Matrix formatting
Adaikalavan Ramasamy
gisar at nus.edu.sg
Thu Feb 13 08:43:04 CET 2003
I do something like this as part of a missing value plot where I 'flush'
the rows with most missing spots to the bottom and the columns with the
most missing values to the right. In this manner I can visually see how
much data I will retain (by looking at the area) if I decide to truncate
rows/columns by missing value criterion. You might find this code useful
in rearranging your matrix:
na.mat <- 1*is.na(data)
spots.na.per.row <- rowSums(na.mat)/ncol(data) # calculates the
percentage missing by row
spots.na.per.column <- colSums(na.mat)/nrow(data)
data.re <- data[ order(spots.na.per.row), order(spots.na.per.column) ]
-----Original Message-----
From: Jerome Asselin [mailto:jerome at hivnet.ubc.ca]
Sent: Thursday, February 13, 2003 5:06 AM
To: Tim Sharac; r-help at stat.math.ethz.ch
Subject: Re: [R] Matrix formatting
Would this do what you want?
Cheers,
Jerome
> foo.dat <- cbind(c(NA, 1, 2, 3, 4, 5), c(NA, NA, 0, 10 ,20, 30))
> apply(foo.dat,2,function(x) x[order(as.logical(x))])
[,1] [,2]
[1,] 1 0
[2,] 2 10
[3,] 3 20
[4,] 4 30
[5,] 5 NA
[6,] NA NA
On Wednesday 12 February 2003 12:42, Tim Sharac wrote:
> Content-Length: 918
> Status: R
> X-Status: N
>
> Hi R-users:
>
> I have a data formatting question. I have a data set that looks
> something like this:
>
> foo.dat <- cbind(c(NA, 1, 2, 3, 4, 5), c(NA, NA, 0, 10 ,20, 30))
>
> What I have:
>
> [,1] [,2]
> [1,] NA NA
> [2,] 1 NA
> [3,] 2 0
> [4,] 3 10
> [5,] 4 20
> [6,] 5 30
>
>
> I want to line up the columns by the first value that is not NA. Like
> so:
>
> [,1] [,2]
> [1,] 1 0
> [2,] 2 10
> [3,] 3 20
> [4,] 4 30
> [5,] 5 NA
> [6,] NA NA
>
> Question is: Is there an elegant way to do this without a for loop?
>
> I tried doing this with na.omit and na.exclude without success.
>
> The real data is many hundreds of columns and many thousands of rows.
>
> Thanks in advance, Tim
>
> Sign up for Internet Service under $10 dollars a month, at
> http://isp.BlueLight.com
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> http://www.stat.math.ethz.ch/mailman/listinfo/r-help
______________________________________________
R-help at stat.math.ethz.ch mailing list
http://www.stat.math.ethz.ch/mailman/listinfo/r-help
More information about the R-help
mailing list