Ajustando o modelo polinomial aos dados em R

Li as respostas a esta pergunta e são bastante úteis, mas preciso de ajuda principalmente na R. Tenho um exemplo de conjunto de dados em R da seguinte forma: x <- c(32,64,96,118,126,144,152.5,158) y <- c(99.5,104.8,108.5,100,86,64,35.3,15) Quero ajustar um modelo a esses dados para...