R/F5.R

In GlenDu/ODBbB: Operating on Data block by block

#' SD of transposed Data block by block
#'@description Takes any data fame and transpose the blocks predeterminated in the function, than it give
#' the standard deviation of the column or the standard deviation of the the rows of each block.
#' It also gives the total standard deviation of each transposed block. ' the hole block transposed.
#'
#' @param Data Data frame that will be imported.
#' @param lengthBlock Nomber of rows that will be transposed for any row.
#' @param cls Columns of the  data frame that we want to transpose.
#' @param rs The rows of the first block. Based on this rows the function will select identically the rows of every other block.
#' @param time_col The 'time parameter' is optional. We can ignore it if our data has no time values.
#' @param SD Character that can take the values Row, Col or Block. Based on the option.
#'that we'll choose and the function will return the mean of the Row, Col of each block or the general mean of all the block.
#'
#' @return tibble
#' @export
#'
BlockTrSD<-function(Data=data(),lengthBlock=integer(),cls=c(),rs=c(),time_col=NULL,SD=character()){
  dist_val=lengthBlock-1
  endT<-nrow(Data)-lengthBlock

  Tdat=NULL
  a1=b=c=rw=bl=cl =NULL
  time<-Data[,time_col]
  time<-(unique(time))
  for(j in seq(from=0, to=endT,by=lengthBlock)){
    rowSel<-rs+j
    a<-t(Data[rowSel,cls])
    if(SD=="Row"){
      rw<-apply(a,1,sd)
    }
    if(SD=="Col"){
      cl<-apply(a,2,sd)
    }
    if(SD=="Block"){
      bl<-sd(a)
    }

    a1<-rbind(a1,rw)
    b<-rbind(b,cl)
    c<-rbind(c,bl)
  }
  if (SD==" "){
    stop("Specify Row, Col or Block")}
  if((SD=="Row") & (is.null(time_col))){return(a1)}
  else if((SD=="Row") & (!is.null(time_col))){  a1t<-cbind(time,a1)
  return(a1t)}
  if((SD=="Col") & (is.null(time_col))){ return(b)}
  else if((SD=="Col") & (!is.null(time_col))){  bt<-cbind(time,b)
  return(bt)}
  if((SD=="Block") & (is.null(time_col))){ return(c)}
  else if((SD=="Block") & (!is.null(time_col))){  ct<-cbind(time,c)
  return(ct)}


}