Создание функции для проверки независимости двух переменных от совместной таблицы распределения вероятностей

#r #function #logic #probability-distribution

Вопрос:

Предположим, что у нас есть приведенная ниже совместная таблица распределения вероятностей двух переменных X и Y .

 # you can use this table as an example
distr_table <- data.frame(
  X = rep(c(0, 1), each = 2),
  Y = rep(c(1, 2), times = 2),
  pr = c(0.3, 0.25, 0.15, 0.3)
)
col<-length(unique(distr_table$Y))
row<-length(unique(distr_table$X))
p <- matrix(c(distr_table$pr),ncol=col,byrow=TRUE)
check_independence <- function(distr_table) {
  # write your solution here
    for (i in 1:col){
    for (j in 1:row){
      if(p[i,j]==sum(p[i,])*sum(p[,j])){
      TRUE
      }else{
      FALSE
      }
    }
    }
return(are_independent)
}  

) 

 

Я хочу создать функцию, которая будет возвращать переменную (с именем check_independence), независимо от того, являются ли x и y независимыми или нет(ИСТИНА/ЛОЖЬ). Как я могу это сделать? Код, который я предоставил, где я пытался получить check_independence(TRUE/FALSE), который, кажется, не работает.