После вычисления распределения для максимального абсолютного коэффициента корреляции между v (длиной n = 20) и КАЖДОЙ из других m переменных (где m принимает значения от 1 до 100), как я могу построить m на x-ось и вероятность получения максимальной абсолютной корреляции выше 0,4 по оси у?(предпочтительнее использовать ggplot)
с использованием исторических данных, я получаю частоту корреляций, а не то, что я ищу, и получаю следующее сообщение об ошибке:
Ошибка в xy.coords (x, y, xlabel, ylabel, log): длины 'x' и 'y' различаются
это то, что я пробовал: мне известно только об использовании plot или hist для построения графиков.
set.seed(42) # default set by the question
v <- rnorm (20) # length of variable v
nrun <- 100 # to run m from 1 to 100 to generate variables
corr_vm <- numeric(nrun)
# I want to calculate the max abs correlation coefficient between v
# and each of the other m variables (where m takes on values from 1
# 100)
for ( i in 1: nrun) {
m_norm <- rnorm(length(v))
corr_vm[i] <- cor(v, m_norm)
}
# Find the frequency of the correlation coefficients > 0.4
sum(abs(max(corr_vm)) > 0.4) / nrun
# Plot m on the x axis and frequency of correlation
# coefficients > 0.4 on the y axis
plot(1:100, abs(max(corr_vm)) > 0.4)
Предполагается, что я получу график m на оси x и вероятность получения максимальной абсолютной корреляции выше 0,4 на оси y.Я предполагаю, что m относится к числам 1,2,3,4 ... 100 и вероятность должна быть на оси y (правильно ли я интерпретирую вопрос?)