Comment modifier cette Matrice de Corrélation de l'intrigue?

J'ai les codes ci-dessous pour afficher une matrice de corrélation,

panel.cor <- function(x, y, digits=2, prefix="", cex.cor)
{
    usr <- par("usr"); on.exit(par(usr))
    par(usr = c(0, 1, 0, 1))
    r <- abs(cor(x, y))
    txt <- format(c(r, 0.123456789), digits=digits)[1]
    txt <- paste(prefix, txt, sep="")
    if(missing(cex.cor)) cex <- 0.8/strwidth(txt)

    test <- cor.test(x,y)
    # borrowed from printCoefmat
    Signif <- symnum(test$p.value, corr = FALSE, na = FALSE,
                  cutpoints = c(0, 0.001, 0.01, 0.05, 0.1, 1),
                  symbols = c("***", "**", "*", ".", " "))

    text(0.5, 0.5, txt, cex = cex * r)
    text(.8, .8, Signif, cex=cex, col=2)
}
pairs(USJudgeRatings[,c(2:3,6,1,7)],
  lower.panel=panel.smooth, upper.panel=panel.cor)

Je veux modifier l'intrigue comme:

  1. Avoir de petits points bleus comme 
    pairs(USJudgeRatings[,c(2:3,6,1,7)],
      main="xxx",
      pch=18,
      col="blue",
      cex=0.8)
  2. Comprennent un histogramme des entrées sur la diagonale (comme on le voit dans entrez description du lien ici)
  3. Affichage de la corrélation et de la p-valeur en 
    r=0.9;
p=0.001;

avec des valeurs qui ne sont pas des étoiles.

Il y a un côté de la ligne affichée pour le nuage de points de la paire de données. Quelle est la méthode utilisée pour le montage? La ligne à laquelle est défini le raccord comme les codes ci-dessus? Et comment faire pour modifier la méthode de l'ajustement?

Vous demandez beaucoup, mais vous ne montrez pas ce que vous avez essayé. Je pense que vous avez plus de chance de le faire dans les treillis de package. voir ?splom.
Désolé, je suis assez nouveau pour R la langue. Je ne suis pas sûr de la façon de le faire. J'ai essayé de paires(USJudgeRatings[,c(2:3,6,1,7)], plus bas.panel=panneau.lisse supérieure.panel=panneau.cor, pch=18, col="blue"), mais a eu quelques erreurs.
Il y a un côté de la ligne affichée pour le nuage de points de la paire de données. Quelle est la méthode utilisée pour le montage? La ligne à laquelle est défini le raccord comme les codes ci-dessus? Et comment faire pour modifier la méthode de l'ajustement?

OriginalL'auteur Samo Jerom | 2013-03-07

  1. 35

    Page d'aide pour la fonction pairs() vous donne un exemple, comment définir des panneaux de la parcelle.

    Pour votre cas particulier:

    Changé panel.cor() fonction pour afficher les lignes de texte - p-valeurs et les coefficients de corrélation.

    panel.cor <- function(x, y, digits=2, cex.cor)
{
  usr <- par("usr"); on.exit(par(usr))
  par(usr = c(0, 1, 0, 1))
  r <- abs(cor(x, y))
  txt <- format(c(r, 0.123456789), digits=digits)[1]
  test <- cor.test(x,y)
  Signif <- ifelse(round(test$p.value,3)<0.001,"p<0.001",paste("p=",round(test$p.value,3)))  
  text(0.5, 0.25, paste("r=",txt))
  text(.5, .75, Signif)
}

    Pour panel.smooth() fonction définie cex=, col= et pch= arguments.

    panel.smooth<-function (x, y, col = "blue", bg = NA, pch = 18, 
                        cex = 0.8, col.smooth = "red", span = 2/3, iter = 3, ...) 
{
  points(x, y, pch = pch, col = col, bg = bg, cex = cex)
  ok <- is.finite(x) & is.finite(y)
  if (any(ok)) 
    lines(stats::lowess(x[ok], y[ok], f = span, iter = iter), 
          col = col.smooth, ...)
}

    Pour ajouter des histogrammes, des panel.hist() les fonctions doivent être définies (prises à partir du fichier d'aide de pairs())

    panel.hist <- function(x, ...)
{
  usr <- par("usr"); on.exit(par(usr))
  par(usr = c(usr[1:2], 0, 1.5) )
  h <- hist(x, plot = FALSE)
  breaks <- h$breaks; nB <- length(breaks)
  y <- h$counts; y <- y/max(y)
  rect(breaks[-nB], 0, breaks[-1], y, col="cyan", ...)
}

    Intrigue finale:

    pairs(USJudgeRatings[,c(2:3,6,1,7)],
          lower.panel=panel.smooth, upper.panel=panel.cor,diag.panel=panel.hist)

    Comment modifier cette Matrice de Corrélation de l'intrigue?

    merci à vous pour cette !
    Personne ne sait comment passer le cex.cor variable dans le pairs() appel? Je suppose que cela est utilisé dans le panel.cor() fonction pour text(), tout en n'étant pas. Mais en ajoutant qu'il génère beaucoup de mises en garde!

    OriginalL'auteur Didzis Elferts

  2. 0

    Modifié Nuage De Points De La Matrice

    1. %% Modifiés en fonction de l'histogramme;

      panel.hist <- function(x, ...)
{
usr <- par("usr"); on.exit(par(usr))
par(usr = c(usr[1:2], 0, 1.5) )
par(cex.axis=2, family="Times New Roman", face="bold", size=12, cex.lab=1, cex.main=1, cex.sub=1)
h <- hist(x, plot = FALSE)
breaks <- h$breaks; nB <- length(breaks)
y <- h$counts; y <- y/max(y)
rect(breaks[-nB], 0, breaks[-1], y, col="cyan", ...)

}

    2. %% Modifié Fonction de Régression avec panel.smooth;

      panel.smooth<-function (x, y, col = "black", bg = NA, pch = 16, 
                cex = 2, col.smooth = "red", span = 2/3, iter = 3, ...) 
{
points(x, y, pch = pch, col = col, bg = bg, cex = cex)
ok <- is.finite(x) & is.finite(y)
if (any(ok)) 
lines(stats::lowess(x[ok], y[ok], f = span, iter = iter), 
      col = col.smooth, ...)
}

    3. %% Modifié Fonction de Corrélation avec panel.cor;

      panel.cor <- function(x, y, digits=2, cex.cor)
{
usr <- par("usr"); on.exit(par(usr))
par(usr = c(0, 1, 0, 1))
r <- abs(cor(x, y))
txt <- format(c(r, 0.123456789), digits=digits)[1]
test <- cor.test(x,y)
Signif <- ifelse(round(test$p.value,3)<0.001,"p < 0.001",paste("p = ",round(test$p.value,3)))  
text(0.5, 0.25, paste("r = ",txt), cex = 2.5, family="Times New Roman", face="bold", size=12)
text(.5, .75, Signif, cex = 2.5, family="Times New Roman", face="bold", size=12)
}

    Pour être en mesure de tracer le nuage de points de la matrice, vous devez également installer "Times New Roman" de la police. Pour le faire, suivez les étapes ci-dessous;

    1. %% Installer toutes les polices dans RStudio. Ceci est important afin d'améliorer la qualité de l'intrigue!

      install.packages("extrafont") # Install fonts
library(extrafont)            # Install library
font_import()                 # Import all fonts  
loadfonts(device="win")       # Register fonts for Windows bitmap output
fonts()                       # Finish the process

    2. %% Enfin, l'intrigue de votre figure avec pairs fonction;

       pairs(qq1, lower.panel=panel.smooth, upper.panel=panel.cor ,diag.panel=panel.hist, cex = 2, cex.labels = 2, cex.main = 2)

    3. %% Vérifier le produit final;
      entrez la description de l'image ici

    OriginalL'auteur cgtydmrc