Не думаю, что мне нужно точно объяснять, что делает код. Дело в том, что, выполняя chisq.test вне цикла, я получаю такой результат (ожидаемый):
Chi-squared test for given probabilities
data: observed
X-squared = 185912, df = 5, p-value < 2.2e-16
но когда я пытаюсь выполнить тест в цикле, ожидаемый результат не появляется
total <- dim(crs$dataset_init)[1]
expected.fr <- cl.popul / total
for (i in 1:dim(cl.vs.Onerall)[1] ) {
if (cl.vs.Onerall[i,1] > 0) {
observed <- cl.vs.Onerall[i,2:(clust_no + 1)]
print(rownames(cl.vs.Onerall)[i])
chisq.test(observed, p=expected.fr)
print("------------------------------")
}
}
Любые идеи будут с благодарностью!