Греческая бета в легенде помимо коэффициента OLS - PullRequest
Новая
       19

Греческая бета в легенде помимо коэффициента OLS

0 голосов
/ 01 октября 2018

Я хочу отобразить три линейных подгонки на графике.В легенде я хочу отобразить коэффициент корреляции помимо метки соответствующей группы.В идеале, я бы имел греческую бета и = вместо "rank corr:" перед фактическим коэффициентом в метке.Кто-нибудь знает, как этого добиться?Большое спасибо. 

n <- 100
x <- rnorm(n)
y <- rnorm(n, mean = x)
grp <- sample(LETTERS[1:3], replace = TRUE, size = n)

# Do linear regressions for each group
lin_regs <- lapply(unique(grp), function(g) lm(y ~ x, subset = g==grp))

b1 <- format(summary(lin_regs[[1]])$coefficients[2, 1], digits=3)
b2 <- format(summary(lin_regs[[2]])$coefficients[2, 1], digits=3)
b3 <- format(summary(lin_regs[[3]])$coefficients[2, 1], digits=3)

plot1 <- sprintf("Basic, rank corr: %s", b1)
plot2 <- sprintf("Intermediate, rank corr: %s", b2)
plot3 <- sprintf("High, rank corr: %s", b3)

# Plot lines
plot(x, y, type = "n")  # Initialize plot
for (i in seq_along(lin_regs)) {
   abline(lin_regs[[i]], col = i) # Plot each linear reg. with abline()
}
legend('bottom', legend = c(plot1,plot2,plot3),
       col=c('black', 'red', 'green'), lty=1)
...