У меня есть фрейм данных с 9 столбцами, я хочу сформировать таблицу сопряженности и критерий хи-квадрат для каждой комбинации двух столбцов и собрать значения и интерактивные столбцы.
Пример: если у меня есть 4 столбца X, Y, Z и P. Я хочу создать таблицу сопряженности между X-Y, X-Z, Y-P, P-X, P-Y.
В моих данных у меня 9 столбцов.
Это мой код:
pvalues=numeric(0)
interactions=character(0)
for(i in 1:9)
{
{
for(z in i+1:9)
contingencytable=table(data1[,i],data1[,z])
print(contingencytable)
output=chisq.test(contingencytable)
pvalues=c(pvalues,output$p.value)
interactions=c(interactions,paste(colnames(data1[i]),colnames(data1[z])))
}
}
Однако в выходных данных у меня есть только 8 значений и 8 таблиц сопряженности. У меня должно быть намного больше восьми, потому что это комбинация.
Что я должен исправить в своем коде?
Заранее спасибо