我有以下代碼來顯示相關矩陣,如何修改此相關矩陣圖?
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)
我想修改類的情節:
有小藍點作爲
pairs(USJudgeRatings[,c(2:3,6,1,7)], main="xxx", pch=18, col="blue", cex=0.8)
包括一個對角線條目的直方圖(如enter link description here所示)
顯示相關和p值作爲
r=0.9; p=0.001;
與值不星。
顯示配對數據散點圖的擬合線。接頭使用什麼方法?以上顯示的代碼將哪一行定義爲擬合?以及如何改變擬合方法?
你問了很多,但你沒有顯示你已經嘗試過。我認爲你有更多的運氣可以在格子包中做到這一點。見'?splom'。 – agstudy 2013-03-07 12:38:07
@agstudy對不起,我對R語言很新。我不知道如何做到這一點。我嘗試了對(USJudgeRatings [,c(2:3,6,1,7)], lower.panel = panel.smooth,upper.panel = panel.cor,pch = 18,col =「blue」),但得到一些錯誤。 – 2013-03-07 12:48:45
顯示配對數據散點圖的擬合線。接頭使用什麼方法?以上顯示的代碼將哪一行定義爲擬合?以及如何改變擬合方法? – 2013-03-07 14:34:34