Tenho o seguinte conjunto de dados que estabelece uma relação entre duas variáveis “X” e “Y”:
df <- data.frame(X=c(25,25,25,25,25,25,50,50,50,50,50,50,
75,75,75,75,75,75,100,100,100,100,100,100,
125,125,125,125,125,125,150,150,150,150,150,150),
Y=c(2457524,2391693,2450828,2391252,2444638,2360293,
4693194,4844527,4835596,4878092,4809226,4722253,
7142763,7182769,7135550,7173920,7216871,7076359,
9496553,9537788,9405825,9439201,9609870,9707734,
12031958,12027037,11935594,11930086,12154132,
12096462,14298064,14396607,13964716,14221039,
14283992,14042220))
Considere o seguinte problema:
“Ajustar um modelo linear ponderado usando a função “lm” e, como fator de ponderação, o inverso da variância de “Y” para cada nível de “X””. Ou seja, o modelo linear deverá ser ponderado pelo inverso da variância de cada nível de “X”. Neste caso, como podemos especificar a relação funcional ponderada? Existe alguma função específica para ser entrada como argumento em “weights”?
Detalhe técnico: Só vale ajustar pela função “lm”. Não vale ajustar por nenhum outro método (gls, glm, etc.).