d <- read.csv(file = ‘https://stats.dip.jp/01_ds/data/dirty_data2.csv’)

datatable(d)

d1 <- na.omit(d)

summary(d1) # 要約統計値　NAがなくなった。

各気象変数の許容範囲に入っているか否かを示す論理式を作成

is.tp <- -20 <= d1\(気温 & d1\)気温 <= 50 is.hm <- 0 <= d1\(湿度 & d1\)湿度 <= 100 is.ws <- 0 <= d1\(風速 & d1\)風速 <= 70 is.pp <- -1 <= d1\(全天日射量MJ&d1\)全天日射量MJ <= 4 summary(is.tp) summary(is.hm) summary(is.ws) summary(is.pp) # 適切な範囲のレコードだけを抽出する。 d2 <- d1[is.tp & is.hm & is.ws & is.pp, ]

datatable(d2)

d3 <- d2 d3\(気温 <- round(d2\)気温, 1) d3\(湿度 <- round(d2\)湿度, 0) d3\(風速 <- round(d2\)風速, 0) d3\(全天日射量MJ <- round(d2\)全天日射量MJ, 1) datatable(d3)

d3\(会社名 LAB.TOSHIBA <- c("東芝" ,"Toshiba","TOSHIBA" ,"toshiba" ,"東芝株式会社", "東芝（株）" ,"東芝(株)" ) LAB.TOYOTA <- c("トヨタ","Toyota" ,"トヨタ自動車株式会社","トヨタ(株)","toyota" ) is.toshiba <- d3\)会社名 %in% LAB.TOSHIBA is.toyota <- d3$会社名 %in% LAB.TOYOTA

d4 <- d3

d4\(会社名[is.toshiba] <- '東芝株式会社' d4\)会社名[is.toyota] <- ‘トヨタ自動車株式会社’

d4$会社名

datatable(d4)