R/maize_multible_loc_analysis.R

In dengfei2013/met:

#' Return anova table to a norm way
#'
#' @description
#' \code{maize_met_analysis} Return anova and LSD for the MET analysis
#' @param MET data: data: Loc, Rep, Cul, yield
#'
#' @examples
#' library(met)
#' data("maize")
#' head(maize)
#' dat = maize
#' re = maize_met_analysis(dat,json=F)
#' re
#' maize_met_analysis(dat)
#' data(hch)
#' head(hch)
#' maize_met_analysis(hch)
#' maize_met_analysis(hch,json=F)
#' tt = maize_met_analysis(hch,json=F)



maize_multible_loc_analysis = function(dat,json=TRUE){
  suppressMessages(require(jsonlite))
  names(dat) = c("Loc","Rep","Cul","yield")
  dat$Loc.Rep = as.factor(paste0(dat$Loc,dat$Rep))
  dat$Loc.Cul = as.factor(paste0(dat$Loc,dat$Cul))
  mod = aov(yield ~ Loc + Loc.Rep + Cul + Loc.Cul, data = dat)
  aa = anova_tab(mod)

  nn = length(aa[[1]][,3])
  mse = aa[[1]][,3][(nn-1)]
  mean = mean(dat$yield,na.rm=TRUE)
  cv = round(sqrt(mse)*100/mean,4)

  cul_lsd0.05 = LSD_test(mod,"Cul",alpha = 0.05)$groups
  cul_lsd0.05_Name = rownames(cul_lsd0.05)
  cul_lsd0.05$Name =  rownames(cul_lsd0.05)
  cul_lsd0.05 = cul_lsd0.05[,c(3,1,2)]

  cul_lsd0.01 = LSD_test(mod,"Cul",alpha = 0.01)$groups
  cul_lsd0.01_Name = rownames(cul_lsd0.01)
  cul_lsd0.01$Name =  rownames(cul_lsd0.01)
  cul_lsd0.01 = cul_lsd0.01[,c(3,1,2)]

  cul_lsd_value0.05 = LSD_value(mod,"Cul",alpha = 0.05);cul_lsd_value0.01 = LSD_value(mod,"Cul",alpha = 0.01)

  loc_lsd0.05 = LSD_test(mod,"Loc",alpha=0.05)$groups
  loc_lsd0.05_Name = rownames(loc_lsd0.05)
  loc_lsd0.05$Name = rownames(loc_lsd0.05)
  loc_lsd0.05 = loc_lsd0.05[,c(3,1,2)]

  loc_lsd0.01 = LSD_test(mod,"Loc",alpha=0.01)$groups
  loc_lsd0.01_Name = rownames(loc_lsd0.01)
  loc_lsd0.01$Name = rownames(loc_lsd0.01)
  loc_lsd0.01 = loc_lsd0.01[,c(3,1,2)]

  loc_lsd_value0.05 = LSD_value(mod,"Loc",alpha = 0.05);loc_lsd_value0.01 = LSD_value(mod,"Loc",alpha = 0.01)

  loc.cul_lsd0.05 = LSD_test(mod,"Loc.Cul",alpha = 0.05)$groups
  loc.cul_lsd0.05_Name = rownames(loc.cul_lsd0.05)
  loc.cul_lsd0.05$Name = rownames(loc.cul_lsd0.05)
  loc.cul_lsd0.05 = loc.cul_lsd0.05[,c(3,1,2)]

  loc.cul_lsd0.01 = LSD_test(mod,"Loc.Cul",alpha = 0.01)$groups
  loc.cul_lsd0.01_Name = rownames(loc.cul_lsd0.01)
  loc.cul_lsd0.01$Name = rownames(loc.cul_lsd0.01)
  loc.cul_lsd0.01 = loc.cul_lsd0.01[,c(3,1,2)]

  loc.cul_lsd_value0.05 = LSD_value(mod,"Loc.Cul",alpha = 0.05);loc.cul_lsd_value0.01 = LSD_value(mod,"Loc.Cul",alpha = 0.01)

  re = c(aa,cul_lsd0.05,cul_lsd0.01,cul_lsd_value0.05,cul_lsd_value0.01,
         loc_lsd0.05,loc_lsd0.01,loc_lsd_value0.05,loc_lsd_value0.01,
         loc.cul_lsd0.05,loc.cul_lsd0.01,loc.cul_lsd_value0.05,loc.cul_lsd_value0.01,
         cv)

  if(json){
    re = lapply(re, toJSON)
  }else{
    re = list(aa,cul_lsd0.05,cul_lsd0.01,cul_lsd_value0.05,cul_lsd_value0.01,
           loc_lsd0.05,loc_lsd0.01,loc_lsd_value0.05,loc_lsd_value0.01,
           loc.cul_lsd0.05,loc.cul_lsd0.01,loc.cul_lsd_value0.05,loc.cul_lsd_value0.01,
           cv)
  }
  return(re)
}