Hi, thanks so much: pairs(sapply(MAR.omitindep[,c("poldis" , "polres" , "ecdis" , "culres" , "gcc1" )],jitter,amount=1), lower.panel=panel.smooth, upper.panel=panel.cor)
Mehmet. On Mon, 25 Feb 2008 09:05:28 +0100 Petr PIKAL <[EMAIL PROTECTED]> wrote: > Hi > > Add > > sapply(any.data.frame, jitter) > > into your function. Either sapply(na.omit(.....), jitter) or > pairs(sapply(..., jitter), ...) > > Regards > > Petr > [EMAIL PROTECTED] > > [EMAIL PROTECTED] napsal dne 25.02.2008 00:26:27: > > > Hi, > > > > I am just starting to use R for a graduate course, and I like how > > the correlation matrix at > > http://addictedtor.free.fr/graphiques/RGraphGallery.php?graph=137 > > > > I did something similar by copying from the examples(pairs), but it > > seems that I need to jitter the bottom panel... and I have no idea > > how to do that, and I mean no idea at all. I'd appreciate any > > help... > > > > Here are the graphs: > > http://socy602.pbwiki.com/f/dep_correlationmatrix.jpeg > > http://socy602.pbwiki.com/f/indep_correlationmatrix.jpeg > > > > And this is the code I used to produce the graphs (not that I > > understand it, but): > > > > ###################### > > # Primary Component Analyses > > ###################### > > > > ###################### > > # Dependent variables > > ###################### > > > > # List of these variables: > > # langfamr , ethdifxx , catness , gc7 , gc8r , gc12 , culdifxx , > poldifxx , ecdifxx > > ###################### > > # Look at correlation matrix for these > > cor(MAR[,c("langfamr" , "ethdifxx" , "catness" , "gc7" , "gc8r" , > > "gc12" > , > > "culdifxx" , "poldifxx" , "ecdifxx" )], use="complete.obs") > > # Big time correlation matrix > > 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) > > text(0.5, 0.5, txt, cex = cex * r) > > } > > MAR.omitdep <- na.omit(MAR[,c("langfamr" , "ethdifxx" , "catness" , > "gc7" , > > "gc8r" , "gc12" , "culdifxx" , "poldifxx" , "ecdifxx" )]) > > pairs(MAR.omitdep[,c("langfamr" , "ethdifxx" , "catness" , "gc7" , > "gc8r" , > > "gc12" , "culdifxx" , "poldifxx" , "ecdifxx" )], > lower.panel=panel.smooth, > > upper.panel=panel.cor) > > # Save it > > dev.copy(jpeg,filename="dep_correlationmatrix.jpeg",height=600, > width=800,bg="white") > > dev.off() > > > > ###################### > > # Independent variables > > ###################### > > > > # List of variables > > # poldis , polres , ecdis , culres , gcc1 > > ###################### > > # Look at the correlation matrix > > cor(MAR[,c("poldis" , "polres" , "ecdis" , "culres" , "gcc1" )], > use="complete.obs") > > # Big time correlation matrix > > 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) > > text(0.5, 0.5, txt, cex = cex * r) > > } > > MAR.omitindep <- na.omit(MAR[,c("poldis" , "polres" , "ecdis" , > > "culres" > , "gcc1" )]) > > pairs(MAR.omitindep[,c("poldis" , "polres" , "ecdis" , "culres" , > > "gcc1" > )], > > lower.panel=panel.smooth, upper.panel=panel.cor) > > # Save it > > dev.copy(jpeg,filename="indep_correlationmatrix.jpeg",height=600, > width=800,bg="white") > > dev.off() > > > > Thanks a lot in advance, > > Sincerely, > > Mehmet. > > > > ______________________________________________ > > R-help@r-project.org mailing list > > https://stat.ethz.ch/mailman/listinfo/r-help > > PLEASE do read the posting guide > http://www.R-project.org/posting-guide.html > > and provide commented, minimal, self-contained, reproducible code. > ______________________________________________ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.