Eu tenho um conjunto de dados de referências e subamostras em cada referência. Eu corro esses benchmarks e suas subamostras emmáquinas sujeitas. O 'indivíduo' estudado pelas subamostras é o mesmo para cada máquina de assunto e os benchmarks são os mesmos para cada máquina de assunto.
Como realizo uma ANOVA em R nesta situação?
Quero principalmente calcular a média total e os intervalos de confiança. Não me importo com os meios de subamostra, mas quero reconhecer a replicação na confiança e nos meios finais. Talvez eu me importe com os meios de referência. Não consigo descobrir como configurar esta anova em R. Quero poder replicar os meios pelo cálculo manual.
Eu tentei glm
, anova
, aov
, e lme
, mas eu estou totalmente confuso. Eu acho que os resultados da ANOVA devem ser equivalentes para duas máquinas sujeitas à média aninhada da máquina / ponto de referência / ponto de verificação, mas os meios não saem da mesma forma quando eu os experimento.
Editar:
Estou começando a ter uma pista de http://zoonek2.free.fr/UNIX/48_R/13.html
fonte
Respostas:
A principal diferença entre o desenho de parcelas subdivididas e outros desenhos, como desenho completamente aleatório e variações de desenhos de blocos, é a estrutura de nidificação dos sujeitos, ou seja, quando as observações são obtidas do mesmo sujeito (unidade experimental) mais de uma vez. Isso leva a uma estrutura de correlação dentro de um sujeito no design de plotagem dividida, diferente da estrutura de correlação em um bloco.
Vamos tirar uma foto de exemplo do conjunto de dados a partir de um design simples de plotagem dividida (abaixo). Este é um estudo da composição da dieta em saúde, quatro dietas foram distribuídas aleatoriamente em 12 indivíduos, todos com status de saúde semelhante. A pressão arterial basal foi estabelecida e uma medida de saúde foi a alteração da pressão arterial após duas semanas. A pressão arterial foi medida de manhã e à noite. (O exemplo é copiado do exemplo do livro Statistical Design da Casella, 5.1)
Algumas coisas importantes a serem observadas:
cada um dosx s aqui são assuntos diferentes. Isso ilustra o conceito de aninhamento. Ou seja, os sujeitos 1, 2, 3 estão aninhados na dieta 1. - As parcelas completas, as unidades experimentais no nível de plotagem (dieta) (os sujeitos) atuam como blocos para o tratamento de plotagem dividida (manhã e noite)
O modelo para esse design de plotagem dividida é:
R
aov
formulário é trivial:fonte
Diet
eTime*Diet
naaov
ligação? DeveTime:Diet
corresponder à sua fórmula matemática.