Senhores, bom dia!
Alexandre, você falou em 6 áreas, mas seu código está definido 5. Sugiro
verificar acrescentando as linhas que seguem ao seu código:
areas - as.vector(unique(dados$ab)); areas
table(dados$proj.talhao, dados$n.talhao)
Tendo verificado, adicione também as linhas que seguem e
Caros,
Estou tentando resolver um problema que nem sei por onde começar.
Tenho dois data.frame()
data1
name grupo
1Banker, 1986, V32, P1613 1
2 Banker, 1986, V32, P30 1
3 Macmillan, 1986, V60, P41 1
4 Charnes, 1985, V30, P91 1
5
Caros Amigos,
Estou refazendo os comandos apresentado em um artigo de meta-análise.
Acontece que o autor está utilizando o SAS e gostaria de refazer o
comando utilizando o R:
O código utilizado pelo autor no SAS foi:
PROC GLM DATA=Dataregs;
CLASS Study;
MODEL Y = X Study X*Study/SOLUTION;
Caros Amigos,
Estou refazendo os comandos apresentado em um artigo de meta-análise.
Acontece que o autor está utilizando o SAS e gostaria de refazer o
comando utilizando o R:
O código utilizado pelo autor no SAS foi:
PROC GLM DATA=Dataregs;
CLASS Study;
MODEL Y = X Study X*Study/SOLUTION;
Basicamente não existe segredo para reproduzir a mesma análise. O modo mais
direto é declarar o modelo para estimar os parâmetros separados para cada
nível do fator, depois montar as matrizes correspondentes às funções
lineares e submetê-las à glht() ou qualquer outra função que permita