计算分布(的长度n = 20)v和每个其它m个变量的(其中m取数值上从1到100)之间的最大绝对相关系数,如何在x做我情节米后轴和在y轴上获得最大绝对相关性高于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)
我期望在x轴上绘制出m的图,并在y轴上获得大于0.4的最大绝对相关性的概率。我假设米指数1,2,3,4 ... 100和概率应该是在y轴上(我是正确解释的问题?)