Comment trouvez-vous la taille des échantillons utilisés dans les calculs de r?

Je suis en cours d'exécution des corrélations entre les variables, dont certaines ont des données manquantes, de sorte que la taille de l'échantillon pour chaque corrélation sont probablement différentes. J'ai essayé d'imprimer et d'un résumé, mais aucun de ces montre-moi comment taille de mon n est pour chaque corrélation. C'est assez simple problème que je ne trouve pas la réponse à n'importe où.

  • êtes-vous à l'aide de cor(...,use="pairwise.complete.obs")) ... ?
  • suite à @BenBolker commentaire que vous aimeriez voir tinyurl.com/afcglqv
InformationsquelleAutor robert wilson | 2013-01-01
  1. 3

    comme ça..? 

    x <- c(1:100,NA)
length(x)
length(x[!is.na(x)])

    vous pouvez également obtenir les degrés de liberté comme ça...

    y <- c(1:100,NA)
x <- c(1:100,NA)

cor.test(x,y)$parameter

    Mais je pense qu'il serait mieux si vous montrer le code pour la façon dont vous êtes l'estimation de la corrélation exacte de l'aide.

    • ..et ne serait pas la taille de l'échantillon juste être (degrés de liberté)+2 pour cet? 🙂 donc la réponse est.. cor.test( x , y )$parameter + 2
    • oui , l'estimation de deux paramètres de d.f. dans mon cas, c'est N-2, j'ai pensé à cette solution était implicite
    • bien que pour method="spearman" ou "kendall" la cor.test(x,y)$parameter sera de retour NULL
    InformationsquelleAutor user1317221_G
  2. 0

    Voici un exemple de la façon de trouver les paires d'échantillons de tailles entre les colonnes d'une matrice. Si vous voulez l'appliquer à (certains) des colonnes numériques d'une trame de données, de les combiner en conséquence, de contraindre l'objet résultant de la matrice et d'appliquer la fonction. 

    # Example matrix:
xx <- rnorm(3000)
# Generate some NAs
vv <- sample(3000, 200)
xx[vv] <- NA
# reshape to a matrix
dd <- matrix(xx, ncol = 3)
# find the number of NAs per column
apply(dd, 2, function(x) sum(is.na(x)))
# tack on some column names
colnames(dd) <- paste0("x", seq(3))

# Function to find the number of pairwise complete observations 
# among all pairs of columns in a matrix. It returns a data frame
# whose first two columns comprise all column pairs

pairwiseN <- function(mat)
{
    u <- if(is.null(colnames(mat))) paste0("x", seq_len(ncol(mat))) else colnames(mat)
    h <- expand.grid(x = u, y = u)

    f <- function(x, y)
           sum(apply(mat[, c(x, y)], 1, function(z) !any(is.na(z))))
    h$n <- mapply(f, h[, 1], h[, 2])
    h
}

# Call it
pairwiseN(dd)

    La fonction peut être facilement amélioré; par exemple, vous pouvez définir h <- expand.grid(x = u[-1], y = u[-length(u)]) à réduire le nombre de calculs, vous pouvez retourner un n x n de la matrice au lieu de trois colonnes de données d'image, etc.

    InformationsquelleAutor Dennis
  3. -1

    Si vos variables sont des vecteurs nommé a et b, serait quelque chose comme sum(is.na(a) | is.na(b)) vous aider?

    InformationsquelleAutor Theodore Lytras