R/mulamba_mock_funcition.R

In AgronomiaR/IS:

#' Indice de selecao de Mulamba e Mock simplificado
#'
#' @param IS objeto de classe mulamba
#' @param dados data.frame com medias ou saida da PCA
#' @param ps porcentagem de selecao (default e 15)
#' @export
#' @example
#'
#' resp1=c(10,11,12,7,6,5,4,4,5,12,7,16,5,4,4,5)
#' resp2=c(10,11,12,7,6,5,4,5,5,12,7,6,15,4,4,5)
#' resp3=c(10,11,12,7,6,5,4,4,5,12,7,6,10,4,7,5)
#' resp4=c(6,1,2,7,6,3,4,3,5,12,7,6,5,14,4,15)
#' trat=rep(paste("T",1:16))
#' dados=data.frame(trat,resp1,resp2,resp3,resp4)
#' library(IS)
#' mulamba_mock(dados,orde = c("1","1","1","d"),ps=30,"no")
#' mulamba_mock(dados,orde = c("1","1","1","d"),ps=30,"yes")

mulamba_mock=function(dados,orde=NULL,ps=15,variable="no"){
  a=mulamba(dados,orde=orde)
  if(variable=="no"){plot_mulamba(a,dados,ps=ps)}
  if(variable=="yes"){plot_mulamba2(a,dados,ps=ps)}
}