Gurus do R,
Tenho o seguinte data frame (Df) que estabelece a relação entre as variáveis X e Y:
X Y
1 25 2457524
2 25 2391693
3 25 2450828
4 25 2391252
5 25 2444638
6 25 2360293
7 50 4693194
8 50 4844527
9 50 4835596
10 50 4878092
11 50 4809226
12 50 4722253
13 75 7142763
14 75 7182769
15 75 7135550
16 75 7173920
17 75 7216871
18 75 7076359
19 100 9496553
20 100 9537788
21 100 9405825
22 100 9439201
23 100 9609870
24 100 9707734
25 125 12031958
26 125 12027037
27 125 11935594
28 125 11930086
29 125 12154132
30 125 12096462
31 150 14298064
32 150 14396607
33 150 13964716
34 150 14221039
35 150 14283992
36 150 14042220
(Note que temos 7 níveis para a variável X com 6 pontos em cada nível)
Se ajustamos um modelo polinomial de 2º grau para estes dados obteremos o seguinte modelo:
Model<-lm(formula = Y ~ X + I(X^2))
print(Model)
Call:
lm(formula = Y ~ X + I(X^2))
Coefficients:
(Intercept) X I(X^2)
-26588.12 97310.61 -14.02
A representação gráfica deste modelo, que mais parece uma reta, é a seguinte:
Se quisermos utilizar o modelo para fazer a predição dos valores de "Y" a partir dos valores da variável "X" basta executar esta linha de código:
>predicted.intervals <- predict(Model,data.frame(x=X),interval='confidence',
+ level=0.95)
>predicted.intervals
fit lwr upr
1 2397413 2315346 2479481
2 2397413 2315346 2479481
3 2397413 2315346 2479481
4 2397413 2315346 2479481
5 2397413 2315346 2479481
6 2397413 2315346 2479481
7 4803887 4753705 4854070
8 4803887 4753705 4854070
9 4803887 4753705 4854070
10 4803887 4753705 4854070
11 4803887 4753705 4854070
12 4803887 4753705 4854070
13 7192834 7137649 7248019
14 7192834 7137649 7248019
15 7192834 7137649 7248019
16 7192834 7137649 7248019
17 7192834 7137649 7248019
18 7192834 7137649 7248019
19 9564252 9509067 9619438
20 9564252 9509067 9619438
21 9564252 9509067 9619438
22 9564252 9509067 9619438
23 9564252 9509067 9619438
24 9564252 9509067 9619438
25 11918144 11867961 11968326
26 11918144 11867961 11968326
27 11918144 11867961 11968326
28 11918144 11867961 11968326
29 11918144 11867961 11968326
30 11918144 11867961 11968326
31 14254507 14172440 14336574
32 14254507 14172440 14336574
33 14254507 14172440 14336574
34 14254507 14172440 14336574
35 14254507 14172440 14336574
36 14254507 14172440 14336574
A pergunta que não quer calar:
Qual seria o linha(s) de código para fazermos a predição inversa, ou seja, neste modelo, prever "X" a partir dos dados da variável "Y"? Pesquisando no google já tentei vários pacotes e funções específicas mas infelizmente não obtive sucesso (talvez por falta de familiaridade com as funções). Será que algum de vocês poderiam me ajudar a desvendar este mistério? Grande abraço à todos.
lm(X ~Y + I(Y^2))
e prevêX
a partir dey
diretamente? Estatisticamente, fazer isso que você quer é estranho pois X é considerada uma variável observada, não teria pq você querer prevê-la...summary(Model)
, que é uma resposta mais interessante estatisticamente do queprint(Model)
. Com ela tu vai testar as hipóteses dos três coeficientes do teu modelo serem iguais a zero. Dado o gráfico que estou vendo, eu aposto que o termo quadrático não é significante (i.e., p-valor > 0,05). Ou seja, tu tem dados que seguem um modelo linear, não quadrático.