setwd("mt")
library("Hmisc") 

datamt <- read.csv('spring11.csv',header=T)
datamt[is.na(datamt)] <- 0

#######################################
#### MMI
#######################################

mt.score <- function(task,columns,data)
{

	data$total<-0
	for(i in 1:length(columns))
	{
		data[,columns[i]]<-data[,columns[i]]/100
		data$total<-data$total+(data[,columns[i]]*data[,task])
	}
	return(data$total)
}

datamt$b1<-mt.score("Q1.1",c("Q1.2_1","Q1.2_2","Q1.2_3","Q1.2_4","Q1.2_5","Q1.2_6","Q1.2_7","Q1.2_8","Q1.2_9","Q1.2_10"),datamt)
datamt$b2<-mt.score("Q2.1",c("Q2.2_2","Q2.2_3","Q2.2_4","Q2.2_5","Q2.2_6","Q2.2_7","Q2.2_8","Q2.2_9","Q2.2_10"),datamt)
datamt$b3<-mt.score("Q3.1",c("Q3.2_2","Q3.2_3","Q3.2_4","Q3.2_5","Q3.2_6","Q3.2_7","Q3.2_8","Q3.2_9","Q3.2_10"),datamt)
datamt$b4<-mt.score("Q4.1",c("Q4.2_2","Q4.2_3","Q4.2_4","Q4.2_5","Q4.2_6","Q4.2_7","Q4.2_8","Q4.2_9","Q4.2_10"),datamt)
datamt$b5<-mt.score("Q5.1",c("Q5.2_2","Q5.2_3","Q5.2_4","Q5.2_5","Q5.2_6","Q5.2_7","Q5.2_8","Q5.2_9","Q5.2_10"),datamt)
datamt$b6<-mt.score("Q6.1",c("Q6.2_2","Q6.2_3","Q6.2_4","Q6.2_5","Q6.2_6","Q6.2_7","Q6.2_8","Q6.2_9","Q6.2_10"),datamt)
datamt$b7<-mt.score("Q7.1",c("Q7.2_2","Q7.2_3","Q7.2_4","Q7.2_5","Q7.2_6","Q7.2_7","Q7.2_8","Q7.2_9","Q7.2_10"),datamt)
datamt$b8<-mt.score("Q8.1",c("Q8.2_2","Q8.2_3","Q8.2_4","Q8.2_5","Q8.2_6","Q8.2_7","Q8.2_8","Q8.2_9","Q8.2_10"),datamt)
datamt$b9<-mt.score("Q9.1",c("Q9.2_2","Q9.2_3","Q9.2_4","Q9.2_5","Q9.2_6","Q9.2_7","Q9.2_8","Q9.2_9","Q9.2_10"),datamt)
datamt$b10<-mt.score("Q10.1",c("Q10.2_2","Q10.2_3","Q10.2_4","Q10.2_5","Q10.2_6","Q10.2_7","Q10.2_8","Q10.2_9","Q10.2_10"),datamt)


#### Totals
#face-to-face
datamt$ffTotal<-datamt$Q1.1
#homework
datamt$homeworkTotal <-+datamt$Q8.1
#media use
datamt$mediaUseTotal<-datamt$Q2.1+datamt$Q3.1+datamt$Q4.1+datamt$Q5.1+datamt$Q6.1+datamt$Q7.1+datamt$Q8.1+datamt$Q9.1+datamt$Q10.1


#### Simultanious task totals
datamt$ffSimul<-datamt$b1
datamt$homeworkSimul<-datamt$b8
datamt$mediaUseSimul<-datamt$b2+datamt$b3+datamt$b4+datamt$b5+datamt$b6+datamt$b7+datamt$b9+datamt$b10

#### MT ratios
datamt$mtFF<- datamt$b1/datamt$ffTotal
datamt$mtHomework<-datamt$b8/datamt$homeworkTotal
datamt$mt<-datamt$mediaUseSimul/datamt$mediaUseTotal




datamt$mtthirds <- cut2(datamt$mt, g=3)

levels(datamt$mtthirds) <- c("low", "med", "high")
#low