UQM paper - marking times
KZ
11 July 2015
Import, clean, transform, export
BIOL1040 and BIOM2011, 2013, Sem 1 and 2
#clean workspace
rm(list=ls())
#packages
require(reshape2)## Loading required package: reshape2require(dplyr)## Loading required package: dplyr
##
## Attaching package: 'dplyr'
##
## The following objects are masked from 'package:stats':
##
## filter, lag
##
## The following objects are masked from 'package:base':
##
## intersect, setdiff, setequal, union#functions
source("UQM functions.R")
source("HelperFunctions.R")
#load data
ProjectID = read.csv("./Data Stage 1/ProjectIDs.csv")
projects()
folder = "./Data Stage 2/"
SubID.index = uqm.csv("SubIDindex")
folder = "./Data Stage 1/"
mtimes = function(projects, project.names)
{
file.names <- lapply(projects, function(x){paste0(folder, "mt", x, ".csv")})
mt <- lapply(file.names, read.csv, header=TRUE, stringsAsFactors=FALSE)
names(mt) <- project.names
return(mt)
}
mt.ls = mtimes(project.num, project.names.formative) #creates a list for the 11 reports
ls = mt.ls
for (l in 1:length(ls))
print(lapply(ls, nrow)[[l]])## [1] 738
## [1] 758
## [1] 729
## [1] 729
## [1] 973
## [1] 981
## [1] 910
## [1] 257
## [1] 256
## [1] 89
## [1] 82ls = mt.ls
for (l in 1:length(ls))
print(tail(ls[[l]]))## Submission.ID Marking_Time..seconds. Text.... Drawing....
## 733 3395 2308 53 30
## 734 3406 407 0 0
## 735 3407 811 0 0
## 736 5454 250 0 0
## 737 5492 0 0 0
## 738 5493 0 0 0
## Recording....
## 733 17
## 734 100
## 735 100
## 736 0
## 737 0
## 738 0
## Submission.ID Marking_Time..seconds. Text.... Drawing....
## 753 5456 356 0 0
## 754 5473 197 50 50
## 755 5474 0 0 0
## 756 5496 0 0 0
## 757 5497 186 100 0
## 758 5499 554 0 0
## Recording....
## 753 100
## 754 0
## 755 100
## 756 0
## 757 0
## 758 0
## Submission.ID Marking_Time..seconds. Text.... Drawing....
## 724 5469 884 88 12
## 725 5470 97 0 0
## 726 5471 209 67 0
## 727 5472 698 86 14
## 728 5498 0 0 0
## 729 5502 0 0 0
## Recording....
## 724 0
## 725 0
## 726 33
## 727 0
## 728 0
## 729 0
## Submission.ID Marking_Time..seconds. Text.... Drawing....
## 724 5505 0 0 0
## 725 5506 0 0 0
## 726 5507 0 0 0
## 727 5508 0 0 0
## 728 5509 0 0 0
## 729 5510 0 0 0
## Recording....
## 724 0
## 725 0
## 726 0
## 727 0
## 728 0
## 729 0
## Submission.ID Marking_Time..seconds. Text.... Drawing....
## 968 8452 236 0 0
## 969 8453 725 25 0
## 970 8798 282 50 17
## 971 8799 325 90 0
## 972 8806 446 100 0
## 973 9157 509 100 0
## Recording....
## 968 100
## 969 75
## 970 33
## 971 10
## 972 0
## 973 0
## Submission.ID Marking_Time..seconds. Text.... Drawing....
## 976 9158 770 70 22
## 977 9364 632 81 13
## 978 9943 105 0 0
## 979 9944 23 33 0
## 980 9945 925 0 0
## 981 9946 688 50 23
## Recording....
## 976 9
## 977 6
## 978 100
## 979 67
## 980 100
## 981 28
## Submission.ID Marking_Time..seconds. Text.... Drawing....
## 905 10686 0 0 0
## 906 10687 0 0 0
## 907 10692 10 100 0
## 908 10694 1175 64 18
## 909 10695 0 0 0
## 910 10696 194 100 0
## Recording....
## 905 0
## 906 0
## 907 0
## 908 18
## 909 0
## 910 0
## Submission.ID Marking_Time..seconds. Text.... Drawing....
## 252 3877 1978 5 44
## 253 3878 1780 11 22
## 254 4061 3046 16 63
## 255 4062 1785 20 24
## 256 4419 1875 54 47
## 257 5459 1772 56 43
## Recording....
## 252 51
## 253 67
## 254 22
## 255 56
## 256 0
## 257 0
## Submission.ID Marking_Time..seconds. Text.... Drawing....
## 251 5515 2535 66 35
## 252 5516 721 80 20
## 253 5517 842 56 44
## 254 5928 2266 59 41
## 255 5929 2212 61 39
## 256 5935 1208 62 39
## Recording....
## 251 0
## 252 0
## 253 0
## 254 0
## 255 0
## 256 0
## Submission.ID Marking_Time..seconds. Text.... Drawing.... Recording....
## 84 8780 455 0 8 92
## 85 8781 261 0 0 100
## 86 8797 877 29 14 57
## 87 9152 994 0 0 100
## 88 9365 3335 11 61 28
## 89 9366 1428 29 43 29
## Submission.ID Marking_Time..seconds. Text.... Drawing.... Recording....
## 77 10690 4136 0 73 27
## 78 10691 1485 19 30 52
## 79 10693 1882 10 52 38
## 80 10707 630 0 15 85
## 81 10708 292 0 0 100
## 82 10711 590 18 54 27mt = ls.df(mt.ls, 0)
mt = rename(mt, 1:5, c("SubmissionID", "marking.sec", "Text", "Drawing", "Recording"))
mt = mt %>% mutate(marking.min = marking.sec/60)
mt[,10] = round(mt[,10],2)
head(mt)## SubmissionID marking.sec Text Drawing Recording project
## 1 1487 392 94 0 6 BIOL1040Sem1Report 0
## 2 1486 2931 42 8 50 BIOL1040Sem1Report 0
## 3 1488 1587 50 6 44 BIOL1040Sem1Report 0
## 4 1489 524 0 0 100 BIOL1040Sem1Report 0
## 5 1490 1029 0 0 100 BIOL1040Sem1Report 0
## 6 1491 4203 52 40 7 BIOL1040Sem1Report 0
## course sem report marking.min
## 1 BIOL1040 Sem1 Report 0 6.53
## 2 BIOL1040 Sem1 Report 0 48.85
## 3 BIOL1040 Sem1 Report 0 26.45
## 4 BIOL1040 Sem1 Report 0 8.73
## 5 BIOL1040 Sem1 Report 0 17.15
## 6 BIOL1040 Sem1 Report 0 70.05tapply(mt[,10], mt[,6], mean)## BIOL1040Sem1Report 0 BIOL1040Sem1Report 1 BIOL1040Sem1Report 2
## 17.04890 16.20677 13.62438
## BIOL1040Sem1Report 3 BIOL1040Sem2Report 1 BIOL1040Sem2Report 2
## 11.12753 13.89927 12.77069
## BIOL1040Sem2Report 3 BIOM2011Sem1Report 1 BIOM2011Sem1Report 2
## 10.97477 37.00658 21.87836
## BIOM2011Sem2Report 1 BIOM2011Sem2Report 2
## 26.79775 15.78561tapply(mt[,10], mt[,6], range)## $`BIOL1040Sem1Report 0`
## [1] 0.00 146.83
##
## $`BIOL1040Sem1Report 1`
## [1] -0.02 160.88
##
## $`BIOL1040Sem1Report 2`
## [1] 0.00 71.38
##
## $`BIOL1040Sem1Report 3`
## [1] 0.00 135.12
##
## $`BIOL1040Sem2Report 1`
## [1] 0.00 107.23
##
## $`BIOL1040Sem2Report 2`
## [1] 0.00 121.68
##
## $`BIOL1040Sem2Report 3`
## [1] 0.00 173.37
##
## $`BIOM2011Sem1Report 1`
## [1] 0.00 304.02
##
## $`BIOM2011Sem1Report 2`
## [1] 0.0 111.4
##
## $`BIOM2011Sem2Report 1`
## [1] 0.00 64.72
##
## $`BIOM2011Sem2Report 2`
## [1] 3.10 68.93mt.f = mt[which(mt[,1] %in% SubID.index[,1]),]
head(mt.f)## SubmissionID marking.sec Text Drawing Recording project
## 1 1487 392 94 0 6 BIOL1040Sem1Report 0
## 2 1486 2931 42 8 50 BIOL1040Sem1Report 0
## 3 1488 1587 50 6 44 BIOL1040Sem1Report 0
## 4 1489 524 0 0 100 BIOL1040Sem1Report 0
## 5 1490 1029 0 0 100 BIOL1040Sem1Report 0
## 6 1491 4203 52 40 7 BIOL1040Sem1Report 0
## course sem report marking.min
## 1 BIOL1040 Sem1 Report 0 6.53
## 2 BIOL1040 Sem1 Report 0 48.85
## 3 BIOL1040 Sem1 Report 0 26.45
## 4 BIOL1040 Sem1 Report 0 8.73
## 5 BIOL1040 Sem1 Report 0 17.15
## 6 BIOL1040 Sem1 Report 0 70.05table(mt.f$project)##
## BIOL1040Sem1Report 0 BIOL1040Sem1Report 1 BIOL1040Sem1Report 2
## 682 674 671
## BIOL1040Sem1Report 3 BIOL1040Sem2Report 1 BIOL1040Sem2Report 2
## 667 897 870
## BIOL1040Sem2Report 3 BIOM2011Sem1Report 1 BIOM2011Sem1Report 2
## 854 239 241
## BIOM2011Sem2Report 1 BIOM2011Sem2Report 2
## 85 80tapply(mt.f[,10], mt.f[,6], mean)## BIOL1040Sem1Report 0 BIOL1040Sem1Report 1 BIOL1040Sem1Report 2
## 17.38937 16.58795 14.00405
## BIOL1040Sem1Report 3 BIOL1040Sem2Report 1 BIOL1040Sem2Report 2
## 11.42303 14.08667 13.47653
## BIOL1040Sem2Report 3 BIOM2011Sem1Report 1 BIOM2011Sem1Report 2
## 10.99412 39.00205 22.24365
## BIOM2011Sem2Report 1 BIOM2011Sem2Report 2
## 27.69329 15.97213tapply(mt.f[,10], mt.f[,6], range)## $`BIOL1040Sem1Report 0`
## [1] 1.35 146.83
##
## $`BIOL1040Sem1Report 1`
## [1] 0.00 160.88
##
## $`BIOL1040Sem1Report 2`
## [1] 0.50 71.38
##
## $`BIOL1040Sem1Report 3`
## [1] 1.15 135.12
##
## $`BIOL1040Sem2Report 1`
## [1] 0.23 107.23
##
## $`BIOL1040Sem2Report 2`
## [1] 0.00 121.68
##
## $`BIOL1040Sem2Report 3`
## [1] 0.00 173.37
##
## $`BIOM2011Sem1Report 1`
## [1] 4.42 304.02
##
## $`BIOM2011Sem1Report 2`
## [1] 2.42 111.40
##
## $`BIOM2011Sem2Report 1`
## [1] 4.35 64.72
##
## $`BIOM2011Sem2Report 2`
## [1] 3.10 68.93which(mt[,2] == 0)## [1] 176 189 282 520 542 629 671 737 738 746 747 748 758 759
## [15] 766 873 901 915 943 950 953 1000 1002 1040 1082 1140 1198 1203
## [29] 1229 1255 1281 1288 1344 1371 1374 1380 1386 1390 1401 1473 1485 1493
## [43] 1494 1723 1882 1895 2043 2059 2177 2179 2180 2181 2182 2183 2184 2185
## [57] 2186 2187 2205 2218 2224 2225 2293 2296 2316 2545 2553 2596 2626 2638
## [71] 2746 2844 2927 2939 2947 2948 2949 2950 2951 2952 2953 2954 2960 3063
## [85] 3198 3389 3465 3519 3539 3623 3779 3784 3785 3874 3896 3914 3918 4131
## [99] 4138 4142 4159 4166 4168 4169 4181 4193 4194 4199 4202 4204 4208 4211
## [113] 4215 4219 4223 4224 4225 4226 4238 4243 4253 4256 4265 4269 4284 4287
## [127] 4288 4322 4341 4409 4437 4527 4559 4593 4628 4640 4655 4715 4726 4748
## [141] 4760 4787 4826 4843 4855 4860 4861 4862 4877 4891 4892 4894 4954 4964
## [155] 5384 5497 5655 5803 5811 5812 5813 5814 5817 5819 5825 5890 5927 5930
## [169] 5937 5961 5962 5974 6014 6170 6207 6254 6313 6321 6380 6391which(mt.f[,2] == 0)## [1] 1354 3950 4044 4127 4188 4506empty = which(mt.f[,2] == 0)
mt.f[empty,]## SubmissionID marking.sec Text Drawing Recording project
## 1493 5474 0 0 0 100 BIOL1040Sem1Report 1
## 4341 8149 0 33 0 67 BIOL1040Sem2Report 2
## 4437 8245 0 0 0 100 BIOL1040Sem2Report 2
## 4527 8335 0 0 0 100 BIOL1040Sem2Report 2
## 4593 8401 0 0 0 100 BIOL1040Sem2Report 2
## 4954 8884 0 100 0 0 BIOL1040Sem2Report 3
## course sem report marking.min
## 1493 BIOL1040 Sem1 Report 1 0
## 4341 BIOL1040 Sem2 Report 2 0
## 4437 BIOL1040 Sem2 Report 2 0
## 4527 BIOL1040 Sem2 Report 2 0
## 4593 BIOL1040 Sem2 Report 2 0
## 4954 BIOL1040 Sem2 Report 3 0mt.f = mt.f[-empty,]
head(mt.f)## SubmissionID marking.sec Text Drawing Recording project
## 1 1487 392 94 0 6 BIOL1040Sem1Report 0
## 2 1486 2931 42 8 50 BIOL1040Sem1Report 0
## 3 1488 1587 50 6 44 BIOL1040Sem1Report 0
## 4 1489 524 0 0 100 BIOL1040Sem1Report 0
## 5 1490 1029 0 0 100 BIOL1040Sem1Report 0
## 6 1491 4203 52 40 7 BIOL1040Sem1Report 0
## course sem report marking.min
## 1 BIOL1040 Sem1 Report 0 6.53
## 2 BIOL1040 Sem1 Report 0 48.85
## 3 BIOL1040 Sem1 Report 0 26.45
## 4 BIOL1040 Sem1 Report 0 8.73
## 5 BIOL1040 Sem1 Report 0 17.15
## 6 BIOL1040 Sem1 Report 0 70.05table(mt.f$project)##
## BIOL1040Sem1Report 0 BIOL1040Sem1Report 1 BIOL1040Sem1Report 2
## 682 673 671
## BIOL1040Sem1Report 3 BIOL1040Sem2Report 1 BIOL1040Sem2Report 2
## 667 897 866
## BIOL1040Sem2Report 3 BIOM2011Sem1Report 1 BIOM2011Sem1Report 2
## 853 239 241
## BIOM2011Sem2Report 1 BIOM2011Sem2Report 2
## 85 80tapply(mt.f[,10], mt.f[,6], mean)## BIOL1040Sem1Report 0 BIOL1040Sem1Report 1 BIOL1040Sem1Report 2
## 17.38937 16.61260 14.00405
## BIOL1040Sem1Report 3 BIOL1040Sem2Report 1 BIOL1040Sem2Report 2
## 11.42303 14.08667 13.53878
## BIOL1040Sem2Report 3 BIOM2011Sem1Report 1 BIOM2011Sem1Report 2
## 11.00701 39.00205 22.24365
## BIOM2011Sem2Report 1 BIOM2011Sem2Report 2
## 27.69329 15.97213tapply(mt.f[,10], mt.f[,6], range)## $`BIOL1040Sem1Report 0`
## [1] 1.35 146.83
##
## $`BIOL1040Sem1Report 1`
## [1] 0.75 160.88
##
## $`BIOL1040Sem1Report 2`
## [1] 0.50 71.38
##
## $`BIOL1040Sem1Report 3`
## [1] 1.15 135.12
##
## $`BIOL1040Sem2Report 1`
## [1] 0.23 107.23
##
## $`BIOL1040Sem2Report 2`
## [1] 0.32 121.68
##
## $`BIOL1040Sem2Report 3`
## [1] 0.05 173.37
##
## $`BIOM2011Sem1Report 1`
## [1] 4.42 304.02
##
## $`BIOM2011Sem1Report 2`
## [1] 2.42 111.40
##
## $`BIOM2011Sem2Report 1`
## [1] 4.35 64.72
##
## $`BIOM2011Sem2Report 2`
## [1] 3.10 68.93default.plot()
hist(log(mt.f[,10]))
my.logs(-3, 6, 3, 9)## -3 -2 -1 0 1 2 3 4
## "3 sec" "8 sec" "22 sec" "1 min" "3 min" "7 min" "20 min" "55 min"
## 5 6
## "2 hr" "7 hr"tapply(log(mt.f[,10]), mt.f[,6], hist)
## $`BIOL1040Sem1Report 0`
## $breaks
## [1] 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0
##
## $counts
## [1] 1 1 17 88 198 205 112 39 17 4
##
## $density
## [1] 0.002932551 0.002932551 0.049853372 0.258064516 0.580645161
## [6] 0.601173021 0.328445748 0.114369501 0.049853372 0.011730205
##
## $mids
## [1] 0.25 0.75 1.25 1.75 2.25 2.75 3.25 3.75 4.25 4.75
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOL1040Sem1Report 1`
## $breaks
## [1] -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5
##
## $counts
## [1] 1 0 3 10 83 218 199 109 34 13 1 2
##
## $density
## [1] 0.002971768 0.000000000 0.008915305 0.029717682 0.246656761
## [6] 0.647845468 0.591381872 0.323922734 0.101040119 0.038632987
## [11] 0.002971768 0.005943536
##
## $mids
## [1] -0.25 0.25 0.75 1.25 1.75 2.25 2.75 3.25 3.75 4.25 4.75
## [12] 5.25
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOL1040Sem1Report 2`
## $breaks
## [1] -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
##
## $counts
## [1] 1 1 1 4 24 134 187 193 94 30 2
##
## $density
## [1] 0.002980626 0.002980626 0.002980626 0.011922504 0.071535022
## [6] 0.399403875 0.557377049 0.575260805 0.280178838 0.089418778
## [11] 0.005961252
##
## $mids
## [1] -0.75 -0.25 0.25 0.75 1.25 1.75 2.25 2.75 3.25 3.75 4.25
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOL1040Sem1Report 3`
## $breaks
## [1] 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0
##
## $counts
## [1] 3 12 46 184 220 127 61 13 0 1
##
## $density
## [1] 0.008995502 0.035982009 0.137931034 0.551724138 0.659670165
## [6] 0.380809595 0.182908546 0.038980510 0.000000000 0.002998501
##
## $mids
## [1] 0.25 0.75 1.25 1.75 2.25 2.75 3.25 3.75 4.25 4.75
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOL1040Sem2Report 1`
## $breaks
## [1] -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0
##
## $counts
## [1] 3 5 2 2 5 52 178 294 192 101 47 14 2
##
## $density
## [1] 0.006688963 0.011148272 0.004459309 0.004459309 0.011148272
## [6] 0.115942029 0.396878484 0.655518395 0.428093645 0.225195095
## [11] 0.104793757 0.031215162 0.004459309
##
## $mids
## [1] -1.25 -0.75 -0.25 0.25 0.75 1.25 1.75 2.25 2.75 3.25 3.75
## [12] 4.25 4.75
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOL1040Sem2Report 2`
## $breaks
## [1] -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0
##
## $counts
## [1] 1 1 1 0 7 42 178 276 212 112 30 5 1
##
## $density
## [1] 0.002309469 0.002309469 0.002309469 0.000000000 0.016166282
## [6] 0.096997691 0.411085450 0.637413395 0.489607390 0.258660508
## [11] 0.069284065 0.011547344 0.002309469
##
## $mids
## [1] -1.25 -0.75 -0.25 0.25 0.75 1.25 1.75 2.25 2.75 3.25 3.75
## [12] 4.25 4.75
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOL1040Sem2Report 3`
## $breaks
## [1] -3 -2 -1 0 1 2 3 4 5 6
##
## $counts
## [1] 4 17 17 34 286 400 88 6 1
##
## $density
## [1] 0.004689332 0.019929660 0.019929660 0.039859320 0.335287222 0.468933177
## [7] 0.103165299 0.007033998 0.001172333
##
## $mids
## [1] -2.5 -1.5 -0.5 0.5 1.5 2.5 3.5 4.5 5.5
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOM2011Sem1Report 1`
## $breaks
## [1] 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 5.5 6.0
##
## $counts
## [1] 1 5 19 45 61 55 41 10 1 1
##
## $density
## [1] 0.008368201 0.041841004 0.158995816 0.376569038 0.510460251
## [6] 0.460251046 0.343096234 0.083682008 0.008368201 0.008368201
##
## $mids
## [1] 1.25 1.75 2.25 2.75 3.25 3.75 4.25 4.75 5.25 5.75
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOM2011Sem1Report 2`
## $breaks
## [1] 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0
##
## $counts
## [1] 1 10 17 42 57 72 34 6 2
##
## $density
## [1] 0.008298755 0.082987552 0.141078838 0.348547718 0.473029046 0.597510373
## [7] 0.282157676 0.049792531 0.016597510
##
## $mids
## [1] 0.75 1.25 1.75 2.25 2.75 3.25 3.75 4.25 4.75
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOM2011Sem2Report 1`
## $breaks
## [1] 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
##
## $counts
## [1] 1 2 7 22 24 26 3
##
## $density
## [1] 0.02352941 0.04705882 0.16470588 0.51764706 0.56470588 0.61176471
## [7] 0.07058824
##
## $mids
## [1] 1.25 1.75 2.25 2.75 3.25 3.75 4.25
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"
##
## $`BIOM2011Sem2Report 2`
## $breaks
## [1] 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5
##
## $counts
## [1] 8 14 18 15 18 6 1
##
## $density
## [1] 0.200 0.350 0.450 0.375 0.450 0.150 0.025
##
## $mids
## [1] 1.25 1.75 2.25 2.75 3.25 3.75 4.25
##
## $xname
## [1] "X[[i]]"
##
## $equidist
## [1] TRUE
##
## attr(,"class")
## [1] "histogram"writing
folder = "./Data Stage 2/"
write.csv(mt.f, file=paste0(folder, "MarkingTimes.csv"), row.names=FALSE)