Colegas,

Tenho um conjunto de dados conforme modelo abaixo:


Genotipo        Tratamento      Repeticao       Composto1       Composto2       Composto3       Composto4
C13     Resistente      SI      R1      0,04    8,77    0
C13     Resistente      SI      R2      0       8,67    0
C13     Resistente      SI      R3      0       7,99    0
C13     Resistente      CI      R1      0,04    33,83   0
C13     Resistente      CI      R2      0,02    6,68    0
C13     Resistente      CI      R3      0       7,32    0
C08     Resistente      SI      R1      0,01    30,48   0,04
C08     Resistente      SI      R2      0       22,83   0
C08     Resistente      SI      R3      0       30,66   0
C08     Resistente      CI      R1      0       30,19   0
C08     Resistente      CI      R2      0       26,69   0
C08     Resistente      CI      R3      0       33,55   0
C17     Suscetível      SI      R1      0,08    16,94   0
C17     Suscetível      SI      R2      0,06    22,3    0
C17     Suscetível      SI      R3      0,07    17,19   0
C17     Suscetível      CI      R1      0,05    17,72   0
C17     Suscetível      CI      R2      0,05    19,46   0
C17     Suscetível      CI      R3      0,13    30,15   0
C11     Suscetível      SI      R1      0,04    30,03   0
C11     Suscetível      SI      R2      0       38,64   0
C11     Suscetível      SI      R3      0,06    42,59   0
C11     Suscetível      CI      R1      0,04    36,96   0
C11     Suscetível      CI      R2      0       49,7    0
C11     Suscetível      CI      R3      0       37,67   0,06



Esse exemplo está disponível no endereço:
https://www.dropbox.com/s/izccqoifwnjkdlu/Exemplo.csv

Genotipo - diferentes genótipos, incluindo plantas resistentes e
suscetíveis

Tratamento - SI = sem inoculação; CI = com inoculação

Repeticao - R1 = planta 1; R2 = planta 2; R3 = planta 3

Composto1 a Composton = valor aferido (medido) nas folhas de cada
planta para o respectivo composto

Interesse:

1. verificar quais compostos são produzidos em função da inoculação (CI
vs SI)

2. verificar quais compostos são produzidos em função do Estado
(Resistente vs Suscetível)

3. verificar se o Genotipo interfere na produção de determinado
composto (composto específico a um dado Genotipo)

4. verificar se os demais genótipos diferem do E. camaldulensis.

O valor 0 (zero) para um determinado composto não significa zero, mas,
significa que aquele composto não foi encontrado naquela planta
(repetição). Logo, 0 significa NA.

A minha ideia é analisar composto a composto separadamente. Tenho 117
compostos.

Modelo que tentei usar:

dados02 <- read.table(url("https://www.dropbox.com/s/izccqoifwnjkdlu/Exemplo.csv?dl=1"), sep="\t", header=TRUE, dec=",")

dados02

fit02 <- lm(Composto1 ~ Genotipo * Estado * Tratamento, data=dados02)

summary(fit02)

Aparece o seguinte erro:

Coefficients: (9 not defined because of singularities)

Pelo que li, parece que a variável é colinear ou possui correlação. Não
sei como resolver.

Outras perguntas:

A minha abordagem está apropriada? Terei que rodar as 117?

Existe uma maneira mais adequada de responder às quatro perguntas
acima? Principalmente a 1 e 2? Gráficos?

Obrigado!

--
Marcelo
_______________________________________________
R-br mailing list
R-br@listas.c3sl.ufpr.br
https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br
Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.