Se vc tiver interesse em salvar o código tex desse gráfico pode usar o pacote
tikz. Ele vai gerar um arquivo tex que vc pode dar \input no latex. Não testei
essa solução, é apenas um caminho.
Sent from my iPhone
___
R-br mailing list
Obrigado Walmes, mas os contrastes não pareceram ficar muito corretos,
sendo meu CRM corrigido:
##-
require(multcomp)
require(doBy)
require(latticeExtra)
require(wzRfun)
É uma questão complicada de ser comparada. Por exemplo, se um pesquisador
consegue trabalhar com:
Latex
R +
tikz = gerar códigos tex das figuras
tables = gerar tabelas com dados que estão no R
Hmisc::latex = converter as tabelas geradas com o tables em tex
ou knitr
e consegue ainda editar tudo
WALMES, muito obrigado!!
Agora a saída do desdobramento está correta.
adaptado de Walmes ---# Fatiar solo/nc.
m1 - aov(mspa~bloco+np*solo + solo/nc + np:solo:nc, data=a1)
anova(m1)
m1$assign
names(coef(m1))[m1$assign==5] # ir mudando o nº ate achar a inter solo:nc
# buscando as
Bem, talvez não seja o seu caso, mas se esse for um script longo e isso
estiver acontecendo lá do meio em diante é sinal de que algo abriu e não
fechou em algum lugar acima, um ( [ { '. Tente encontrar e
fechar/corrigir. A indentação volta ao normal.
À disposição.
Walmes.
Segue exemplos que podem ser úteis.
##-
## CMR.
## Fatorial 4x5 em parcela subdividida.
## Calagem na parcela, adubação na subparcela.
## DIC.
da - expand.grid(cal=gl(4,1), adub=gl(5,1), rept=gl(4,1))
da$y -
Você pode obter o quadro de anova com SQ fatiadas para solo/nc e com com a
linha do resíduo do modelo completo, o seu avi1. Para isso você não pode
abandonar nenhum termo. Veja o exemplo.
da - expand.grid(bloc=gl(4,1),
A=gl(2,1), B=gl(3,1), C=gl(4,1))
da$y - rnorm(nrow(da))
m0
Alexandre,
No seu CMR tem-se um fator apenas. Você não precisa de a e a.int, eles
são iguais. Será que não é a versão da apc que você tem? Ela tem sido
aprimorada, está no meu pacote wzRfun, disponível do
www.github.com/walmes/wzRfun. A mais recente tá assim:
## wzRfun::apc
apc - function (lfm,