prep_data.R

library(psych)
library(plyr)
setwd("~/Dropbox/r-code/SUS/")

extractQ = function(x){
    temp = as.factor(as.character(x[-1]))
    return(temp)
}
extractQ2 = function(x){
    temp = as.character(x[-1])
    return(temp)
}
fac2num = function(x){
    return(as.numeric(as.character(x)))
}

############################
###### Med Stats Data  #####
############################

dat=read.csv("raw_data/medstats_250813.csv")
dat=read.csv("raw_data/medstats_280813.csv")
dat=read.csv("raw_data/medstats_06092013.csv")
names(dat)


# Remove duplicates and test responses
dat = subset(dat, a!="")
dat = dat[!duplicated(dat$a),]

# Correlates
oe=extractQ(dat$Q1.1)
oe=fac2num(factor(oe, labels=c(NA,0,6,2,4,3,1,5)))

difficulty=extractQ(dat$Q1.2)
difficulty=fac2num(factor(difficulty, labels=c(NA,0,6,3,2,4,1,5)))

redo=extractQ(dat$Q1.3)
redo=6-fac2num(factor(redo, labels=c(NA,0,6,2,4,3,1,5)))

learnsat=extractQ(dat$Q1.4)
learnsat=fac2num(factor(learnsat, labels=c(NA,0,6,3,2,4,1,5)))

cor(cbind(oe,difficulty,redo,learnsat), use="pair")

browser=cbind(extractQ2(dat$Q59_1_TEXT), extractQ2(dat$Q76_1_TEXT), extractQ2(dat$Q62_1_TEXT))
os=cbind(extractQ2(dat$Q59_3_TEXT), extractQ2(dat$Q76_3_TEXT), extractQ2(dat$Q62_3_TEXT))
res=cbind(extractQ2(dat$Q59_4_TEXT), extractQ2(dat$Q76_4_TEXT), extractQ2(dat$Q62_4_TEXT))
for(i in 1:nrow(browser)){
	browser[i,1] = ifelse( sum(browser[i,]!="")==1 , browser[i,(browser[i,]!="")], "")
	os[i,1] = ifelse( sum(os[i,]!="")==1 , os[i,(os[i,]!="")], "")
	res[i,1] = ifelse( sum(res[i,]!="")==1 , res[i,(res[i,]!="")], "")
}
browser=browser[,1]
os=os[,1]
res=res[,1]

# try to get timing which didn't export right

# How easy or difficult was it to learn how to use the system?
q1=extractQ(dat$Q45)
q1=6-fac2num(factor(q1, labels=c(NA,6,0,5,1,3,4,2)))
# How much more is there to learn about the system?
q2=extractQ(dat$Q46)
q2=4-fac2num(factor(q2, labels=c(NA,4,1,3,2,0)))
# How confident were you using the system?
q3=extractQ(dat$Q47)
q3=fac2num(factor(q3, labels=c(NA,4,2,0,1,3)))
# How easy or difficult is it to use the system?
q4=extractQ(dat$Q50)
q4=6-fac2num(factor(q4, labels=c(NA,6,0,5,1,3,4,2)))
# How consistent is the system?
q5=extractQ(dat$Q51)
q5=fac2num(factor(q5, labels=c(NA,4,2,0,1,3)))

sus2=data.frame(q1,q2,q3,q4,q5)

sus2_all=data.frame(sus2, browser, os, res, oe, difficulty, redo, learnsat)[rowSums(is.na(sus2))==0,]
write.csv(sus2_all, "sus2_medstat.csv", row.names=F)



# REV SUS
# I needed to learn very few things before I could get going with this system
# I think that I would not need any support of a technical person to be able to use this system
# I did not feel very confident using the system
# I found the system very manageable to use
# I would imagine that most people would learn to use this system very slowly
# I found the system appropriately simple
# I thought the system was hard to use
# I found the various functions in this system were not well integrated
# I did not think there was too much inconsistency in this system
# I do not think that I would like to use this system frequently
labAD = function(x, rev=0){
	x=as.character(x[-1])
	x[x=="Strongly disagree (1)"]=0
	x[x=="(2)"]=1
	x[x=="(3)"]=2
	x[x=="(4)"]=3
	x[x=="Strongly agree (5)"]=4
	x=abs((-4*rev)+as.numeric(x))
	return(x)
}
revsus=subset(dat, select=Q75_2:Q75_15)
revsus[1,]
rev_sus=data.frame(
	q1=labAD(revsus[,10],1),
	q2=labAD(revsus[,6]),
	q3=labAD(revsus[,7],1),
	q4=labAD(revsus[,2]),
	q5=labAD(revsus[,8],1),
	q6=labAD(revsus[,9]),
	q7=labAD(revsus[,5],1),
	q8=labAD(revsus[,4]),
	q9=labAD(revsus[,3],1),
	q10=labAD(revsus[,1]))

rev_sus_all=data.frame(rev_sus, browser, os, res, oe, difficulty, redo, learnsat)[rowSums(is.na(rev_sus))==0,]
write.csv(rev_sus_all, "rev_sus_medstat.csv", row.names=F)



# Original SUS
# I needed to learn a lot of things before I could get going with this system
# I think that I would need the support of a technical person to be able to use this system
# I felt very confident using the system
# I found the system very cumbersome to use
# I would imagine that most people would learn to use this system very quickly
# I found the system unnecessarily complex
# I thought the system was easy to use
# I found the various functions in this system were well integrated
# I thought there was too much inconsistency in this system
# I think that I would like to use this system frequently

orisus=subset(dat, select=Q59_2:Q59_15)
orisus[1,]
ori_sus=data.frame(
	q1=labAD(orisus[,10]),
	q2=labAD(orisus[,6],1),
	q3=labAD(orisus[,7]),
	q4=labAD(orisus[,2],1),
	q5=labAD(orisus[,8]),
	q6=labAD(orisus[,9],1),
	q7=labAD(orisus[,5]),
	q8=labAD(orisus[,4],1),
	q9=labAD(orisus[,3]),
	q10=labAD(orisus[,1],1))

ori_sus_all=data.frame(ori_sus, browser, os, res, oe, difficulty, redo, learnsat)[rowSums(is.na(ori_sus))==0,]
write.csv(ori_sus_all, "ori_sus_medstat.csv", row.names=F)


############################
###### Algorithms Data #####
############################

dat=read.csv("raw_data/algo_06092013.csv")
dat=read.csv("raw_data/algo_08092013.csv")
names(dat)

# Correlates

#How good or bad was your overall experience with the course?
dat$Q1.1 #extremely bad=1, extremely good=7
oe=fac2num(dat$Q1.1[-1])-1

#How easy or difficulty was the course?
dat$Q1.2 #extremely difficult=1, extremely easy=7
difficulty=fac2num(dat$Q1.2[-1])-1

#How likely or unlikely are you take another course with the same format?
dat$Q1.3 #extremely unlikely=1, extremely likely=7
redo=fac2num(dat$Q1.3[-1])-1

#How satisfied or dissatisfied are you with the amount you learnt in the course?
dat$Q1.4  #1=Extremely Dissatisfied,2=Very Dissatisfied,3=Somewhat Dissatisfied,4=Neutral,5=Somewhat Satisfied,8=Very Satisfied,9=Extremely Satisfied
learnsat=fac2num(dat$Q1.4[-1])-1
learnsat[!is.na(learnsat) & learnsat>6]=learnsat[!is.na(learnsat) & learnsat>6]-2

cor(cbind(oe,difficulty,redo,learnsat), use="pair")

# browser=cbind(extractQ2(dat$Q59_1_TEXT), extractQ2(dat$Q76_1_TEXT), extractQ2(dat$Q62_1_TEXT))
# os=cbind(extractQ2(dat$Q59_3_TEXT), extractQ2(dat$Q76_3_TEXT), extractQ2(dat$Q62_3_TEXT))
# res=cbind(extractQ2(dat$Q59_4_TEXT), extractQ2(dat$Q76_4_TEXT), extractQ2(dat$Q62_4_TEXT))
# for(i in 1:nrow(browser)){
#     browser[i,1] = ifelse( sum(browser[i,]!="")==1 , browser[i,(browser[i,]!="")], "")
#     os[i,1] = ifelse( sum(os[i,]!="")==1 , os[i,(os[i,]!="")], "")
#     res[i,1] = ifelse( sum(res[i,]!="")==1 , res[i,(res[i,]!="")], "")
# }
# browser=browser[,1]
# os=os[,1]
# res=res[,1]
names(dat)

#How easy or difficult was it to learn how to use the system?
dat$Q54
# The scale for this question was displayed incorrectly! Cannot use it!

#How much more is there to learn about the system?
dat$Q56 #0=nothing at all, 4=a great deal
q2=4-fac2num(dat$Q56[-1])

#How confident were you using the system?
dat$Q48 # 0=not at all confident, 4=extremely confident
q3=fac2num(dat$Q48[-1])

#How easy or difficult is it to use the system?
dat$Q52 #0=extremely diff, 6=extremely easy
q4=fac2num(dat$Q52[-1])

#How consistent is the system?
dat$Q50 # 0=not at all consistent, 4=extremely consistent
q5=fac2num(dat$Q50[-1])

#How likely are you to need support of a technical person to be able to use the system?
dat$Q60 #0=extremely unlikely, 6=extremely likely
q6=6-fac2num(dat$Q60[-1])

#How much did you like or dislike the system?
dat$Q63 #0=extremely dislike, 6=extremely like
q7=fac2num(dat$Q63[-1])

#How complex is the system?
dat$Q64 #0=not at all, 4=extremely complex
q8=4-fac2num(dat$Q64[-1])

#How cumbersome is it to use the system?
dat$Q65 #0=not at all cumbersome, 4=extremely cumbersome
q9=4-fac2num(dat$Q65[-1])

#How integrated were the system’s various functions?
dat$Q61 #0=not at all integrated, 4=extremely integrated
q10=fac2num(dat$Q61[-1])

sus2=data.frame(q1=NA,q2=q2,q3=q3,q4=q4,q5=q5,q6=q6,q7=q7,q8=q8,q9=q9,q10=q10)

sus2_all=data.frame(sus2, oe, difficulty, redo, learnsat)[rowSums(is.na(sus2))==1,]
write.csv(sus2_all, "sus2_algo.csv", row.names=F)





# REV SUS
# I needed to learn very few things before I could get going with this system
# I think that I would not need any support of a technical person to be able to use this system
# I did not feel very confident using the system
# I found the system very manageable to use
# I would imagine that most people would learn to use this system very slowly
# I found the system appropriately simple
# I thought the system was hard to use
# I found the various functions in this system were not well integrated
# I did not think there was too much inconsistency in this system
# I do not think that I would like to use this system frequently
rev_sus=data.frame(
    q1=5-fac2num(dat$Q63_15[-1]),
    q2=fac2num(dat$Q63_11[-1])-1,
    q3=5-fac2num(dat$Q63_12[-1]),
    q4=fac2num(dat$Q63_3[-1])-1,
    q5=5-fac2num(dat$Q63_13[-1]),
    q6=fac2num(dat$Q63_14[-1])-1,
    q7=5-fac2num(dat$Q63_6[-1]),
    q8=fac2num(dat$Q63_5[-1])-1,
    q9=5-fac2num(dat$Q63_4[-1]),
    q10=fac2num(dat$Q63_2[-1])-1
)

rev_sus_all=data.frame(rev_sus, oe, difficulty, redo, learnsat)[rowSums(is.na(rev_sus))==0,]
write.csv(rev_sus_all, "rev_sus_algo.csv", row.names=F)



# Original SUS
# I needed to learn a lot of things before I could get going with this system
# I think that I would need the support of a technical person to be able to use this system
# I felt very confident using the system
# I found the system very cumbersome to use
# I would imagine that most people would learn to use this system very quickly
# I found the system unnecessarily complex
# I thought the system was easy to use
# I found the various functions in this system were well integrated
# I thought there was too much inconsistency in this system
# I think that I would like to use this system frequently
orig_sus=data.frame(
    q10=fac2num(dat$Q56_15[-1])-1,
    q6=5-fac2num(dat$Q56_11[-1]),
    q7=fac2num(dat$Q56_12[-1])-1,
    q2=5-fac2num(dat$Q56_3[-1]),
    q8=fac2num(dat$Q56_13[-1])-1,
    q9=5-fac2num(dat$Q56_14[-1]),
    q5=fac2num(dat$Q56_6[-1])-1,
    q4=5-fac2num(dat$Q56_5[-1]),
    q3=fac2num(dat$Q56_4[-1])-1,
    q1=5-fac2num(dat$Q56_2[-1])
)

orig_sus_all=data.frame(orig_sus, oe, difficulty, redo, learnsat)[rowSums(is.na(orig_sus))==0,]
write.csv(orig_sus_all, "ori_sus_algo.csv", row.names=F)