期中作業

Statistic <- c(68,85,74,88,63,78,90,80,58,63) Math <- c(85,91,74,100,82,84,78,100,51,70)

Statistic <- c(68,85,74,88,63,78,90,80,58,63)
Math <- c(85,91,74,100,82,84,78,100,51,70)

plot(Statistic,Math,
     pch = 25,
     col= "green",
     main ="統計和數學分數",
     xlab ="統計",
     ylab ="數學")

hist(Math,
     col= "lightgreen",
     main ="數學成績",
     xlab ="分數",
     ylab ="人數")

# Load ggplot2
library(ggplot2)
# Load ggplot2
library(ggplot2)

# Create data
data <- data.frame(
  name=c("公益活動","知識閱讀","科技創新","娛樂休閒","體育競技") ,  
  value=c(25,82,28,185,36)
  )

# Barplot
ggplot(data, aes(x=name, y=value)) + 
  geom_bar(stat = "identity", width=0.2,  fill="skyblue",col="orange") 

data2<- c(25,82,28,185,36)
labels <- c("公益活動","知識閱讀","科技創新","娛樂休閒","體育競技")

pie(data2,labels,main ="喜愛社團比例", col=heat.colors(length(data2)))

# Library
library(tidyverse)

## ── Attaching core tidyverse packages ──────────────────────── tidyverse 2.0.0 ──
## ✔ dplyr     1.1.4     ✔ readr     2.1.5
## ✔ forcats   1.0.0     ✔ stringr   1.5.1
## ✔ lubridate 1.9.4     ✔ tibble    3.2.1
## ✔ purrr     1.0.4     ✔ tidyr     1.3.1
## ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ──
## ✖ dplyr::filter() masks stats::filter()
## ✖ dplyr::lag()    masks stats::lag()
## ℹ Use the conflicted package (<http://conflicted.r-lib.org/>) to force all conflicts to become errors

# plot
ggplot(data, aes(x=name, y=value)) +
  geom_segment( aes(x=name, xend=name, y=0, yend=value)) +
  geom_point( size=5, color="pink", fill=alpha("yellow", 1.0), alpha=0.7, shape=21, stroke=2)

v<-c(84,63,61,49,89,51,59,53,79,91)
stem(v)

## 
##   The decimal point is 1 digit(s) to the right of the |
## 
##   4 | 9
##   5 | 139
##   6 | 13
##   7 | 9
##   8 | 49
##   9 | 1

japanes <- c(84,63,61,49,89,51,59,53,79,91)
mean(japanes)

## [1] 67.9

median(japanes)

## [1] 62

as.numeric(names(table(japanes)))[which.max(table(japanes))]

## [1] 49

sd(japanes)

## [1] 16.25115

var(japanes)

## [1] 264.1

Q1 <- quantile(japanes, 1 / 4) 
Q2 <- quantile(japanes, 2 / 4) 
Q3 <- quantile(japanes, 3 / 4) 

P1 <- quantile(japanes, 1 / 10) 
P2 <- quantile(japanes, 2 / 10) 
P3 <- quantile(japanes, 3 / 10) 

P5 <- quantile(japanes, 5 / 10) 

Q1

##  25% 
## 54.5

Q2

## 50% 
##  62

Q3

##   75% 
## 82.75