library(dplyr)
library(readr)
library(ggplot2)
?mtcars
head(mtcars)
tail(mtcars)
summary(mtcars)
      mpg             cyl             disp             hp             drat      
 Min.   :10.40   Min.   :4.000   Min.   : 71.1   Min.   : 52.0   Min.   :2.760  
 1st Qu.:15.43   1st Qu.:4.000   1st Qu.:120.8   1st Qu.: 96.5   1st Qu.:3.080  
 Median :19.20   Median :6.000   Median :196.3   Median :123.0   Median :3.695  
 Mean   :20.09   Mean   :6.188   Mean   :230.7   Mean   :146.7   Mean   :3.597  
 3rd Qu.:22.80   3rd Qu.:8.000   3rd Qu.:326.0   3rd Qu.:180.0   3rd Qu.:3.920  
 Max.   :33.90   Max.   :8.000   Max.   :472.0   Max.   :335.0   Max.   :4.930  
       wt             qsec             vs               am              gear      
 Min.   :1.513   Min.   :14.50   Min.   :0.0000   Min.   :0.0000   Min.   :3.000  
 1st Qu.:2.581   1st Qu.:16.89   1st Qu.:0.0000   1st Qu.:0.0000   1st Qu.:3.000  
 Median :3.325   Median :17.71   Median :0.0000   Median :0.0000   Median :4.000  
 Mean   :3.217   Mean   :17.85   Mean   :0.4375   Mean   :0.4062   Mean   :3.688  
 3rd Qu.:3.610   3rd Qu.:18.90   3rd Qu.:1.0000   3rd Qu.:1.0000   3rd Qu.:4.000  
 Max.   :5.424   Max.   :22.90   Max.   :1.0000   Max.   :1.0000   Max.   :5.000  
      carb       cylFactor
 Min.   :1.000   4:11     
 1st Qu.:2.000   6: 7     
 Median :2.000   8:14     
 Mean   :2.812            
 3rd Qu.:4.000            
 Max.   :8.000            
mean(mtcars$cyl)
[1] 6.1875

Scatter Plot

# scatter plot shape 1
ggplot(mtcars, aes(x=mpg, y=wt)) + 
  geom_point(shape=1)

# scatter plot shape 19
ggplot(mtcars, aes(x=mpg, y=wt)) + 
  geom_point(shape=19)

# shape identity
ggplot(mtcars,aes(x=mpg, y=wt, shape=cyl)) + 
  geom_point() + 
  scale_shape_identity()

mtcars$cylFactor <- factor(mtcars$cyl)
ggplot(mtcars,aes(x=mpg, y=wt, shape=cylFactor)) + 
  geom_point()

# color parameter
ggplot(mtcars, aes(x=mpg, y=wt)) + 
  geom_point(shape=19, color="blue")

# shape 19
ggplot(mtcars, aes(x=mpg, y=wt, color=cyl)) + 
  geom_point(shape=19)

ggplot(mtcars,aes(x=mpg, y=wt, color=cylFactor)) + 
  geom_point(shape=19)

# rename the legend
ggplot(mtcars, aes(x=mpg, y=wt, color=cylFactor)) + 
  geom_point(shape=19) + 
  labs(color="Cylinders")

# add new labels and title
ggplot(mtcars, aes(x=mpg, y=wt, color=cylFactor)) + 
  geom_point(shape=19) +
  xlab("Miles per Gallon") +
  ylab("Weight") +
  labs(color="Cylinders") +
  ggtitle("Scatterplot")

Line Graph EuStockMarkets Dataset

?EuStockMarkets
EuStockDF <- as.data.frame(EuStockMarkets)
head(EuStockDF)
# DAX stock price
ggplot(EuStockDF, aes(x=c(1:nrow(EuStockDF)), y=DAX)) + 
  geom_line()

# line width
ggplot(EuStockDF, aes(x=c(1:nrow(EuStockDF)), y=DAX)) + 
  geom_line(size=1.5) +
  labs(x="Stocks")

# color
ggplot(EuStockDF,aes(x=c(1:nrow(EuStockDF)), y=DAX)) + 
  geom_line(size=1.5, colour="light blue") + 
  labs(x="Time", y="Stocks")

# plot the 2 stocks
dax_smi_plot <- ggplot() +
  geom_line(data=EuStockDF, aes(x=c(1:nrow(EuStockDF)), y=DAX), size=1.5, colour="light blue") +
  geom_line(data=EuStockDF, aes(x=c(1:nrow(EuStockDF)), y=SMI), size=1.5, colour="red") +
  labs(x="Time", y="Stocks")
print(dax_smi_plot)

# all stocks
all_stocks <- ggplot() +
  geom_line(data = EuStockDF, aes(x=c(1:nrow(EuStockDF)), y=DAX), size=1, colour="light blue") +
  geom_line(data = EuStockDF, aes(x=c(1:nrow(EuStockDF)), y=SMI), size=1, colour="red") + 
  geom_line(data = EuStockDF, aes(x=c(1:nrow(EuStockDF)), y=CAC), size=1, colour="purple") + 
  geom_line(data = EuStockDF, aes(x=c(1:nrow(EuStockDF)), y=FTSE), size=1, colour="green") +
  labs(x="Time", y="Stocks")
print(all_stocks)

# rename the labels
legend_stocks <- all_stocks + 
  xlab("Days") + 
  ylab("Price") + 
  ggtitle("Eu Stocks")
print(legend_stocks)

Linear Regression

ggplot(mtcars, aes(x=mpg, y=wt))  + 
  geom_point(shape=19) +
  geom_smooth(method="lm", se=FALSE, color="red")

#se = TRUE -> confidence interval appear (default = true)

ggplot(mtcars, aes(x=mpg, y=wt, color=cylFactor)) + 
  geom_point(shape=19) +
  geom_smooth(method="lm", se=TRUE, color="red")

# add the labels
ggplot(mtcars, aes(x=mpg, y=wt, color=cylFactor)) + 
  geom_point(shape=19) + 
  geom_smooth(method="lm", se=TRUE, color="red") + 
  xlab("Miles per Gallon") + 
  ylab("Weight") +
  labs(color="Cylinders") +
  ggtitle("Linear Regression")

Gaussian Regression

ggplot(mtcars, aes(x=mpg, y=wt, color=cylFactor)) +
  geom_point(shape=19) + 
  geom_smooth(method="auto", se=TRUE, color="red") +
  xlab("Miles per Gallon") + 
  ylab("Weight") +
  labs(colour="Cylinders") + 
  ggtitle("Gaussian Regression")

End Of Report

LS0tCnRpdGxlOiAiREFUQSBWSVNVQUxJWkFUSU9OIFBST0pFQ1QiCmF1dGhvcjogIk5pY2EgViAoaHR0cHM6Ly93d3cuY3JlZGx5LmNvbS91c2Vycy9uaWNhLXYpIgpvdXRwdXQ6IGh0bWxfbm90ZWJvb2sKLS0tCgoKCgpgYGB7cn0KbGlicmFyeShkcGx5cikKbGlicmFyeShyZWFkcikKbGlicmFyeShnZ3Bsb3QyKQpgYGAKCgoKYGBge3J9Cj9tdGNhcnMKYGBgCgoKCmBgYHtyfQpoZWFkKG10Y2FycykKYGBgCgoKCmBgYHtyfQp0YWlsKG10Y2FycykKYGBgCgoKCmBgYHtyfQpzdW1tYXJ5KG10Y2FycykKYGBgCgoKCmBgYHtyfQptZWFuKG10Y2FycyRjeWwpCmBgYAoKCgojICoqU2NhdHRlciBQbG90KioKCgoKYGBge3J9CiMgc2NhdHRlciBwbG90IHNoYXBlIDEKZ2dwbG90KG10Y2FycywgYWVzKHg9bXBnLCB5PXd0KSkgKyAKICBnZW9tX3BvaW50KHNoYXBlPTEpCmBgYAoKCgoKYGBge3J9CiMgc2NhdHRlciBwbG90IHNoYXBlIDE5CmdncGxvdChtdGNhcnMsIGFlcyh4PW1wZywgeT13dCkpICsgCiAgZ2VvbV9wb2ludChzaGFwZT0xOSkKYGBgCgoKCgpgYGB7cn0KIyBzaGFwZSBpZGVudGl0eQpnZ3Bsb3QobXRjYXJzLGFlcyh4PW1wZywgeT13dCwgc2hhcGU9Y3lsKSkgKyAKICBnZW9tX3BvaW50KCkgKyAKICBzY2FsZV9zaGFwZV9pZGVudGl0eSgpCgpgYGAKCgoKCmBgYHtyfQptdGNhcnMkY3lsRmFjdG9yIDwtIGZhY3RvcihtdGNhcnMkY3lsKQpgYGAKCgoKCmBgYHtyfQpnZ3Bsb3QobXRjYXJzLGFlcyh4PW1wZywgeT13dCwgc2hhcGU9Y3lsRmFjdG9yKSkgKyAKICBnZW9tX3BvaW50KCkKYGBgCgoKCgpgYGB7cn0KIyBjb2xvciBwYXJhbWV0ZXIKZ2dwbG90KG10Y2FycywgYWVzKHg9bXBnLCB5PXd0KSkgKyAKICBnZW9tX3BvaW50KHNoYXBlPTE5LCBjb2xvcj0iYmx1ZSIpCmBgYAoKCgoKCmBgYHtyfQojIHNoYXBlIDE5CmdncGxvdChtdGNhcnMsIGFlcyh4PW1wZywgeT13dCwgY29sb3I9Y3lsKSkgKyAKICBnZW9tX3BvaW50KHNoYXBlPTE5KQpgYGAKCgoKCgpgYGB7cn0KZ2dwbG90KG10Y2FycyxhZXMoeD1tcGcsIHk9d3QsIGNvbG9yPWN5bEZhY3RvcikpICsgCiAgZ2VvbV9wb2ludChzaGFwZT0xOSkKYGBgCgoKCgoKYGBge3J9CiMgcmVuYW1lIHRoZSBsZWdlbmQKZ2dwbG90KG10Y2FycywgYWVzKHg9bXBnLCB5PXd0LCBjb2xvcj1jeWxGYWN0b3IpKSArIAogIGdlb21fcG9pbnQoc2hhcGU9MTkpICsgCiAgbGFicyhjb2xvcj0iQ3lsaW5kZXJzIikKYGBgCgoKCgoKYGBge3J9CiMgYWRkIG5ldyBsYWJlbHMgYW5kIHRpdGxlCmdncGxvdChtdGNhcnMsIGFlcyh4PW1wZywgeT13dCwgY29sb3I9Y3lsRmFjdG9yKSkgKyAKICBnZW9tX3BvaW50KHNoYXBlPTE5KSArCiAgeGxhYigiTWlsZXMgcGVyIEdhbGxvbiIpICsKICB5bGFiKCJXZWlnaHQiKSArCiAgbGFicyhjb2xvcj0iQ3lsaW5kZXJzIikgKwogIGdndGl0bGUoIlNjYXR0ZXJwbG90IikKYGBgCgoKCgojICoqTGluZSBHcmFwaCBFdVN0b2NrTWFya2V0cyBEYXRhc2V0KioKCgoKYGBge3J9Cj9FdVN0b2NrTWFya2V0cwpgYGAKCgoKYGBge3J9CkV1U3RvY2tERiA8LSBhcy5kYXRhLmZyYW1lKEV1U3RvY2tNYXJrZXRzKQpgYGAKCgoKYGBge3J9CmhlYWQoRXVTdG9ja0RGKQpgYGAKCgoKCgpgYGB7cn0KIyBEQVggc3RvY2sgcHJpY2UKZ2dwbG90KEV1U3RvY2tERiwgYWVzKHg9YygxOm5yb3coRXVTdG9ja0RGKSksIHk9REFYKSkgKyAKICBnZW9tX2xpbmUoKQpgYGAKCgoKCgpgYGB7cn0KIyBsaW5lIHdpZHRoCmdncGxvdChFdVN0b2NrREYsIGFlcyh4PWMoMTpucm93KEV1U3RvY2tERikpLCB5PURBWCkpICsgCiAgZ2VvbV9saW5lKHNpemU9MS41KSArCiAgbGFicyh4PSJTdG9ja3MiKQpgYGAKCgoKCgpgYGB7cn0KIyBjb2xvcgpnZ3Bsb3QoRXVTdG9ja0RGLGFlcyh4PWMoMTpucm93KEV1U3RvY2tERikpLCB5PURBWCkpICsgCiAgZ2VvbV9saW5lKHNpemU9MS41LCBjb2xvdXI9ImxpZ2h0IGJsdWUiKSArIAogIGxhYnMoeD0iVGltZSIsIHk9IlN0b2NrcyIpCmBgYAoKCgoKCmBgYHtyfQojIHBsb3QgdGhlIDIgc3RvY2tzCmRheF9zbWlfcGxvdCA8LSBnZ3Bsb3QoKSArCiAgZ2VvbV9saW5lKGRhdGE9RXVTdG9ja0RGLCBhZXMoeD1jKDE6bnJvdyhFdVN0b2NrREYpKSwgeT1EQVgpLCBzaXplPTEuNSwgY29sb3VyPSJsaWdodCBibHVlIikgKwogIGdlb21fbGluZShkYXRhPUV1U3RvY2tERiwgYWVzKHg9YygxOm5yb3coRXVTdG9ja0RGKSksIHk9U01JKSwgc2l6ZT0xLjUsIGNvbG91cj0icmVkIikgKwogIGxhYnMoeD0iVGltZSIsIHk9IlN0b2NrcyIpCnByaW50KGRheF9zbWlfcGxvdCkKYGBgCgoKCgoKYGBge3J9CiMgYWxsIHN0b2NrcwphbGxfc3RvY2tzIDwtIGdncGxvdCgpICsKICBnZW9tX2xpbmUoZGF0YSA9IEV1U3RvY2tERiwgYWVzKHg9YygxOm5yb3coRXVTdG9ja0RGKSksIHk9REFYKSwgc2l6ZT0xLCBjb2xvdXI9ImxpZ2h0IGJsdWUiKSArCiAgZ2VvbV9saW5lKGRhdGEgPSBFdVN0b2NrREYsIGFlcyh4PWMoMTpucm93KEV1U3RvY2tERikpLCB5PVNNSSksIHNpemU9MSwgY29sb3VyPSJyZWQiKSArIAogIGdlb21fbGluZShkYXRhID0gRXVTdG9ja0RGLCBhZXMoeD1jKDE6bnJvdyhFdVN0b2NrREYpKSwgeT1DQUMpLCBzaXplPTEsIGNvbG91cj0icHVycGxlIikgKyAKICBnZW9tX2xpbmUoZGF0YSA9IEV1U3RvY2tERiwgYWVzKHg9YygxOm5yb3coRXVTdG9ja0RGKSksIHk9RlRTRSksIHNpemU9MSwgY29sb3VyPSJncmVlbiIpICsKICBsYWJzKHg9IlRpbWUiLCB5PSJTdG9ja3MiKQpwcmludChhbGxfc3RvY2tzKQpgYGAKCgoKCmBgYHtyfQojIHJlbmFtZSB0aGUgbGFiZWxzCmxlZ2VuZF9zdG9ja3MgPC0gYWxsX3N0b2NrcyArIAogIHhsYWIoIkRheXMiKSArIAogIHlsYWIoIlByaWNlIikgKyAKICBnZ3RpdGxlKCJFdSBTdG9ja3MiKQpwcmludChsZWdlbmRfc3RvY2tzKQpgYGAKCgoKCiMgKipMaW5lYXIgUmVncmVzc2lvbioqCgoKCmBgYHtyfQpnZ3Bsb3QobXRjYXJzLCBhZXMoeD1tcGcsIHk9d3QpKSAgKyAKICBnZW9tX3BvaW50KHNoYXBlPTE5KSArCiAgZ2VvbV9zbW9vdGgobWV0aG9kPSJsbSIsIHNlPUZBTFNFLCBjb2xvcj0icmVkIikKYGBgCgoKCgoKI3NlID0gVFJVRSAgLT4gY29uZmlkZW5jZSBpbnRlcnZhbCBhcHBlYXIgKGRlZmF1bHQgPSB0cnVlKQoKCgpgYGB7cn0KZ2dwbG90KG10Y2FycywgYWVzKHg9bXBnLCB5PXd0LCBjb2xvcj1jeWxGYWN0b3IpKSArIAogIGdlb21fcG9pbnQoc2hhcGU9MTkpICsKICBnZW9tX3Ntb290aChtZXRob2Q9ImxtIiwgc2U9VFJVRSwgY29sb3I9InJlZCIpCgpgYGAKCgoKCmBgYHtyfQojIGFkZCB0aGUgbGFiZWxzCmdncGxvdChtdGNhcnMsIGFlcyh4PW1wZywgeT13dCwgY29sb3I9Y3lsRmFjdG9yKSkgKyAKICBnZW9tX3BvaW50KHNoYXBlPTE5KSArIAogIGdlb21fc21vb3RoKG1ldGhvZD0ibG0iLCBzZT1UUlVFLCBjb2xvcj0icmVkIikgKyAKICB4bGFiKCJNaWxlcyBwZXIgR2FsbG9uIikgKyAKICB5bGFiKCJXZWlnaHQiKSArCiAgbGFicyhjb2xvcj0iQ3lsaW5kZXJzIikgKwogIGdndGl0bGUoIkxpbmVhciBSZWdyZXNzaW9uIikKYGBgCgoKCiMgKipHYXVzc2lhbiBSZWdyZXNzaW9uKioKCgpgYGB7cn0KZ2dwbG90KG10Y2FycywgYWVzKHg9bXBnLCB5PXd0LCBjb2xvcj1jeWxGYWN0b3IpKSArCiAgZ2VvbV9wb2ludChzaGFwZT0xOSkgKyAKICBnZW9tX3Ntb290aChtZXRob2Q9ImF1dG8iLCBzZT1UUlVFLCBjb2xvcj0icmVkIikgKwogIHhsYWIoIk1pbGVzIHBlciBHYWxsb24iKSArIAogIHlsYWIoIldlaWdodCIpICsKICBsYWJzKGNvbG9yPSJDeWxpbmRlcnMiKSArIAogIGdndGl0bGUoIkdhdXNzaWFuIFJlZ3Jlc3Npb24iKQpgYGAKCgoKIyAqKkVuZCBPZiBSZXBvcnQqKgoK