Предположим, что мои данные выглядят так (это просто данные, иллюстрирующие мой вопрос):
set.seed(1234)
x=runif(12, 22,28);x
y=runif(12,88,120);y
y1=y[1:4];y1
y2=y[5:8];y2
y3=y[9:12];y3
x1=x[1:4];x1
x2=x[5:8];x2
x3=x[9:12];x3
dat= cbind(y1,x1,y2,x2,y3,x3);dat
y1 x1 y2 x2 y3 x3
[1,] 118.72718 27.79948 100.56455 23.59740 95.37221 25.51741
[2,] 90.46189 26.12509 94.38618 27.96536 103.42347 24.44862
[3,] 92.48808 25.67335 108.15713 24.08918 108.85686 26.52390
[4,] 103.06760 25.84996 108.88237 26.37924 103.26130 22.05003
Предположим, что (y1,x1) представляет лечение 1, (y2,x2) представляет лечение 2, а (y3,x3) представляет лечение 3. Моя цель состоит в том, чтобы оценить точки пересечения и наклоны (используя пошаговый план). Затем определите матрицу гипотезы (нулевая гипотеза о том, что все наклоны одинаковы).
Я думал, что могу использовать индикаторы. Например, для x1= 1 при лечении и 0 в противном случае. Тогда x2= 1, если лечение 2, и 0 в противном случае. Вот иначе относится к третьему обращению.
Вот что я придумал:
X< as.matrix(cbind(rep(1,12),c(rep(1,4),rep(0,8)),c(rep(0,4),rep(1,4),rep(0,4))))
Во всяком случае, я не знаю, куда идти отсюда. Мне было интересно, если кто-нибудь даст мне что-нибудь, чтобы начать. Пожалуйста, все обозначения должны быть в матричных формах. Благодарю вас!
Редактировать: я пытаюсь сделать это в матричной нотации, вот что я получил до сих пор, хотя у меня есть некоторые коэффициенты как NA:
set.seed(1234)
x=runif(12, 22,28);x
y=runif(12,88,120);y
y1=y[1:4];y1
y2=y[5:8];y2
y3=y[9:12];y3
x1=x[1:4];x1
x2=x[5:8];x2
x3=x[9:12];x3
dat= cbind(y1,x1,y2,x2,y3,x3);dat
X1<- cbind( dat[,2],rep(0,4),rep(0,4), dat[,2],rep(1,4),rep(1,4));X1
X2 <- cbind(dat[,4],rep(1,4),dat[,4], rep(0,4),rep(1,4),rep(1,4));X2
X3 <- cbind(dat[,6],rep(1,4),dat[,6], rep(0,4),rep(0,4),rep(0,4));X3
X <- rbind(X1,X2,X3)
model <- lm(formula = y ~ X);summary(model)