![\"Dopasowanie](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/courberaccordr1.png\")
<\/p>\n
Poni\u017cszy przyk\u0142ad wyja\u015bnia krok po kroku, jak dopasowa\u0107 krzywe do danych w R za pomoc\u0105 funkcji Poly()<\/strong> i jak okre\u015bli\u0107, kt\u00f3ra krzywa najlepiej pasuje do danych.<\/span><\/p>\n Zacznijmy od utworzenia fa\u0142szywego zbioru danych, a nast\u0119pnie utw\u00f3rz wykres rozrzutu w celu wizualizacji danych:<\/span> <\/p>\n Dopasujmy nast\u0119pnie do danych kilka modeli regresji wielomianowej i zwizualizujmy krzyw\u0105 ka\u017cdego modelu na tym samym wykresie:<\/span> <\/p>\n Aby okre\u015bli\u0107, kt\u00f3ra krzywa najlepiej pasuje do danych, mo\u017cemy spojrze\u0107 na skorygowany kwadrat R<\/a> ka\u017cdego modelu.<\/span><\/p>\n Warto\u015b\u0107 ta m\u00f3wi nam, jaki procent zmienno\u015bci zmiennej odpowiedzi mo\u017cna wyja\u015bni\u0107 za pomoc\u0105 zmiennych predykcyjnych w modelu, skorygowanych o liczb\u0119 zmiennych predykcyjnych.<\/span><\/p>\n Z wyniku wida\u0107, \u017ce model z najwy\u017cszym skorygowanym wsp\u00f3\u0142czynnikiem R-kwadrat jest wielomianem czwartego stopnia, kt\u00f3rego skorygowany wsp\u00f3\u0142czynnik R-kwadrat wynosi 0,959<\/strong> .<\/span><\/p>\n Na koniec mo\u017cemy utworzy\u0107 wykres punktowy z krzyw\u0105 modelu wielomianowego czwartego stopnia:<\/span> <\/p>\n R\u00f3wnanie dla tej linii mo\u017cemy r\u00f3wnie\u017c uzyska\u0107 za pomoc\u0105 funkcji podsumowania()<\/strong> :<\/span><\/span><\/p>\n R\u00f3wnanie krzywej jest nast\u0119puj\u0105ce:<\/span><\/p>\n y = -0,0192x 4<\/sup> + 0,7081x 3<\/sup> \u2013 8,3649x 2<\/sup> + 35,823x \u2013 26,516<\/span><\/p>\n Mo\u017cemy u\u017cy\u0107 tego r\u00f3wnania do przewidzenia warto\u015bcizmiennej odpowiedzi<\/a> na podstawie zmiennych predykcyjnych w modelu. Na przyk\u0142ad, je\u015bli x<\/em> = 4, to przewidywaliby\u015bmy, \u017ce y<\/em> = 23,34<\/strong> :<\/span><\/p>\n y = -0,0192(4) 4<\/sup> + 0,7081(4) 3<\/sup> \u2013 8,3649(4) 2<\/sup> + 35,823(4) \u2013 26,516 = 23,34<\/span><\/p>\n Wprowadzenie do regresji wielomianowej Cz\u0119sto mo\u017cesz chcie\u0107 znale\u017a\u0107 r\u00f3wnanie, kt\u00f3re najlepiej pasuje do krzywej R. Poni\u017cszy przyk\u0142ad wyja\u015bnia krok po kroku, jak dopasowa\u0107 krzywe do danych w R za pomoc\u0105 funkcji Poly() i jak okre\u015bli\u0107, kt\u00f3ra krzywa najlepiej pasuje do danych. Krok 1: Utw\u00f3rz i wizualizuj dane Zacznijmy od utworzenia fa\u0142szywego zbioru danych, a nast\u0119pnie utw\u00f3rz wykres rozrzutu w […]<\/p>\n","protected":false},"author":1,"featured_media":0,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[3],"tags":[],"class_list":["post-1640","post","type-post","status-publish","format-standard","hentry","category-przewodnik"],"yoast_head":"\n Krok 1: Utw\u00f3rz i wizualizuj dane<\/strong><\/span><\/h3>\n
#create data frame\n<\/span>df <- data. frame<\/span> (x=1:15,\n y=c(3, 14, 23, 25, 23, 15, 9, 5, 9, 13, 17, 24, 32, 36, 46))\n\n#create a scatterplot of x vs. y\n<\/span>plot(df$x, df$y, pch= 19<\/span> , xlab=' x<\/span> ', ylab=' y<\/span> ')<\/strong> <\/pre>\n![\"\"](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/courbeajustementr0.png\")
<\/p>\n Krok 2: Dostosuj wiele krzywych<\/strong><\/span><\/h3>\n
#fit polynomial regression models up to degree 5\n<\/span>fit1 <- lm(y~x, data=df)\nfit2 <- lm(y~poly(x,2,raw= TRUE<\/span> ), data=df)\nfit3 <- lm(y~poly(x,3,raw= TRUE<\/span> ), data=df)\nfit4 <- lm(y~poly(x,4,raw= TRUE<\/span> ), data=df)\nfit5 <- lm(y~poly(x,5,raw= TRUE<\/span> ), data=df)<\/span>\n\n#create a scatterplot of x vs. y\n<\/span>plot(df$x, df$y, pch=19, xlab=' x<\/span> ', ylab=' y<\/span> ')\n\n#define x-axis values\n<\/span>x_axis <- seq(1, 15, length= 15<\/span> )\n\n#add curve of each model to plot\n<\/span>lines(x_axis, predict(fit1, data. frame<\/span> (x=x_axis)), col=' green<\/span> ')\nlines(x_axis, predict(fit2, data. frame<\/span> (x=x_axis)), col=' red<\/span> ')\nlines(x_axis, predict(fit3, data. frame<\/span> (x=x_axis)), col=' purple<\/span> ')\nlines(x_axis, predict(fit4, data. frame<\/span> (x=x_axis)), col=' blue<\/span> ')\nlines(x_axis, predict(fit5, data. frame<\/span> (x=x_axis)), col=' orange<\/span> ')\n<\/strong><\/pre>\n![\"\"](\"https:\/\/statorials.org\/wp-content\/uploads\/2023\/08\/courbeajustementr2.png\")
<\/p>\n #calculated adjusted R-squared of each model\n<\/span>summary(fit1)$adj. r<\/span> . squared<\/span>\nsummary(fit2)$adj. r<\/span> . squared<\/span>\nsummary(fit3)$adj. r<\/span> . squared<\/span>\nsummary(fit4)$adj. r<\/span> . squared<\/span>\nsummary(fit5)$adj. r<\/span> . squared<\/span>\n\n[1] 0.3144819\n[1] 0.5186706\n[1] 0.7842864\n[1] 0.9590276\n[1] 0.9549709\n<\/strong><\/pre>\n Krok 3: Wizualizuj ostateczn\u0105 krzyw\u0105<\/strong><\/span><\/h3>\n
#create a scatterplot of x vs. y\n<\/span>plot(df$x, df$y, pch=19, xlab=' x<\/span> ', ylab=' y<\/span> ')\n\n#define x-axis values\n<\/span>x_axis <- seq(1, 15, length= 15<\/span> )\n\n#add curve of fourth-degree polynomial model\n<\/span>lines(x_axis, predict(fit4, data. frame<\/span> (x=x_axis)), col=' blue<\/span> ')\n<\/strong><\/pre>\n<\/p>\n summary(fit4)\n\nCall:\nlm(formula = y ~ poly(x, 4, raw = TRUE), data = df)\n\nResiduals:\n Min 1Q Median 3Q Max \n-3.4490 -1.1732 0.6023 1.4899 3.0351 \n\nCoefficients:\n Estimate Std. Error t value Pr(>|t|) \n(Intercept) -26.51615 4.94555 -5.362 0.000318 ***\npoly(x, 4, raw = TRUE)1 35.82311 3.98204 8.996 4.15e-06 ***\npoly(x, 4, raw = TRUE)2 -8.36486 0.96791 -8.642 5.95e-06 ***\npoly(x, 4, raw = TRUE)3 0.70812 0.08954 7.908 1.30e-05 ***\npoly(x, 4, raw = TRUE)4 -0.01924 0.00278 -6.922 4.08e-05 ***\n---\nSignificant. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1\n\nResidual standard error: 2.424 on 10 degrees of freedom\nMultiple R-squared: 0.9707, Adjusted R-squared: 0.959 \nF-statistic: 82.92 on 4 and 10 DF, p-value: 1.257e-07\n<\/strong><\/pre>\n
Dodatkowe zasoby<\/strong><\/span><\/h3>\n
<\/a> Regresja wielomianowa w R (krok po kroku)<\/a>
Jak korzysta\u0107 z funkcji seq w R<\/a><\/p>\n","protected":false},"excerpt":{"rendered":"