R/RMttest.R

In SSDforR: Functions to Analyze Single System Data

RMttest <-
function(behavior,phaseX, v1,v2){
 options (scipen=999)
  options(warn=-1)
  
  t1<-table(phaseX)
  tmaxA<-t1[names(t1)==v1]
  startA<-match(v1,phaseX)
  endA<-tmaxA+startA-1
  A<-behavior[startA:endA]
  meanA=mean(A,na.rm=T)
 # x1=rep(1:tmaxA, each = v1)
  
  
  tmaxB<-t1[names(t1)==v2]
  startB<-match(v2,phaseX)
  endB<-tmaxB+startB-1

  B=(behavior[startB:endB])
  meanB=mean(B,na.rm=T)

t1<-t.test(A,B)
var1<-var.test(A,B)
vartestp<-var.test(A,B)$p.value
t2<-t.test(A,B,var.equal=TRUE) 
DV<-( paste(substitute(behavior)) )

writeLines(" ")
l1<-c("t-test analysis for behavior",'"',DV,'"')
cat(sprintf(l1),"\n")
writeLines(" ")
  print(t2)
  print(var1)
  print(t1)
 Means<-c(meanA,meanB) 
  
  
  
 #layout(rbind(1,2), heights=c(4,1))
 layout(rbind(1,2), heights=c(6,1))
 yt<-c("mean",DV)
 
 
  barplot(Means,ylab=yt,names.arg=c(v1,v2),xlab="phase",main="Difference of Means")
 
  if (vartestp>.05) { 
  note1<-c("Note: Your variances are equal (F test to compare two variances p-value is greater than 0.05) between phases",v1, "and", v2,".","Use the Two Sample t-test")
  cat(sprintf(note1),"\n")
  }
  
  if (vartestp<=.05) { 
    note2<-c("Note: Your variances are unequal (F test to compare two variances p-value is less than or equal to 0.05) between phases ",v1, "and", v2,".","Use the Welch Two Sample t-test")
    cat(sprintf(note2),"\n")
  }
 
  
  options(warn=-0)  }