#read_data

data <- read.csv("C:/Users/sabih/RWork/data0157.txt", stringsAsFactors = FALSE)

#sample_selection
data1<-cbind(data[6:15])
midtermdata<-rbind(data1[20:29,])

midtermdatamatrix<-as.matrix(midtermdata)

colnames(midtermdatamatrix)<-c("v1","v2","v3","v4","v5","v6","v7","v8","v9","v10")

#control
is.matrix(midtermdatamatrix)

#n and k 
k<-ncol(midtermdatamatrix)
n<-nrow(midtermdatamatrix)

#dolap:)
output<-matrix(c(rep(0,10*6)),10,6)
colnames(output)<-c("p","q","p*q","sd","total","rel")
rownames(output)<-c(colnames(midtermdatamatrix))

#p values
p.stats<-colMeans(midtermdatamatrix)

#q values
q.stats<-1-p.stats

#total scores
total<-rowSums(midtermdatamatrix)
summary(total)

#filling in the matrix
output[,1]<-round(p.stats,3)
output[,2]<-round(q.stats,3)
output[,3]<-round(p.stats*q.stats,3)
output[,4]<-round(sqrt(p.stats*q.stats),3)
output[,5]<-round(rowSums(midtermdatamatrix),3)
output[,6]<-round(k/(k-1))*(sum(p.stats*q.stats))


#creating tables
values1<-tapply(midtermdata[,1],midtermdata[,10],mean)
values2<-tapply(midtermdata[,2],midtermdata[,10],mean)
values3<-tapply(midtermdata[,3],midtermdata[,10],mean)
values4<-tapply(midtermdata[,4],midtermdata[,10],mean)
values5<-tapply(midtermdata[,5],midtermdata[,10],mean)
values6<-tapply(midtermdata[,6],midtermdata[,10],mean)
values7<-tapply(midtermdata[,7],midtermdata[,10],mean)
values8<-tapply(midtermdata[,8],midtermdata[,10],mean)
values9<-tapply(midtermdata[,9],midtermdata[,10],mean)

#graphs up to now
hist(rowSums(midtermdatamatrix))
hist(p.stats)
hist(q.stats)
hist(total)
hist(p.stats*q.stats)
hist(sqrt(p.stats*q.stats))
hist(rowSums(midtermdatamatrix))

plot(values1,ylim<-c(0,1))
lines(values2)
lines(values3)
lines(values4)
lines(values5)
lines(values6)
lines(values7)
lines(values8)
lines(values9)

plot(p.stats)
plot(q.stats)
plot(total)
plot(p.stats*q.stats)
plot(sqrt(p.stats*q.stats))
plot(rowSums(midtermdatamatrix))

boxplot(total)

pie_stats<-rbind(p.stats,q.stats)
pie(pie_stats[,1],main="p/q ratio", col=rainbow(length(pie_stats[,1])))
pie(pie_stats[,2],main="p/q ratio", col=rainbow(length(pie_stats[,1])))
pie(pie_stats[,3],main="p/q ratio", col=rainbow(length(pie_stats[,1])))
pie(pie_stats[,4],main="p/q ratio", col=rainbow(length(pie_stats[,1])))
pie(pie_stats[,5],main="p/q ratio", col=rainbow(length(pie_stats[,1])))
pie(pie_stats[,6],main="p/q ratio", col=rainbow(length(pie_stats[,1])))
pie(pie_stats[,7],main="p/q ratio", col=rainbow(length(pie_stats[,1])))
pie(pie_stats[,8],main="p/q ratio", col=rainbow(length(pie_stats[,1])))
pie(pie_stats[,9],main="p/q ratio", col=rainbow(length(pie_stats[,1])))
pie(pie_stats[,10],main="p/q ratio", col=rainbow(length(pie_stats[,1])))

install.packages("ggplot2")
library(ggplot2)