R/faS.R

In sontron/madis: Manipulating and Analyzing Data Interactively with Shiny

#' faS
#' 
#' factor analysis used in shiny app.
#' 
#' @export


faS<-function(
  data,
  vars,
  nfcts=2,
  Rotate=c('none','varimax','quartimax','bentlerT',
           'varmin','equamax','geominT','bifactor','promax','oblimin',
           'simplimax','bentlerQ','geominQ','biquartimin','cluster')[1],
  Scores=c('regression','Thurstone','tenBerge','Anderson','Barlett')[1],
  FM=c('minres','uls','ols','wls','gls','pa','ml','minchi','minrank')[1],
  subset='all',
  faVarName='',
  addVar=T
){
  require('psych')
  require('stringi')
  
  if(is.character(data)) data=eval(as.name(data))
  if(subset=='all'){
    data<-data
  } else {
    subset(data,eval(parse(text=subset)))->data
  }
  
  # data[,vars]->dat
  # na.omit(dat)->dat
  unlist(stri_split_fixed(vars,','))->vars
  data[,vars]->dt
  which(sapply(dt,class)%in%c('numeric','integer'))->indVar
  dt[,indVar]->dt2
  fa(dt2,nfactors = nfcts,rotate=Rotate,scores=Scores,fm=FM)->res
  
  if(addVar){
    if(faVarName==''){
      cbind(data,res$scores)->data
    } else {
      colnames(res$scores)<-paste(colnames(res$scores),faVarName,sep='_')
      cbind(data,res$scores)->data
    }
  } else {
    NULL
  }
  
  return(list(resFA=res,dtFA=data,dataScree=dt2,cumVar=res$Vaccounted))
  
  
}