Я хочу проверить, есть ли различие между диетами и счетом 1, счетом 2 и счетом 3 по ненормально распределенным данным. Это дизайн с повторными измерениями, поэтому я думаю, что буду использовать тест Фридмана. Интересно, возможно ли включить в модель ковариаты?
set.seed(8)
id <- rep(1:3,each=2)
diet <- rep(c("A","B"),3)
period <- rep(c(1,2),3)
score1 <- sample(1:100,6)
score2 <- sample(1:100,6)
score3 <- sample(1:100,6)
df <- data.frame(id,diet,period,score1, score2, score3)
df
library(tidyverse)
library(ggpubr)
library(rstatix)
colnames(df)
df <- df[,c(1,2,4)]
df %>% friedman_test(score1 ~ diet |id)
Хотелось бы узнать, есть ли разница между периодами, возможно ли это? Если это так, как это может быть включено в модель? Это не работает:
df <- df[,c(1:4)]
df %>% friedman_test(score1 ~ period + diet |id)
Error in friedman.test.formula(formula, data = data, ...) :
incorrect specification for 'formula'