#!/usr/bin/Rscript 
#
# R-experiment http://tukachev.flogiston.ru/exp.html
# Ðåçóëüòàòû ëîíãèòþäíîãî ýêñïåðèìåíòà â ðåæèìå ðåàëüíîãî âðåìåíè
 
# Èçìåðåíèå âðåìåíè âûïîëíåíèÿ ñêðèïòà
time1 <- Sys.time()
 
# Óâåäîìëåíèå î çàïóñêå ñêðèïòà
system("DISPLAY=:0.0 notify-send -i /home/aelita/Rlogo.png R-stat \"Çàïóñê ñêðèïòà. Îáðàáîòêà äàííûõ ýêñïåðèìåíòà\"")
 
# Çàãðóæàåì áèáëèîòåêè
require(RCurl)
library(ggplot2)
library(grid)
library(psych)
 
labels = c("Ñàìî÷óâñòâèå","Àêòèâíîñòü","Íàñòðîåíèå","Òåìïåðàòóðà")
 
# Çàãðóæàåì äàííûå
myCsv <- getURL("çäåñòü URL íà äîñòóï ê òàáëèöå â Google.Docs &single=true&gid=0&output=csv")
data <- read.table(textConnection(myCsv), header = T, sep = ",", stringsAsFactors=FALSE)
 
# Çàäàåì äàòû äëÿ ãåíåðàöèè ãðàôèêîâ
text1 <- "Ïåðâîå íàáëþäåíèå:"
text2 <- "Ïîñëåäíåå íàáëþäåíèå:"
text3 <- "Äàòà ñîçäàíèÿ:"
data1 <- paste(data$dt[1])
data2 <- paste(data$dt[length(data$dt)])
data3 <- paste(format(Sys.time(), "%d.%m.%Y %H:%M:%S"))
ann1 <- paste(text1,data1)
ann2 <- paste(text2,data2)
ann3 <- paste(text3,data3)
l <- length(data$dt)
 
# Ïðîñòàÿ ëèíåéíàÿ ðåãðåññèÿ
model1 <- lm(data$act~data$temp)
model2 <- lm(data$sam~data$temp)
model3 <- lm(data$nas~data$temp)
m1 <- summary.lm(model1)
m2 <- summary.lm(model2)
m3 <- summary.lm(model3)
r1 <-round((m1$r.squared),3)
r2 <-round((m2$r.squared),3)
r3 <-round((m3$r.squared),3)
rsq1 <- paste("R2=",r1) #Ñîõðàíÿåì êîýôôèöèåíò äåòåðìèíàöèè R2
rsq2 <- paste("R2=",r2) #Ñîõðàíÿåì êîýôôèöèåíò äåòåðìèíàöèè R2
rsq3 <- paste("R2=",r3) #Ñîõðàíÿåì êîýôôèöèåíò äåòåðìèíàöèè R2
 
# Óñòàíàâëèâàåì êàòàëîã äëÿ ñîõðàíåíèÿ ãðàôèêîâ
setwd("/home/aelita/Dropbox/Public/exp")
 
# 1 ãðàôèê (àêòèâíîñòü~òåìïåðàòóðà)
png('act~temp.png', width=600,height=450)
p <- qplot(temp,act,data=data,
           xlab='Òåìïåðàòóðà, Ñ', 
           ylab='Àêòèâíîñòü,\nîöåíêà ïî øêàëå îò 1 äî 9 áàëëîâ', 
           main='Óðîâåíü àêòèâíîñòè â çàâèñèìîñòè \n îò òåìïåðàòóðû òåëà')
p <- p + geom_point(colour = "black", size = 4) 
p <- p + annotate("text", x = 35.7, y = 1.5, label = ann1, colour="green", size=4)
p <- p + annotate("text", x = 35.7, y = 1.2, label = ann2, colour="red", size=4)
p <- p + annotate("text", x = 35.7, y = 0.9, label = ann3, size=4)
p <- p + annotate("text", x = 35.7, y = 1.8, label = rsq1, size=4)
p <- p + geom_segment(aes(x = data$temp[l]+0.1, y = data$act[l]+0.5, 
	xend = data$temp[l], yend = data$act[l]), colour="red", arrow = arrow(length = unit(0.5, "cm")))
p <- p + geom_segment(aes(x = data$temp[1]+0.1, y = data$act[1]+0.5, 
	xend = data$temp[1], yend = data$act[1]), colour="green", arrow = arrow(length = unit(0.5, "cm")))
p <- p + theme_bw(base_size = 14, base_family="Georgia") 
p + geom_smooth(aes(temp, act), method="lm", size = 1, se=TRUE)
dev.off()
 
# 2 ãðàôèê (ñàìî÷óâñòâèå~òåìïåðàòóðà)
png('sam~temp.png', width=600,height=450)
p <- qplot(temp,sam,data=data,
           xlab='Òåìïåðàòóðà, Ñ', 
           ylab='Ñàìî÷óâñòâèå,\nîöåíêà ïî øêàëå îò 1 äî 9 áàëëîâ', 
           main='Óðîâåíü ñàìî÷óâñòâèÿ â çàâèñèìîñòè \n îò òåìïåðàòóðû òåëà')
p <- p + geom_point(colour = "black", size = 4) 
p <- p + annotate("text", x = 35.7, y = 1.5, label = ann1, colour="green", size=4)
p <- p + annotate("text", x = 35.7, y = 1.2, label = ann2, colour="red", size=4)
p <- p + annotate("text", x = 35.7, y = 0.9, label = ann3, size=4)
p <- p + annotate("text", x = 35.7, y = 1.8, label = rsq2, size=4)
p <- p + geom_segment(aes(x = data$temp[l]+0.1, y = data$sam[l]+0.5, 
	xend = data$temp[l], yend = data$sam[l]), colour="red", arrow = arrow(length = unit(0.5, "cm")))
p <- p + geom_segment(aes(x = data$temp[1]+0.1, y = data$sam[1]+0.5, 
	xend = data$temp[1], yend = data$sam[1]), colour="green", arrow = arrow(length = unit(0.5, "cm")))
p <- p + theme_bw(base_size = 14, base_family="Georgia") 
p + geom_smooth(aes(temp, sam), method="lm", size = 1, se=TRUE)
dev.off()
 
# 3 ãðàôèê (íàñòðîåíèå~òåìïåðàòóðà)
png('nas~temp.png', width=600,height=450)
p <- qplot(temp,nas,data=data,
           xlab='Òåìïåðàòóðà, Ñ', 
           ylab='Íàñòðîåíèå,\nîöåíêà ïî øêàëå îò 1 äî 9 áàëëîâ', 
           main='Óðîâåíü íàñòðîåíèÿ â çàâèñèìîñòè \n îò òåìïåðàòóðû òåëà')
p <- p + geom_point(colour = "black", size = 4) 
p <- p + annotate("text", x = 35.7, y = 1.9, label = ann1, colour="green", size=4)
p <- p + annotate("text", x = 35.7, y = 1.6, label = ann2, colour="red", size=4)
p <- p + annotate("text", x = 35.7, y = 1.3, label = ann3, size=4)
p <- p + annotate("text", x = 35.7, y = 2.2, label = rsq3, size=4)
p <- p + geom_segment(aes(x = data$temp[l]+0.1, y = data$nas[l]+0.5, 
	xend = data$temp[l], yend = data$nas[l]), colour="red", arrow = arrow(length = unit(0.5, "cm")))
p <- p + geom_segment(aes(x = data$temp[1]+0.1, y = data$nas[1]+0.5, 
	xend = data$temp[1], yend = data$nas[1]), colour="green", arrow = arrow(length = unit(0.5, "cm")))
p <- p + theme_bw(base_size = 14, base_family="Georgia") 
p + geom_smooth(aes(temp, nas), method="lm", size = 1, se=TRUE)
dev.off()
 
# 4 ãðàôèêà ðàñïðåäåëåíèé
 
png('4plots.png', width=458,height=500)
 
layout(matrix(c(1,2,3,4,5,6,7,8), 2, 2, byrow=TRUE))
 
hist(data$temp, breaks = 9, col="skyblue", border="blue",
     main="Òåìïåðàòóðà",
     xlab="Ãðàäóñû",
     ylab="×àñòîòà")
 
hist(data$act, breaks = 8, col="skyblue", border="blue",
     main="Àêòèâíîñòü",
     xlab="Áàëëû",
     ylab="×àñòîòà")
 
hist(data$sam, breaks = 8, col="skyblue", border="blue",
     main="Ñàìî÷óâñòâèå",
     xlab="Áàëëû",
     ylab="×àñòîòà")
 
hist(data$nas, breaks = 8, col="skyblue", border="blue",
     main="Íàñòðîåíèå",
     xlab="Áàëëû",
     ylab="×àñòîòà")
dev.off()
 
# Ñ÷èòàåì êîððåëÿöèþ
r <-cor(data[-1])
colnames(r) = labels
row.names(r) = labels
 
png('corplot.png', width=458,height=500)
cor.plot(mat.sort(r),TRUE,zlim=c(-1,1), main="Ìàòðèöà èíòåðêîððåëÿöèé")
dev.off()
 
# Ôèêñèðóåì âðåìÿ çàâåðøåíèÿ
time2 <- Sys.time()
runtime <- round((time2 - time1), 3)
t <- paste(runtime, "ñåê")
 
# Óâåäîìëåíèå î çàâåðøåíèè ñêðèïòà
system("DISPLAY=:0.0 notify-send -i /home/aelita/Rlogo.png R-stat \"Ñêðèïò âûïîëíåí\"")