Nothing
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)
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.