ANOVA-ошибка смешанного дизайна: объект типа «замыкание» не может быть подмножеством

Question

Я использую смешанный дизайн ANOVA 2X3X3.Я пытался заставить функцию Куэлса ученика-новичка работать в R, но постоянно получаю сообщение об ошибке: объект типа «замыкание» не может быть подмножеством.Любая помощь будет принята с благодарностью!

library(ez)
library(ggplot2)
library(nlme)
library(pastecs)
library(reshape)
library(WRS)
library(GAD)
library(multcomp)
library(psych)
library(lmerTest)
CPPData <- read.table(file = "CPPdatar.csv", header=TRUE, sep = ",")
str(CPPData)
CPPData$Test<-as.factor(CPPData$Test)
CPPData$Exposure<-as.factor(CPPData$Exposure)
CPPData$Dose<-as.factor(CPPData$Dose)
CPPData$Subject<-as.factor(CPPData$Subject)
levels(CPPData$Test)<-list("Habit"=1, "Test of Conditioning"=2)
levels(CPPData$Exposure)<-list("0% HFCS CONT"=1, "50% HFCS CONT"=2, "50% HFCS INT"=3)
levels(CPPData$Dose)<-list("0 OXY"=1, "0.16 OXY"=2, "2.5 OXY"=3)
str(CPPData)
options(contrasts=c("contr.helmert", "contr.poly"))
aov1<-aov(Time~Exposure*Dose*Test+ Error(Subject/(Test)), data=CPPData)
summary(aov1)

aov2 = lme(Time~Exposure + Dose + Test, 
         random = ~1|Subject, 
         data = CPPData, 
         method = "ML")
summary(aov2)

#SNK Test
Test2<-snk.test(lm(Time~Exposure*Test*Dose, data = CPPData))
Test1<-snk.test(lm, term = Exposure*Test*Dose, among = Exposure, within = Test*Dose)
snk.test(aov2, term = "Dose:Exposure:Test", among = "Test", within = "Dose")

Ошибки:

Тест SNK

T est2 <-snk.test (лм (Время ~ Экспозиция <em>Тест *)1011 * Доза, данные = CPPData)) Ошибка в оценках (объект): Дизайн несбалансирован!Эта функция может работать только с сбалансированным дизайном. Test1 <-snk.test (лм, термин = Экспозиция <em>Тест Доза, среди = Экспозиция, в пределах = Тест * Доза) Ошибка: объект типа 'замыкание' не может быть подмножеством > snk.test (aov2, term = "Доза: Экспозиция: Тест", среди = "Тест", в пределах = "Доза) Ошибка в объекте $ model [, 2: (длина (объект $ x) + 1)]: неверное число измерений

Dput

dput (CPPData) структура (список (время = c (476,98, 436,94, 451,79, 514,68, 548,38, 457,96,489,99, 536,7, 517,02, 566,9, 487,15, 553,89, 527,86, 580,41, 310,14, 508,84, 364,7, 456,79, 616,12, 598,43, 447,45, 570,9, 520,35, 812,48, 756,92, 667,5, 603,77, 547,55568,74, 580,91, 526,19, 582,92, 568,23, 576,74, 378,04, 549,38, 548,55, 492,99, 388,72, 581,75, 538,71, 511,51, 895,89, 851,68, 685,85, 741,24, 738,4, 662,57, 638,65, 638,85, 638,85, 638,85, 638,25, 6,25, 6,28, 6,38, 6,35, 6,338, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 6, 5, 6, 6, 5, 6, 6, 5, 5, 6, 5, 5, 9, 5, 5, 5, 9, 5, 5, 9, 5, 5, 8, 8, 5, 8, 5, 6, 5, 5, 6, 5, 5, 5, 5, 5, 5, 5, 6, 5, 5, 6, 5, 5, 5, 6, 5, 5, 5, 5, 5, 6, 5, 5, 6, 5, 6, таков, 5-ых годах видеоблоков видеах измерениях входят:819,65, 475,14, 599,93, 668,34, 486,99, 426,26, 570,74, 482,48, 460,13, 578,08, 541,37, 465,63, 297,13, 543,04, 560,89, 536,04, 419,25, 555,89, 587,75, 530,2, 619,95, 778,61, 602,27, 403,74, 792,63, 815,98, 599,77, 784,28, 543,04, 748,08, 610,94, 700,87, 710,38, 686,02, 735,57, 805,47, 525,02, 524,86, 482,65, 278,28, 547,21, 559,23, 450,95,579,91, 174,17, 312,64, 683,02, 790,46, 422,92, 664, 435,44, 478,14, 427,43, 308,31, 342,84, 517,02, 541,54, 565,73, 474,81, 522,02, 503,34, 321,99, 498,16, 554,55, 527,69, 776,44, 966,47, 867,2,846,68, 554,22, 750,75, 435,44, 769,94, 561,39, 544,71, 404,24, 464,8, 482,48, 550,88, 459,79, 552,05, 389,56, 361,86, 814,15, 907,41, 755,92, 729,9, 945,11, 459,39, 529,88, 6,88,81, 5,88, 6, 8,88, 6, 8, 8, 8, 6, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 6, 8, 8, 7, 8, 8, 8], 5, 8, 8, 8, 6, 8, 6, 8, 6, 8, 6, 8, 6, 8, 8, 6, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 8, 6, 8, 6, 4,2 и 6 раза.584,42, 566,4, 420,75, 465,46, 533,87, 579,58, 492,16, 521,86, 565,07, 594,76, 280,95, 563,23, 537,37, 409,07, 439,44, 547,88, 543,88, 496,66, 474,14, 433,48, 54,77, 541, 54,77, 541, 54,77, 54,75, 54,77, 541, 54,77, 541, 54,77, 54,77, 541,47, 54,77, 54,77, 54,77, 54,77335,5, 405,74, 319,49, 482,15, 501, 409,07, 610,28, 667,17, 560,56, 474,14, 467,3, 561,56, 718,88, 584,25, 689,52, 507,51, 711,04, 520,19, 357,52, 373,04, 504, 360,86, 470,9582,41, 554,05, 481,15, 776,78, 658,66, 721,55, 433,93, 532,03, 625,79, 651,48, 783,78, 729,06, 828,5), Test = структура (c (1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1л, 1л, 1л, 1л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 1л, 1л, 1л, 1л, 1л, 1л, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1л, 1л, 1л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 1л, 1л, 1л, 1л, 1л, 1л, 1л, 1л, 1л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 1л, 1л, 1л, 1л, 1л, 1л, 1л, 1л, 1л, 1л, 1л, 1л, 2л, 2л, 2л, 2л, 2л, 2л, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c («Привычка», «Испытание кондиционирования»), классs = «фактор»), Экспозиция = структура (c (2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л,2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л,2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л,3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л,3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 1л, 1л, 1л, 1л, 1л, 1л,1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c («0% HFCS CONT», «50% HFCS CONT», «50% HFCS INT»), класс = «фактор»), Доза = структура (c (1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 3л, 3л, 3л, 3л, 3л, 3л, 3л,3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л,1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L,2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л, 3л,3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L,1л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л, 2л,3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L), .Label = c ("0 OXY"), "0.16 OXY", "2.5 OXY"), класс = "фактор"), субъект = структура (c (1L, 2L, 3L, 4L, 5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 13L, 14л, 15л, 1л, 2л, 3л, 4л, 5л, 6л, 7л, 8л, 9л, 10л, 11л, 12л, 13л, 14л, 15л, 17л, 18л, 19л, 20л, 21л, 22л, 23л, 24л, 25л, 26л, 27л, 28л, 29л, 30л, 16л, 17л, 18л, 18л, 19л, 20л, 21л, 22л, 23л, 24л, 25л, 27л, 28л, 29л, 30л, 31л, 32л, 33л, 34 л, 35 л, 36 л, 37 л, 38 л, 39 л, 40 л, 41L, 42л, 43л, 44л, 45л, 46л, 31л, 32л, 33л, 34л, 35л, 36л, 37л, 38л, 39л, 40л, 41л, 44л, 45л, 46л, 47л, 48л, 49л,50л, 51л, 52л, 53л, 54л, 55л, 56л, 47л, 48л, 49л, 50л, 51л, 52л, 53л, 54л, 55л, 56л, 57л, 58л, 59л, 60л, 61л, 62л, 63л, 64л,65л, 57л, 58л, 59л, 60л, 61л, 62л, 63л, 64л, 65л, 65л, 66л, 67л, 68л, 69л, 70л, 71л, 72л, 73л, 74л, 66л, 67л, 68л, 69л, 70л, 71л,72л, 73л, 74л, 75л, 76л, 77л, 78л, 79л, 80л, 81л, 82л, 83л, 84л, 75л, 76л, 77л, 79л, 80л, 81л, 82л, 83л, 84л, 85л, 86л,87л, 88л, 89л, 90л, 91л, 92л, 93л, 94л, 95л, 96л, 85л, 86л, 87л, 88л, 89л, 90л, 93л, 94л, 95л, 96л, 97л, 98л, 99л,100л, 101л, 102л, 103л, 104л, 105л, 106л, 97л, 98л, 99л, 100л, 101л, 102л, 103л, 104л, 105л, 106л), .Label = c ("1", "2", "3"," 4 "," 5 "," 6 "," 7 "," 8 "," 9 "," 10 "," 11 "," 12 "," 13 "," 14 "," 15 ",«16», «17», «18», «19», «20», «21», «22», «23», «24», «25», «26», «27», «28"," 29 "," 30 "," 31 "," 32 "," 33 "," 34 "," 35 "," 36 "," 37 "," 38 "," 39 "," 40 ",«41», «42», «43», «44», «45», «46»,«47», «48», «49», «50», «51», «52», «53», «54», «55», «56», «57», «58», «59"," 60 "," 61 "," 62 "," 63 "," 64 "," 65 "," 66 "," 67 "," 68 "," 69 "," 70 "," 71 ",«72», «73», «74», «75», «76», «77», «78», «79», «80», «81», «82», «83», «84"," 85 "," 86 "," 87 "," 88 "," 89 "," 90 "," 91 "," 92 "," 93 "," 94 "," 95 "," 96 ",«97», «98», «99», «100», «101», «102», «103», «104», «105», «106»), класс = «фактор»)), строка.names = c (NA, -212L), класс = "data.frame")

Answer 1

Причину ошибки, о которой вы спрашивали, можно увидеть здесь:

Test1 <- snk.test (lm, term = ExposureTestDose, .... </p>

Буквы «lm» - это имя функции, и вы дали ее функции, которая ожидала объект данных, который может быть подмножеством. Первый вызов дал функции snk.test объект данных, созданный lm,но он был остановлен кодом, который дал другую причину.

Что касается других ошибок, то, похоже, вам нужен какой-то статистический совет. snk.test, похоже, не предназначен для того типа анализа, который выВы можете обратиться за помощью по вопросам дизайна на сайте CrossValidated.com. (Это еще один форум StackExchange, который занимается статистическими вопросами.)