plot(dev_count~pearsonmr,xlab=xl,main="# Developers",ylab="",log="y"); grid()
plot(age/30~pearsonmr,xlab=xl,main="SF.net age (months)",ylab="",log="y"); grid()
a <- cvsage/30
plot(cvsage[a<=91]/30~pearsonmr[a<=91],xlab=xl,main="CVS age (months)",ylab="",log="y"); grid()
plot(size~pearsonmr,xlab=xl,main="Size (SLOC)",ylab="",log="y"); grid()
plot(files~pearsonmr,xlab=xl,main="# Files",ylab="",log="y"); grid()
plot(changes~pearsonmr,xlab=xl,main="# Changes",ylab="",log="y"); grid()
save.history()
savehistory()
save.image()
{
r <- {}
r$pid <- c()
r$r <- c()
r$errors <- c()
r$short <- c()
for(p in pl)
{
query <- paste('select day, nmr from afmodrequest where pid=',p,' order by day asc;')
irregts <- dbGetQuery(con,query)
if(length(irregts$day)>=1)
{
regts <- irregtoregts(irregts$day,irregts$n)
bw=5
regts_s <- ksmooth(c(1:length(regts)),regts,"normal",bandwidth=bw)
if(length(regts_s$y)<364)
{
r$short <- c(r$short,p)
print(paste("Short:",p))
}
else
{
r$pid <- c(r$pid,p)
a <- acf(regts_s$y,plot=FALSE,lag.max=2e5)
y <- a$acf
y <- y[y>0.2]
x <- c(1:length(y))
rpearson <- cor(log(x),log(y))
r$r <- c(r$r,rpearson)
print(paste("Good:",p))
}
}
else
{
print(paste("Error:",p))
r$errors <- c(r$errors,p)
}
}
r
}
modrequest <- globalacfwithr(pids$pid,con)
summary(modrequest)
x<-modrequest$r;plot(density(x))
modrequest$r <- -1*modrequest$r
x<-modrequest$r;plot(density(x,from=0,to=1))
boxplot(x)
globalacfwithr <- function(pl,con)
{
r <- {}
r$pid <- c()
r$r <- c()
r$errors <- c()
r$short <- c()
for(p in pl)
{
query <- paste('select day, nmr from afmodrequest where pid=',p,' order by day asc;')
irregts <- dbGetQuery(con,query)
if(length(irregts$day)>=1)
{
regts <- irregtoregts(irregts$day,irregts$nmr)
bw=5
regts_s <- ksmooth(c(1:length(regts)),regts,"normal",bandwidth=bw)
if(length(regts_s$y)<364)
{
r$short <- c(r$short,p)
print(paste("Short:",p))
}
else
{
r$pid <- c(r$pid,p)
a <- acf(regts_s$y,plot=FALSE,lag.max=2e5)
y <- a$acf
y <- y[y>0.2]
x <- c(1:length(y))
rpearson <- cor(log(x),log(y))
r$r <- c(r$r,rpearson)
print(paste("Good:",p))
}
}
else
{
print(paste("Error:",p))
r$errors <- c(r$errors,p)
}
}
r
}
modrequest <- globalacfwithr(pids$pid,con)
x<-modrequest$r;plot(density(x,from=0,to=1))
summary(modrequest)
summary(modrequest$r)
modrequest$r <- -1*modrequest$r
summary(modrequest$r)
x<-modrequest$r;plot(density(x,from=0,to=1))
x<-fr$pearsonmr;plot(density(x,from=0,to=1))
x<-modrequest$r;plot(density(x,from=0,to=1))
par(mfrow=c(2,1))
x<-fr$pearsonmr;plot(density(x,from=0,to=1))
x<-modrequest$r;plot(density(x,from=0,to=1))
help(cor)
modrequest$r[modrequest$pid=3]
modrequest$r[modrequest$pid==3]
fr[fr$pid==3,]$pearsonmr
fr[fr$pid==3,]$pearson
modrequest$r[modrequest$pid==29779]
fr[fr$pid==29779,]$pearsonmr
fr[fr$pid==29779,]$pearson
globalacfwithr <- function(pl,con)
{
r <- {}
r$pid <- c()
r$r <- c()
r$errors <- c()
r$short <- c()
for(p in pl)
{
query <- paste('select day, nmr from afmodrequest where pid=',p,' order by day asc;')
irregts <- dbGetQuery(con,query)
if(length(irregts$day)>=1)
{
regts <- irregtoregts(irregts$day,irregts$nmr)
bw=5
regts_s <- ksmooth(c(1:length(regts)),regts,"normal",bandwidth=bw)
if(length(regts_s$y)<364)
{
r$short <- c(r$short,p)
print(paste("Short:",p))
}
else
{
r$pid <- c(r$pid,p)
a <- acf(regts_s$y,plot=FALSE,lag.max=2e5)
y <- a$acf
#y <- y[y>0.2]
x <- c(1:length(y))
rpearson <- cor(log(x),log(y))
r$r <- c(r$r,rpearson)
print(paste("Good:",p))
}
}
else
{
print(paste("Error:",p))
r$errors <- c(r$errors,p)
}
}
r
}
modrequest <- globalacfwithr(pids$pid,con)
globalacfwithr <- function(pl,con)
{
r <- {}
r$pid <- c()
r$r <- c()
r$errors <- c()
r$short <- c()
for(p in pl)
{
query <- paste('select day, nmr from afmodrequest where pid=',p,' order by day asc;')
irregts <- dbGetQuery(con,query)
if(length(irregts$day)>=1)
{
regts <- irregtoregts(irregts$day,irregts$nmr)
bw=5
regts_s <- ksmooth(c(1:length(regts)),regts,"normal",bandwidth=bw)
if(length(regts_s$y)<364)
{
r$short <- c(r$short,p)
print(paste("Short:",p))
}
else
{
r$pid <- c(r$pid,p)
a <- acf(regts_s$y,plot=FALSE,lag.max=2e5)
y <- a$acf
y <- y[y>0.1]
x <- c(1:length(y))
rpearson <- cor(log(x),log(y))
r$r <- c(r$r,rpearson)
print(paste("Good:",p))
}
}
else
{
print(paste("Error:",p))
r$errors <- c(r$errors,p)
}
}
r
}
modrequest <- globalacfwithr(pids$pid,con)
x<-modrequest$r;plot(density(x,from=0,to=1))
modrequest$r <- -1*modrequest$r
x<-modrequest$r;plot(density(x,from=0,to=1))
globalacfwithr <- function(pl,con)
{
r <- {}
r$pid <- c()
r$r <- c()
r$errors <- c()
r$short <- c()
for(p in pl)
{
query <- paste('select day, nmr from afmodrequest where pid=',p,' order by day asc;')
irregts <- dbGetQuery(con,query)
if(length(irregts$day)>=1)
{
regts <- irregtoregts(irregts$day,irregts$nmr)
#bw=5
#regts_s <- ksmooth(c(1:length(regts)),regts,"normal",bandwidth=bw)
regts_s$y <- regts
if(length(regts_s$y)<364)
{
r$short <- c(r$short,p)
print(paste("Short:",p))
}
else
{
r$pid <- c(r$pid,p)
a <- acf(regts_s$y,plot=FALSE,lag.max=2e5)
y <- a$acf
y <- y[y>0.1]
x <- c(1:length(y))
rpearson <- cor(log(x),log(y))
r$r <- c(r$r,rpearson)
print(paste("Good:",p))
}
}
else
{
print(paste("Error:",p))
r$errors <- c(r$errors,p)
}
}
r
}
modrequest <- globalacfwithr(pids$pid,con)
globalacfwithr <- function(pl,con)
{
r <- {}
r$pid <- c()
r$r <- c()
r$errors <- c()
r$short <- c()
for(p in pl)
{
query <- paste('select day, nmr from afmodrequest where pid=',p,' order by day asc;')
irregts <- dbGetQuery(con,query)
if(length(irregts$day)>=1)
{
regts <- irregtoregts(irregts$day,irregts$nmr)
#bw=5
#regts_s <- ksmooth(c(1:length(regts)),regts,"normal",bandwidth=bw)
regts <- c()
regts_s$y <- regts
if(length(regts_s$y)<364)
{
r$short <- c(r$short,p)
print(paste("Short:",p))
}
else
{
r$pid <- c(r$pid,p)
a <- acf(regts_s$y,plot=FALSE,lag.max=2e5)
y <- a$acf
y <- y[y>0.1]
x <- c(1:length(y))
rpearson <- cor(log(x),log(y))
r$r <- c(r$r,rpearson)
print(paste("Good:",p))
}
}
else
{
print(paste("Error:",p))
r$errors <- c(r$errors,p)
}
}
r
}
modrequest <- globalacfwithr(pids$pid,con)
globalacfwithr <- function(pl,con)
{
r <- {}
r$pid <- c()
r$r <- c()
r$errors <- c()
r$short <- c()
for(p in pl)
{
query <- paste('select day, nmr from afmodrequest where pid=',p,' order by day asc;')
irregts <- dbGetQuery(con,query)
if(length(irregts$day)>=1)
{
regts <- irregtoregts(irregts$day,irregts$nmr)
#bw=5
#regts_s <- ksmooth(c(1:length(regts)),regts,"normal",bandwidth=bw)
regts_s <- c()
regts_s$y <- regts
if(length(regts_s$y)<364)
{
r$short <- c(r$short,p)
print(paste("Short:",p))
}
else
{
r$pid <- c(r$pid,p)
a <- acf(regts_s$y,plot=FALSE,lag.max=2e5)
y <- a$acf
y <- y[y>0.1]
x <- c(1:length(y))
rpearson <- cor(log(x),log(y))
r$r <- c(r$r,rpearson)
print(paste("Good:",p))
}
}
else
{
print(paste("Error:",p))
r$errors <- c(r$errors,p)
}
}
r
}
modrequest <- globalacfwithr(pids$pid,con)
globalacfwithr <- function(pl,con)
{
r <- {}
r$pid <- c()
r$r <- c()
r$rlog <- c()
r$errors <- c()
r$short <- c()
for(p in pl)
{
query <- paste('select day, nmr from afmodrequest where pid=',p,' order by day asc;')
irregts <- dbGetQuery(con,query)
if(length(irregts$day)>=1)
{
regts <- irregtoregts(irregts$day,irregts$nmr)
#bw=5
#regts_s <- ksmooth(c(1:length(regts)),regts,"normal",bandwidth=bw)
regts_s <- c()
regts_s$y <- regts
if(length(regts_s$y)<364)
{
r$short <- c(r$short,p)
print(paste("Short:",p))
}
else
{
r$pid <- c(r$pid,p)
a <- acf(regts_s$y,plot=FALSE,lag.max=2e5)
y <- a$acf
y <- y[y>0.1]
x <- c(1:length(y))
rpearson <- cor(x,y)
rpearsonlog <- cor(log(x),log(y))
r$r <- c(r$r,rpearson)
r$rlog <- c(r$rlog,pearsonlog)
print(paste("Good:",p))
}
}
else
{
print(paste("Error:",p))
r$errors <- c(r$errors,p)
}
}
r
}
modrequest <- globalacfwithr(pids$pid,con)
globalacfwithr <- function(pl,con)
{
r <- {}
r$pid <- c()
r$r <- c()
r$rlog <- c()
r$errors <- c()
r$short <- c()
for(p in pl)
{
query <- paste('select day, nmr from afmodrequest where pid=',p,' order by day asc;')
irregts <- dbGetQuery(con,query)
if(length(irregts$day)>=1)
{
regts <- irregtoregts(irregts$day,irregts$nmr)
#bw=5
#regts_s <- ksmooth(c(1:length(regts)),regts,"normal",bandwidth=bw)
regts_s <- c()
regts_s$y <- regts
if(length(regts_s$y)<364)
{
r$short <- c(r$short,p)
print(paste("Short:",p))
}
else
{
r$pid <- c(r$pid,p)
a <- acf(regts_s$y,plot=FALSE,lag.max=2e5)
y <- a$acf
y <- y[y>0.1]
x <- c(1:length(y))
rpearson <- cor(x,y)
rpearsonlog <- cor(log(x),log(y))
r$r <- c(r$r,rpearson)
r$rlog <- c(r$rlog,rpearsonlog)
print(paste("Good:",p))
}
}
else
{
print(paste("Error:",p))
r$errors <- c(r$errors,p)
}
}
r
}
modrequest <- globalacfwithr(pids$pid,con)
x<-modrequest$r;plot(density(x,from=0,to=1))
summary(modrequest)
summary(modrequest$r)
modrequest$r<--1*modrequest$r
summary(modrequest$rlog)
modrequest$rlog<--1*modrequest$rlog
x<-modrequest$r;plot(density(x,from=0,to=1))
x<-ommit.na(modrequest$r);plot(density(x,from=0,to=1))
x<-na.omit(modrequest$r);plot(density(x,from=0,to=1))
x<-na.omit(modrequest$rlog);plot(density(x,from=0,to=1))
x<-na.omit(modrequest$r);plot(density(x,from=0,to=1))
savehistory()
x<-na.omit(modrequest$r);plot(density(x,from=0,to=1))
ls()
summary(fr)
x<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1))
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1))
x<-na.omit(r$pearson);plot(density(x,from=0,to=1))
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1))
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1))
lines(fr$pearsonmr)
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1))
lines(fr$pearsonmr,col="blue")
x<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1))
lines(fr$pearson,col="blue")
x<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1))
lines(density(fr$pearson),col="blue")
lines(fr$pearson,col="blue")
lines(density(fr$pearsonmr),col="blue",lty=2)
lines(fr$pearson,col="blue")
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1))
lines(fr$pearson,col="blue")
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1))
lines(density(fr$pearsonmr),col="blue",lty=2)
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1))
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1)); lines(density(fr$pearsonmr,from=0,to=1),col="blue",lty=2)
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1),col="blue"); lines(density(fr$pearsonmr,from=0,to=1),col="red",lwd=2)
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1),col="blue",lwd=3); lines(density(fr$pearsonmr,from=0,to=1),col="red",lwd=3)
x<-na.omit(fr$pearson);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(fr$pearsonmr,from=0,to=1),col="red",lwd=3)
help(legend)
x<-na.omit(fr$pearson);y<-na.ommit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y),from=0,to=1),col="red",lwd=3)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y),from=0,to=1),col="red",lwd=3)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs")
)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05,col=c("blue","red"),lwd=c(3,3))
postscript(pointsize=24)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05,col=c("blue","red"),lwd=c(3,3))
dev.off()
png(width=1024,heigh=768)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05,col=c("blue","red"),lwd=c(3,3))
dev.off()
png(width=800,heigh=600,pointsize=24)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05,col=c("blue","red"),lwd=c(3,3))
dev.off()
png(width=800,heigh=600,pointsize=14)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05,col=c("blue","red"),lwd=c(3,3))
dev.off()
png(width=1000,heigh=650,pointsize=14)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05,col=c("blue","red"),lwd=c(3,3))
dev.off()
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),xlim=c(0.4,1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05,col=c("blue","red"),lwd=c(3,3))
png(width=1000,heigh=650,pointsize=14)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),xlim=c(0.4,1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05,col=c("blue","red"),lwd=c(3,3))
dev.off()
png(width=800,heigh=600,pointsize=12)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),xlim=c(0.4,1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main=""); lines(density(y,from=0,to=1),col="red",lwd=3); legend("topleft",c("Commits","MRs"),inset=.05,col=c("blue","red"),lwd=c(3,3))
dev.off()
png(width=800,heigh=600,pointsize=12)
dev.off()
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),xlim=c(0.4,1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main="")
png(width=800,heigh=600,pointsize=12)
x<-na.omit(fr$pearson);y<-na.omit(fr$pearsonmr);plot(density(x,from=0,to=1),xlim=c(0.4,1),col="blue",lwd=3,xlab="Pearson coefficient",ylab="Density",main="")
dev.off()
x<-na.omit(modrequest$r);plot(density(x,from=0,to=1))
help(quantile)
help(quantile)
quantile(fr$pearson,probs=seq(0,.1,1))
quantile(fr$pearson,probs=seq(0,1,.1))
length(fr$pearson)
quantile(fr$pearsonmr,probs=seq(0,1,.1))
savehistory()