R Kodu (Örnek)
1 R Markdown
Bu bir R Markdown dokümanıdır. Basit bir formatlama sentaksı olan Markdown, HTML, PDF ve MS Word dokümanlarının hazırlanmasında kullanılabilir. Daha fazla bilgi için lütfen linki inceleyin http://rmarkdown.rstudio.com.
2 R Markdown’da Kod Yığını
R Markdown’da, R kod yığınını aşağıdaki gibi yerleştirebilirsiniz: Örnek olarak ortalaması 0, standard sapması 1 ve normal dağılım yapan bir ana kütleden 200 gözlemli bir örneklem çekip, grafiğini çıkartalım.
x <- rnorm(200, mean = 0, sd = 1)
par(mar = c(2, 2, 0.1, 2))
plot(x, dnorm(x))
3 Basit Doğrusal Regresyon Örneği
R’ın içindeki cars
datasının ilk 10 gözlemini kullanarak basit doğrusal bir regresyon oluşturalım…
data <- cars[1:10, ] ## Yeni data seti data olarak adlandırıldı.
summary(data)
#> speed dist
#> Min. : 4.0 Min. : 2.0
#> 1st Qu.: 7.0 1st Qu.:10.0
#> Median : 8.5 Median :16.5
#> Mean : 8.0 Mean :15.9
#> 3rd Qu.:10.0 3rd Qu.:21.0
#> Max. :11.0 Max. :34.0
model <- lm(dist ~ speed, data = data)
summary(model)
#>
#> Call:
#> lm(formula = dist ~ speed, data = data)
#>
#> Residuals:
#> Min 1Q Median 3Q Max
#> -9.3464 -5.8420 -1.4536 4.8232 12.9929
#>
#> Coefficients:
#> Estimate Std. Error t value Pr(>|t|)
#> (Intercept) -4.5286 8.9164 -0.5079 0.62523
#> speed 2.5536 1.0688 2.3893 0.04391 *
#> ---
#> Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
#>
#> Residual standard error: 7.998 on 8 degrees of freedom
#> Multiple R-squared: 0.41642, Adjusted R-squared: 0.34348
#> F-statistic: 5.7085 on 1 and 8 DF, p-value: 0.043907
3.1 Basit Doğrusal Regresyon: Basit Grafik Örneği
Basit doğrusal regresyonun sonuçlarını kullanarak Örneklem Regresyon Fonksiyonu grafiğini basit bir kod kullanarak gösterelim…
par(mar = c(2, 2, 0.1, 2))
plot(data$speed, data$dist)
abline(model)
3.2 Basit Doğrusal Regresyon: Detaylı Grafik Örneği
Daha önce basit bir kod kullanarak Grafik 3.1’de gösterdiğimiz Örneklem Regresyon Fonksiyonu grafiğini şimdi de daha detaylı bir kod kullanarak çizelim…
par(mar = c(2, 2, 1.5, 2))
plot(data$speed, data$dist, type = "p", pch = 16, cex = 1, col = "blue", main = "Örneklem Regresyon Fonksiyonu", xlab = "X", ylab = "Y", xlim = c(0, max(data$speed)), ylim = c(0, max(data$dist)))
abline(model, lwd = 2, lty = 2, col = "red")
text(data$dist ~ data$speed, labels = rownames(data), cex = 1, font = 2, pos = 4)
cf <- round(coef(model), 2)
mtext(latex2exp::TeX(paste("$", "y", "=", cf[[1]], ifelse(sign(cf[[2]]) == 1, " + ", " - "), abs(cf[[2]]), "x", "$")), adj = 1, padj = 0)
4 R Video ve Kodları
Akademi Ekonometri olarak YouTube kanalımızda paylaştığımız R videolarını ve bu videoları hazırlarken kullandığımız data ve R kodlarını da paylaşarak sizleri olabildiğince desteklemek ve kodlamaya teşvik etmeyi amaçlıyoruz. Daha önce yayınladığımız örnek bir videoyu aşağıda inceleyebilirsiniz.
Ayrıca bu videoya ait data ve R kodlarına aşağıdaki linklere tıklayarak erişebilirsiniz.
- test.xls
ve R
kodu
- CPI.csv
ve R
kodu