## Robust mediaiton analysis

## data
dset<-read.table('active.txt', na.string='99999')

sel<-((is.na(dset[,2])) + (is.na(dset[,3])))

comp<-dset[sel==0,]

comp300<-comp[1:300,]

library(moments)

apply(comp300, 2, kurtosis)