R/ABdescrip.R

In SSDforR: Functions to Analyze Single System Data

ABdescrip <-
  function(behavior,PhaseX){
   
    
    

    writeLines("-----------n-------------")
    t1<-table(PhaseX)
    print(t1)
    abmean<-tapply(behavior,PhaseX,mean,na.rm=T)
    pmean<-c(round(abmean,3))
    writeLines("-----------mean-------------")
    print(pmean)
    
    tabmean<-tapply(behavior,PhaseX,mean,trim=.1,na.rm=T)
    tmean<-c(round(tabmean,3))
    writeLines("-----------10% trim mean-------------")
    print(tmean)
    
    abmedian<-tapply(behavior,PhaseX,median,na.rm=T)
    pmedian<-c(round(abmedian,3))
    writeLines("----------median------------")
    print(pmedian)
    absd<-tapply(behavior,PhaseX,sd,na.rm=T)
    psd<-c(round(absd,3))
    writeLines("------------SD--------------")
    print(psd)
    cv<-psd/pmean
    pcv<-c(round(cv,3))
    writeLines("------------CV--------------")
    print(pcv)
    
    writeLines("---------range----------")
    abrange<-tapply(behavior,PhaseX,range,na.rm=T)
    prange<- do.call("rbind", abrange)
    print(prange)
    
    
    writeLines("---------iqr----------")
    abiqr<-tapply(behavior,PhaseX,IQR,na.rm=T)
    print(abiqr)
    
    writeLines("---------quantiles----------")
    abquan<-(tapply(behavior,PhaseX,quantile,na.rm=T))
      
     print(abquan,quote=FALSE)
    
    #do.call("rbind", abquan)
    
    layout(rbind(1,2), heights=c(4,1))
    boxplot(behavior~PhaseX)
    
  }