Re: [R] number of pairwise present data in matrix with missings
Andreas Wolf andreas.wolf at uni-jena.de writes:
:
: is there a smart way of determining the number of pairwise present data
: in a data matrix with missings (maybe as a by-product of some
: statistical function?)
:
: so far, i used several loops like:
:
: for (column1 in 1:99) {
: for (column2 in 2:100) {
: for (row in 1:500) {
: if (!is.na(matrix[row,column1]) !is.na(matrix[row,column2])) {
: pairs[col1,col2] - pairs[col1,col2]+1
: }
: }
: }
: }
:
: but this seems neither the most elegant nor an utterly fast solution.
This is just matrix multiplication of the !na(x) matrix:
R x - matrix(1:12,4,3)
R x[c(1,10)] - NA
R x
[,1] [,2] [,3]
[1,] NA59
[2,]26 NA
[3,]37 11
[4,]48 12
R crossprod(!is.na(x))
[,1] [,2] [,3]
[1,]332
[2,]343
[3,]233
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] number of pairwise present data in matrix with missings
Suppose your matrix is called A (`matrix' is not a good name). Then
crossprod(!is.na(A)) is pretty efficient. Test:
A - matrix(1, 6, 3)
A[1,1] - A[3, 1] - A[2,2] - NA
A
[,1] [,2] [,3]
[1,] NA11
[2,]1 NA1
[3,] NA11
[4,]111
[5,]111
[6,]111
crossprod(!is.na(A))
[,1] [,2] [,3]
[1,]434
[2,]355
[3,]456
On Tue, 23 Nov 2004, Andreas Wolf wrote:
is there a smart way of determining the number of pairwise present data
in a data matrix with missings (maybe as a by-product of some
statistical function?)
so far, i used several loops like:
for (column1 in 1:99) {
for (column2 in 2:100) {
for (row in 1:500) {
if (!is.na(matrix[row,column1]) !is.na(matrix[row,column2])) {
pairs[col1,col2] - pairs[col1,col2]+1
}
}
}
}
but this seems neither the most elegant nor an utterly fast solution.
thanks for suggestions.
andreas wolf
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
--
Brian D. Ripley, [EMAIL PROTECTED]
Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/
University of Oxford, Tel: +44 1865 272861 (self)
1 South Parks Road, +44 1865 272866 (PA)
Oxford OX1 3TG, UKFax: +44 1865 272595
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] number of pairwise present data in matrix with missings
Andreas Wolf wrote:
is there a smart way of determining the number of pairwise present data
in a data matrix with missings (maybe as a by-product of some
statistical function?)
so far, i used several loops like:
for (column1 in 1:99) {
for (column2 in 2:100) {
for (row in 1:500) {
if (!is.na(matrix[row,column1]) !is.na(matrix[row,column2])) {
pairs[col1,col2] - pairs[col1,col2]+1
}
}
}
}
but this seems neither the most elegant nor an utterly fast solution.
thanks for suggestions.
andreas wolf
library(Hmisc)
n - naclus(mydataframe)
plot(n) # show pairwise missingness in a dendogram
naplot(n) # show more details in multiple plots
Frank
--
Frank E Harrell Jr Professor and Chair School of Medicine
Department of Biostatistics Vanderbilt University
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] number of pairwise present data in matrix with missings
Hi Andreas,
maybe something like this could do it:
mat - sample(0:3, 20*2, TRUE); dim(mat) - c(20,2)
mat[sample(1:20, 4),] - NA
mat
sum(rowMeans(mat)==mat[,1], na.rm=TRUE)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/16/336899
Fax: +32/16/337015
Web: http://www.med.kuleuven.ac.be/biostat
http://www.student.kuleuven.ac.be/~m0390867/dimitris.htm
- Original Message -
From: Andreas Wolf [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Sent: Tuesday, November 23, 2004 2:42 PM
Subject: [R] number of pairwise present data in matrix with missings
is there a smart way of determining the number of pairwise present
data
in a data matrix with missings (maybe as a by-product of some
statistical function?)
so far, i used several loops like:
for (column1 in 1:99) {
for (column2 in 2:100) {
for (row in 1:500) {
if (!is.na(matrix[row,column1]) !is.na(matrix[row,column2]))
{
pairs[col1,col2] - pairs[col1,col2]+1
}
}
}
}
but this seems neither the most elegant nor an utterly fast
solution.
thanks for suggestions.
andreas wolf
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] number of pairwise present data in matrix with missings
Sorry my first reply was not relevant, I understood a different thing.
Dimitris
- Original Message -
From: Andreas Wolf [EMAIL PROTECTED]
To: [EMAIL PROTECTED]
Sent: Tuesday, November 23, 2004 2:42 PM
Subject: [R] number of pairwise present data in matrix with missings
is there a smart way of determining the number of pairwise present
data
in a data matrix with missings (maybe as a by-product of some
statistical function?)
so far, i used several loops like:
for (column1 in 1:99) {
for (column2 in 2:100) {
for (row in 1:500) {
if (!is.na(matrix[row,column1]) !is.na(matrix[row,column2]))
{
pairs[col1,col2] - pairs[col1,col2]+1
}
}
}
}
but this seems neither the most elegant nor an utterly fast
solution.
thanks for suggestions.
andreas wolf
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
http://www.R-project.org/posting-guide.html
__
[EMAIL PROTECTED] mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
