#四則演算など

1+2

4-1

3*2

10/5

2^3

#変数への値の代入

a <- 2

a

a <- a + 1

a

#ベクトルの入力

shopa <- c(75,82,63,52,89,73,64,75,85,70,49,85,69,68,70)

shopa

#関数の利用

sqrt(2)

a <- 2

sqrt(a)

#最大値，最少値，合計，平均値，中央値，分散，標準偏差

max(shopa)

min(shopa)

sum(shopa)

mean(shopa)

median(shopa)

var(shopa)

sd(shopa)

#共分散と相関

shopb <- c(77,87,59,56,52,56,76,90,87,79,81,70,76,70,81)

var(shopa, shopb)

cor(shopa , shopb)

#作業用フォルダの確認と設定．設定する場合には下の＃を取り除き，setwd("フォルダの位置")とする．

getwd()

# setwd()

#外部データの読み込み

v08c0.dat <- read.csv("v08c0.csv", fileEncoding="Shift_JIS")

v08c0.dat

#特定のデータの取り出し

v08c0.dat[8,2]

v08c0.dat[1,]

v08c0.dat[,1]

v08c0.dat[,]

#リストの計算

sum(v08c0.dat)

sum(v08c0.dat[,1])

mean(v08c0.dat[,1])

median(v08c0.dat[,1])

var(v08c0.dat[,1])

sd(v08c0.dat[,1])

var(v08c0.dat)

cor(v08c0.dat)

#散布図

plot(v08c0.dat)

#折れ線

plot(v08c0.dat[,1], type = "l")

#Macで日本語が表示できない場合，以下の行を＃を外取り除いて実行．

# par(family = "HiraKakuProN-W3")

#パッケージの読み込み

library("MASS")

#ヘルプの閲覧

?survey

#調査データの閲覧と平均値算出

head(survey)

#平均値を求める(方法1)

mean(survey[,10], na.rm = TRUE)

＃平均値を求める(方法2)

mean(survey$Height, na.rm = TRUE)

#ヒストグラムの作成

hist(survey$Height)

hist(survey$Height, breaks = c(150,160,170,180,190,200))

#クロス集計

table(survey$Sex, survey$Smoke)

#====================================================

#2022.09.08 齋藤朗宏・荘島宏二郎