###
9 / 44
0.6 * 0.4 + 0.3 * 0.6
log(0.6 * 0.4 + 0.3 * 0.6)
9 / 44
0.6 * 0.4 + 0.3 * 0.6
log(0.6 * 0.4 + 0.3 * 0.6)
##
Scores = c(40, 39, 15, 6, 18, 22, 30, 21, 15, 23)
Scores

meancomputation = (40+ 39+ 15+ 6+ 18+ 22+ 30+ 21+ 15+ 23)/10
meancomputation
meaninbuilt = mean(Scores) 
meaninbuilt

Scores2 = Scores # create a copy of Scores
Scores2[4] = 16
Scores2
Scores[c(1,3,5)]

Scores
y = which(Scores == 30) 
y

z = which(Scores <= 20)
z

#####
x = 1:100
x
x[x < 10 | x > 90]

####
data()
?airquality
head(airquality)
tail(airquality)
head(airquality, n = 10)
airquality[148,4]
airquality$Temp[148]
airquality[148,]
airquality[,c(1,4)]
summary(airquality$Temp)
hist(airquality$Temp)
plot(airquality$Temp)
plot(airquality$Ozone, airquality$Temp)
plot(airquality)

####
install.packages("UsingR")
library("UsingR")
install.packages("tidyverse")

####

library("tidyverse")
library("viridis")
mpg
ggplot(data=mpg) + 
    geom_point(mapping=aes(x=displ, y=hwy))

ggplot(data=mpg) + 
    geom_point(mapping=aes(x=displ, y=hwy, colour=class))

ggplot(data=mpg) + 
    geom_point(mapping=aes(x=displ, y=hwy, colour=class))+
    scale_colour_viridis_d()

ggplot(data=mpg) + 
    geom_point(mapping=aes(x=displ, y=hwy, colour=class))+
   scale_colour_viridis_d(option = "plasma")

ggplot(data=mpg) + 
    geom_point(mapping=aes(x=displ, y=hwy))+ 
                   facet_wrap(~ class, nrow=2)

ggplot(data=mpg) + 
    geom_point(mapping=aes(x=displ, y=hwy, colour=class))+
   scale_colour_viridis_d(option = "plasma")+
    facet_wrap(~class, nrow=2)

ggplot(data=mpg) + 
    geom_point(mapping=aes(x=displ, y=hwy, colour=class))+
   scale_colour_viridis_d(option = "plasma")+
    facet_grid(drv~cyl)

ggplot(data=mpg) + 
    geom_bar(mapping=aes(x=displ, colour=class),position = "dodge")+  scale_colour_viridis_d(option = "plasma") 

ggplot(data = mpg, mapping = aes(x = displ, y = hwy)) +
  geom_smooth()

ggplot(data = mpg, mapping = aes(x = class, y = hwy)) + 
  geom_boxplot() +
  coord_flip() 

ggplot(data=mpg) + 
    geom_point(mapping=aes(x=displ, y=hwy, colour=class))+
   scale_colour_viridis_d(option = "plasma")

ggplot(data = mpg) + 
  geom_bar(mapping = aes(x = class))

ggplot(data = mpg) + 
  stat_count(mapping = aes(x = class))

table= as.data.frame(table(mpg$class))
ggplot(data = table) +   
geom_bar(mapping = aes(x = Var1, y=Freq), 
  stat="identity")


ggplot(data = mpg) + 
  stat_summary(
    mapping = aes(x = class, y = hwy),
    fun.min = min,
    fun.max = max,
    fun = median
  )

ggplot(data = mpg) +   
geom_bar(mapping = aes(x = class, colour = class)) 

ggplot(data = mpg) +
  geom_bar(mapping = aes(x = class, fill = class)) + 
  scale_fill_viridis_d()

ggplot(data = mpg) +
  geom_bar(mapping = aes(x = class, fill = class)) + 
  scale_fill_viridis_d() + coord_flip()

ggplot(data = mpg) +
  geom_bar(mapping = aes(x = class, fill = class)) + 
  scale_fill_viridis_d() + coord_polar()

ggplot(data = mpg) +
  geom_bar(mapping = aes(x = class, fill = trans)) + 
  scale_fill_viridis_d()

ggplot(data = mpg) +
  geom_bar(mapping = aes(x = class, fill = trans),
           position="dodge") + 
  scale_fill_viridis_d()