R script for plotting pulse of social activity with pulse function fit

data.csv

plotFit.r

#!/usr/bin/env Rscript

library(ggplot2)
library(gridExtra)



fitData <- read.csv(file="./data.csv", header=TRUE, sep=",")

fitData$dates <- as.POSIXct(fitData$date, format = "%m/%d/%y %H:%M")

summary(fitData)

###

p0 <- qplot(fitData$ts1, fitData$twFit, data=fitData,
        , ylab = "timestamp" 
        , xlab="Number of Tweets/Min"
        , geom="line" ) 

p2 <- qplot(fitData$ts1, fitData$twitter, data=fitData,
        , ylab = "timestamp" 
        , xlab="Number of Tweets/Min" )

p1 <- qplot(fitData$ts1, fitData$ngFit, data=fitData,
        , ylab = "timestamp" 
        , xlab="Number of Articles/Min" 
        , geom="line") 

p3 <- qplot(fitData$ts1, fitData$newsgator, data=fitData,
        , ylab = "timestamp" 
        , xlab="Number of Tweets/Min" )

png(filename = "./fitsraw.png", width = 640, height = 900, units = 'px')
print(
   grid.arrange(p0, p1, p2, p3,  ncol=1)
)
dev.off()

 
p5 <- qplot(fitData$dates, fitData$twitter, data=fitData
         , xlab = "Jan 16, 2012 (GMT)"
         , ylab="Tweets/Min" 
         , geom="bar"
         , fill = "red"
         , main = "Twitter"
         , stat = "identity") + 
        scale_x_datetime(format="%H:%M") +
        opts(legend.position = "none") +
        geom_line(aes(fitData$dates, fitData$twFit), fitData, color="blue")
 
p7 <- qplot(fitData$dates, fitData$newsgator, data=fitData
         , xlab = "Jan 16, 2012 (GMT)"
         , ylab="Articles/Min"
         , geom="bar"
         , fill = "red"
         , main = "Newsgator"
         , stat = "identity") + 
        scale_x_datetime(format="%H:%M") +
        opts(legend.position = "none") +
        geom_line(aes(fitData$dates, fitData$ngFit), fitData, color="blue")
                
png(filename = "./pubfits.png", width = 534, height = 700, units = 'px')
print(
    grid.arrange(p5, p7,  ncol=1)
)  
dev.off()