BIODEX_2-6,25-29peak_detect.Rmd

---
title: "BIODEX_data_sampling"
author: "Yagi_Masato"
date: "2020/1/15"
output: html_document
---

```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE)
```

```{r load}
library(xlsx)
library(dplyr)
library(rlist)
fun_load <- function(files_name){
  all_data <- read.table(files_name, skip=5, header=T)
  torque <<- all_data$N.M 
  x <- rep(0, times=40)
  torque <<- append(torque,x)
  torque <<- append(torque,x,after = 0) 
  # torque <<- append(torque,x,after = 0)
  # fun_detect(torque)
}
```

```{r detect}
fun_detect <- function(torque){
  # 波形データのトップを探す
  num <- length(torque)
  i <- 41
  while (i < num) {
    # 前後40個の値より大きいだとピーク値が引っかからない
    biggest_previos <- all(torque[(i-1):(i-40)] <= torque[i])
    biggest_post <- all(torque[(i+1):(i+40)] <= torque[i])
    if (torque[i] > 30 && biggest_previos && biggest_post){
      peak_torque <<- append(peak_torque, torque[i])
      # ピーク値が見つかれば100個先に移動する
      i <- i + 100
    } else {
      # 見つからなければ１個先のデータへ
        i <- i + 1
    }
  }
}

```

``` {r print info}
fun_info <- function(){
  print(file_name)
  print(torque_length)
  print(peak_torque)
}
```

```{r first and last}
#出力するデータフレーム
# output_data <- c()
# df <- data.frame(matrix(rep(1,11),nrow = 1))#[character(0),] # 空のデータフレーム作成
# data.frame( matrix( vector(), length(files), , dimnames = list( c(), colnames ) ) )
df <- data.frame(FILE=character(),
                 "2"=character(),
                 "3"=character(),
                 "4"=character(),
                 "5"=character(),
                 "6"=character(),
                 "25"=character(),
                 "26"=character(),
                 "27"=character(),
                 "28"=character(),
                 "29"=character(),
                 stringsAsFactors=FALSE)

fun_CreateDf <- function(file_name,peak_torque){
  torque_length <- length(peak_torque)
  if (torque_length == 30) {
    # w <- list(file_name,peak_torque[1:5],peak_torque[26:30]) %>% list.ungroup() 型を保持したいもののunlistするとchrになってしまう
    # dfにトルクのベクトルを追加
    df[(nrow(df) + 1),] <<- c(file_name,peak_torque[2:6],peak_torque[25:29])
    print("↓波形データを検出しました！！")
  } else if (torque_length == 31 && peak_torque[30] / 2 > peak_torque[31]){
    df[(nrow(df) + 1),] <<- c(file_name,peak_torque[2:6],peak_torque[25:29])
    print("↓一回多く検出されました")
  } else if (torque_length == 29){
    print("↓波形データが１つとんでいます！！")
    df[(nrow(df) + 1),] <<- c(file_name,numeric(10))
  } else {
    print("##################↓正しく波形データが検出できませんでした！！################")
    df[(nrow(df) + 1),] <<- c(file_name,numeric(10))
  }
  print(file_name)
  print(torque_length)
  print(peak_torque)
  peak_torque <<- c()
}
  
```

```{r ref}
setwd("C:/Users/yagin/Documents/BIODEX")
#上記のディレクトリ内の全てのCSVファイルを列挙
#list all csv files in above directory
files <- list.files(pattern = "\\.txt$")
peak_torque <- c()
for (f in 1:length(files)){ #fileの数だけ実行
  print(files[f]) #ファイル名を出力
  fun_load(files[f]) #各ファイルからトルクを取り出す
  fun_detect(torque) #波形データの取り出し
  #ファイル名とピーク値を出力
  fun_CreateDf(files[f],peak_torque)
}
write.xlsx(df, "C:/Users/yagin/Desktop/output.xlsx", sheetName="Sheet1",col.names=T,row.names=F)

```

```{r EXCEL Connector}
# setwd('C:/Users/yagin/Desktop')
# library("XLConnect")
# wb <- loadWorkbook("XLConnect.xlsx", create = TRUE)
# createSheet(wb, name = "mtcars")
# mtcars <- data.frame(Group = I(paste0("TESTてすと", 1:20)),Data = 1:20)
# writeWorksheet(wb, mtcars, sheet = "mtcars", startRow = 5, startCol = 5, header = FALSE)
# saveWorkbook(wb)
```





```{r excel}



```