EJEL 2009 - Free Statistics and Forecasting Software (Calculators) v.1.1.23-r7

All rights reserved. The non-commercial (academic) use of this software is free of charge. The only thing that is asked in return is to cite this software when results are used in publications.

This free online software (calculator) computes the computations that are presented in the article: How Reproducible Research Leads to Non-Rote Learning Within a Socially Constructivist E-Learning Environment, which is published in the Electronic Journal of e-Learning, 2009
Click here to edit the underlying code of this R Module.
Source code of R module
if (par1 == '2007') par1 <- 0
if (par1 == '2008') par1 <- 1
zz <- as.data.frame(read.table(file="http://www.wessa.net/download/newdata2.csv",sep=",",header=T))
z <- zz[zz$Year==par1,]
if (par1 == 0) {
z <- as.data.frame(cbind(z$X1,z$X2,z$X3,z$X4,z$X5,z$X6, z$X7,z$X8,z$X9,z$X10,z$X11,z$X12,z$X13,z$X14, z$X15,z$X16,z$X17,z$X18,z$TX, z$nnzfg,z$Bcount, z$Gender,z$Pop,z$nnzfr))
cz <- c("X1","X2","X3","X4","X5", "X6","X7","X8","X9","X10","X11", "X12","X13","X14","X15","X16","X17", "X18", "TX", "nnzfg","Bcount","Gender","Pop","nnzfr")
}
if (par1 == 1) {
z <- as.data.frame(cbind(z$X1,z$X2,z$X3,z$X4, z$X5,z$X6,z$X7,z$X8,z$X9,z$X10,z$X11,z$X12,z$X13,z$X14, z$X15,z$X16,z$X17, z$X18,z$X19,z$X20,z$X21, z$X22,z$X23,z$X24, z$X25,z$TX, z$nnzfg,z$Bcount,z$Gender,z$Pop,z$nnzfr))
cz <- c("X1","X2","X3","X4","X5","X6","X7","X8","X9","X10","X11","X12","X13","X14","X15","X16","X17","X18", "X19","X20","X21","X22","X23","X24","X25","TX", "nnzfg","Bcount","Gender","Pop","nnzfr")
}
colnames(z) <- cz
if (par1 == 0) {
z[,"TX"] <- z$X1 + z$X2 + z$X3 + z$X4 + z$X5 + z$X6 + z$X7 + z$X8 + z$X9 + z$X10 + z$X11 + z$X12 + z$X13 + z$X14 + z$X15 + z$X16 + z$X17 + z$X18
}
if (par1 == 1) {
z[,"TX"] <- z$X1 + z$X2 + z$X3 + z$X4 + z$X5 + z$X6 + z$X7 + z$X8 + z$X9 + z$X10 + z$X11 + z$X12 + z$X13 + z$X14 + z$X15 + z$X16 + z$X17 + z$X18 + z$X19 + z$X20 + z$X21 + z$X22 + z$X23 + z$X24 + z$X25
}
z <- z[complete.cases(z),]
if (par1 == 0) mythreshold <- 100
if (par1 == 1) mythreshold <- 55
if (par1 == 0) mythreshold1s <- 150
if (par1 == 0) mythreshold2s <- 170
if (par1 == 1) mythreshold1s <- 45
if (par1 == 1) mythreshold2s <- 50
bitmap(file="pic1.png")
boxplot(z[z$Pop==0,"TX"] ~ cut(z[z$Pop==0,"nnzfg"], c(0,mythreshold,450)), notch=T, main="# Submitted Verbal Feedback Messages", xlab="# Messages", ylab="Exam Scores (Bachelor)")
dev.off()
bitmap(file="pic2.png")
boxplot(z[z$Pop==0,"TX"] ~ cut(z[z$Pop==0,"nnzfr"], c(0,mythreshold,450)), notch=T, main="# Received Verbal Feedback Messages", xlab="# Messages", ylab="Exam Scores (Bachelor)")
dev.off()
bitmap(file="pic3.png")
boxplot(z[z$Pop==1,"TX"] ~ cut(z[z$Pop==1,"nnzfg"], c(0,mythreshold1s,450)), notch=T, main="# Submitted Verbal Feedback Messages", xlab="# Messages", ylab="Exam Scores (Switching)")
dev.off()
bitmap(file="pic4.png")
boxplot(z[z$Pop==1,"TX"] ~ cut(z[z$Pop==1,"nnzfr"], c(0,mythreshold2s,450)), notch=T, main="# Received Verbal Feedback Messages", xlab="# Messages", ylab="Exam Scores (Switching)")
dev.off()
load(file="createtable")
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,"Reproducible Computations - Descriptive Statistics",2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"Number of Bachelor students",header=TRUE)
a<-table.element(a,length(z[z$Pop==0,"TX"]))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"Number of Switching students",header=TRUE)
a<-table.element(a,length(z[z$Pop==1,"TX"]))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"Number of submitted feedback communications",header=TRUE)
a<-table.element(a,sum(z[,"nnzfg"]))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"Number of received feedback communications",header=TRUE)
a<-table.element(a,sum(z[,"nnzfr"]))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file="mytable.tab")
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,"Reproducible Computations - Two-dimensional Contingency Table - by population",6,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"Bachelor",header=TRUE)
a<-table.element(a,"# Submitted Verbal Feedback Messages",2,header=TRUE)
a<-table.element(a," ",header=TRUE)
a<-table.element(a,"# Received Verbal Feedback Messages",2,header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"Exam Score",header=TRUE)
a<-table.element(a,paste("(0,",mythreshold,"]"),header=TRUE)
a<-table.element(a,paste("(",mythreshold,",450]"),header=TRUE)
a<-table.element(a," ",header=TRUE)
a<-table.element(a,paste("(0,",mythreshold,"]"),header=TRUE)
a<-table.element(a,paste("(",mythreshold,",450]"),header=TRUE)
a<-table.row.end(a)
summary(z[z$Pop==0,'nnzfg'])
if (par1 == 0) {
(mytg <- table(cut(z[z$Pop==0,"TX"],c(-3,4,7,18)),cut(z[z$Pop==0,"nnzfg"],c(0,mythreshold,450))))
(mychsqtg <- chisq.test(mytg))
(mytr <- table(cut(z[z$Pop==0,"TX"],c(-3,4,7,18)),cut(z[z$Pop==0,"nnzfr"],c(0,mythreshold,450))))
(mychsqtr <- chisq.test(mytr))
}
if (par1 == 1) {
(mytg <- table(cut(z[z$Pop==0,"TX"],c(-3,4,8,25)),cut(z[z$Pop==0,"nnzfg"],c(0,mythreshold,450))))
(mychsqtg <- chisq.test(mytg))
(mytr <- table(cut(z[z$Pop==0,"TX"],c(-3,4,8,25)),cut(z[z$Pop==0,"nnzfr"],c(0,mythreshold,450))))
(mychsqtr <- chisq.test(mytr))
}
a<-table.row.start(a)
a<-table.element(a,"(-3,4]",header=TRUE)
a<-table.element(a,mytg[1,1])
a<-table.element(a,mytg[1,2])
a<-table.element(a," ")
a<-table.element(a,mytr[1,1])
a<-table.element(a,mytr[1,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"(4,7]",header=TRUE)
a<-table.element(a,mytg[2,1])
a<-table.element(a,mytg[2,2])
a<-table.element(a," ")
a<-table.element(a,mytr[2,1])
a<-table.element(a,mytr[2,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"(7,18]",header=TRUE)
a<-table.element(a,mytg[3,1])
a<-table.element(a,mytg[3,2])
a<-table.element(a," ")
a<-table.element(a,mytr[3,1])
a<-table.element(a,mytr[3,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"X-squared",header=TRUE)
a<-table.element(a,mychsqtg$statistic,2)
a<-table.element(a," ")
a<-table.element(a,mychsqtr$statistic,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"degrees of freedom",header=TRUE)
a<-table.element(a,mychsqtg$parameter,2)
a<-table.element(a," ")
a<-table.element(a,mychsqtr$parameter,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"p value",header=TRUE)
a<-table.element(a,mychsqtg$p.value,2)
a<-table.element(a," ")
a<-table.element(a,mychsqtr$p.value,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a," ",6)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"Switching",header=TRUE)
a<-table.element(a,"# Submitted Verbal Feedback Messages",2,header=TRUE)
a<-table.element(a," ",header=TRUE)
a<-table.element(a,"# Received Verbal Feedback Messages",2,header=TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"Exam Score",header=TRUE)
a<-table.element(a,paste("(0,",mythreshold1s,"]"),header=TRUE)
a<-table.element(a,paste("(",mythreshold1s,",450]"),header=TRUE)
a<-table.element(a," ",header=TRUE)
a<-table.element(a,paste("(0,",mythreshold2s,"]"),header=TRUE)
a<-table.element(a,paste("(",mythreshold2s,",450]"),header=TRUE)
a<-table.row.end(a)
summary(z[z$Pop==1,'nnzfg'])
if (par1 == 0) {
(mytg <- table(cut(z[z$Pop==1,"TX"],c(-3,4,7,18)),cut(z[z$Pop==1,"nnzfg"],c(0,mythreshold1s,450))))
(mychsqtg <- chisq.test(mytg))
(mytr <- table(cut(z[z$Pop==1,"TX"],c(-3,4,7,18)),cut(z[z$Pop==1,"nnzfr"],c(0,mythreshold2s,450))))
(mychsqtr <- chisq.test(mytr))
}
if (par1 == 1) {
(mytg <- table(cut(z[z$Pop==1,"TX"],c(-3,4,8,25)),cut(z[z$Pop==1,"nnzfg"],c(0,mythreshold1s,450))))
(mychsqtg <- chisq.test(mytg))
(mytr <- table(cut(z[z$Pop==1,"TX"],c(-3,4,8,25)),cut(z[z$Pop==1,"nnzfr"],c(0,mythreshold2s,450))))
(mychsqtr <- chisq.test(mytr))
}
a<-table.row.start(a)
a<-table.element(a,"(-3,4]",header=TRUE)
a<-table.element(a,mytg[1,1])
a<-table.element(a,mytg[1,2])
a<-table.element(a," ")
a<-table.element(a,mytr[1,1])
a<-table.element(a,mytr[1,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"(4,7]",header=TRUE)
a<-table.element(a,mytg[2,1])
a<-table.element(a,mytg[2,2])
a<-table.element(a," ")
a<-table.element(a,mytr[2,1])
a<-table.element(a,mytr[2,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"(7,18]",header=TRUE)
a<-table.element(a,mytg[3,1])
a<-table.element(a,mytg[3,2])
a<-table.element(a," ")
a<-table.element(a,mytr[3,1])
a<-table.element(a,mytr[3,2])
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"X-squared",header=TRUE)
a<-table.element(a,mychsqtg$statistic,2)
a<-table.element(a," ")
a<-table.element(a,mychsqtr$statistic,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"degrees of freedom",header=TRUE)
a<-table.element(a,mychsqtg$parameter,2)
a<-table.element(a," ")
a<-table.element(a,mychsqtr$parameter,2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,"p value",header=TRUE)
a<-table.element(a,mychsqtg$p.value,2)
a<-table.element(a," ")
a<-table.element(a,mychsqtr$p.value,2)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file="mytable1.tab")
::Free Statistics and Forecasting Software::

:: EJEL 2009 - Free Statistics Software (Calculator) ::
