log(Th)
Th <- x/mean(x)
Th
x*log(Th)
y <- (10,10,8,8,8,8,6,6,6,6,6,6,4,4,4,4,2,2)
y <- c(10,10,8,8,8,8,6,6,6,6,6,6,4,4,4,4,2,2)
y
Th1 <- y/mean(y)
Th1
Th1 <- sum(y*log(Th1))
Th1
Th1/sum(y)
sum(y)
mean(y)
mean(y)*18
library(igraph)
library(akima)
library(gam)
library(VGAM)
library(R.methodS3)
library(R.methodsS3)
library(R.oo)
library(digest)
library(R.utils)
library(SparseM)
library(Recompression)
library(Rcompression)
library(R.mathlab)
library(R.matlab)
source("/Users/Markus/Research/SoftwareNature/R-Analysis/plfit.r")
ls()
runtest()
help("plfit")
plfit
plfit()
sample <- c(1.26, 0.34, 0.70, 1.75, 50.57, 1.55, 0.08, 0.42, 0.50, 3.20, 0.15, 0.49, 0.95, 0.24, 1.37, 0.17, 6.98, 0.10, 0.94, 0.38)
sample
unique(sample)
xmins <- sort(unique(sample))
xmins
xmins[-length(xmins)]
xmins <- xmins[-length(xmins)]
xmins
max(sample)
rm(xmin)
rm(xmins)
ls()
sample
mean(sample)
sd(sample)
sample <- sort(sample)
sample
plot(sample)
all <- sum(sample)
all
help("cumsum")
cumsum(sample)
cumsum(sample)/all
cumpercent <- cumsum(sample)/sum
cumpercent <- (cumsum(sample)/sum)
cumsample <- (cumsum(sample)/sum)
cumsum(sample)/all
cumsample <- cumsum(sample)/all
plot(cumsample~sample)
rm(cumsample)
ls
ls()
rep(1/20,times=20)
cumsum(rep(1/20,times=20))
props <- cumsum(rep(1/20,times=20))
plot(props~sample)
plot(props~sample,type=h)
plot(props~sample,type="h")
plot(props~sample,type="l")
plot(props~sample,type="s")
plot(props~sample,type="s",log="x")
dnorm(x)
dnorm(1)
seq(-4,4,length=200)
seq(min(sample),max(sample),length=20)
ttt <- seq(min(sample),max(sample),length=20)
plot(props~ttt,type="s",log="x")
ttt <- seq(mean(sample)-3*sd(sample),mean(sample)+3*sd(sample),length=20)
plot(props~ttt,type="s",log="x")
ls()
props
rep(1/20,times=20)
pnorm(rep(1/20,times=20),mean=mean(sample),sd=sd(sample))
cumsum(pnorm(rep(1/20,times=20),mean=mean(sample),sd=sd(sample)))
cumsum(rep(1,times=20))
qnorm(cumsum(rep(1,times=20)),mean=mean(sample),sd=sd(sample))
rnorm(20,mean=mean(sample),sd=sd(sample))
plot(rnorm(20,mean=mean(sample),sd=sd(sample)))
ttt <- rnorm(20,mean=mean(sample),sd=sd(sample))
plot(props~ttt,type="s")
ttt <- sort(ttt)
plot(props~ttt,type="s")
cumsum(1:10)
cumsum(1:100)
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
load_cluster_data()
ls
ls()
strsplit("C00C100.txt",'.')
unlist(strsplit("C00C100.txt",'.'))
help("strsplit")
unlist(strsplit("C00C100.txt","\."))
unlist(strsplit("C00C100.txt","."))
unlist(strsplit("C00C100.txt",c('.'),fixed=T))
unlist(strsplit("C00C100.txt",c('.'),fixed=T))[1]
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
load_cluster_data()
ls()
systems <- load_cluster_data()
systems
rm(systems)
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
load_cluster_data()
ls
ls()
systems
cluster_data <- load_cluster_data()
ls
ls()
cluster_data
for ( s in systems ) get_system(s)$Clusters
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
cluster_data <- load_cluster_data()
for ( s in systems ) get_system(s,env=cluster_data)$Clusters
for ( s in systems ) get_system(s,cluster_data)$Clusters
get_system("C999C908",cluster_data)
ls()
cluster_data
get("C999C908",cluster_data)
ls(cluster_data)
get("C999C908",cluster_data)
ls(cluster_data)
ls()
ls(cluster_data)
exists("C990C718",cluster_data)
get("C990C718",cluster_data)
cluster_data <- load_cluster_data()
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
cluster_data <- load_cluster_data()
get("C999C908C845C240C935C807C414C979C380C249C138",cluster_data)
exists("C990C718",cluster_data)
get("C990C718",cluster_data)
get("C990C718",cluster_data)$Clusters
get("C999C908C845C240C935C807C414C979C380C249C138",cluster_data)$Clusters
for ( s in systems ) get_system(s,cluster_data)$Clusters
systems
for ( s in systems ) get(s,cluster_data)$Clusters
for ( s in systems ) { print(s) }
for ( s in systems ) { print(s) get(s,cluster_data) }
for ( s in systems ) { print(s); get(s,cluster_data) }
get("C999C908",cluster_data)
get_system("C999C908",cluster_data)
get_system("C999C908",cluster_data)$Clusters
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
gini(get_system("C999C908",cluster_data)$Clusters)
mean(get_system("C999C908",cluster_data)$Clusters)
sd(get_system("C999C908",cluster_data)$Clusters)
mean
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
summarize_cluster("C999C908",cluster_data)
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
summarize_cluster("C999C908",cluster_data)
get_system("C999C908",cluster_data)$Clusters
get_system("C999C908C845C240C935C807C414",cluster_data)$Clusters
summarize_cluster("C999C908C845C240C935C807C414",cluster_data)
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
summarize_all(systems,cluster_data)
c(1,2,3)[1]
c(1,2,3)[2]
c(1,2,3)[3]
c(1,2,3)[0]
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
summarize_all(systems,cluster_data)
all <- summarize_all(systems,cluster_data)
plot(all$Gini~all$N)
all$N
all
all[0]
all[1]
all[2]
source("/Users/Markus/Research/Papers/WCSI2010/Data/Interactions/LoadDataCluster.r")
all <- summarize_all(systems,cluster_data)
all
all$N
plot(all$Gini~all$N)
plot(all$N~all$Gini)
plot(all$Max~all$Min)
plot(all$Max~all$Gini)
plot(all$Min~all$Gini)
plot(all$Mean~all$Gini)
plot(all$Mean~all$N)
plot(all$N~all$Gini)
all$N
plot(all$N~all$Gini,log="y")
plot(all$Max~all$Gini,log="y")
plot(all$Max~all$Gini)
plot(all$Max~all$Gini,log="y")
all$Max
plot(all$Max~all$Gini,log="y")
plot(all$Mean~all$Gini,log="y")
plot(all$StdD~all$Gini,log="y")
plot(all$StdD~all$Gini)
plot(log10(all$StdD)~all$Gini)
x <- c(1,1,1,1,1,2,3,4,1,2,2,3,3,4,5,6,7)
hist(x)
help(hist)
help(hist,freq=F)
help(hist,freq=FALSE)
hist(x,freq=F)
table(x)
sum(x)
table(x)/sum(x)
length(x)
x <- c(1,3,4,7,7,12,14,20,23,43,67,88)
x
table(x)
table(x)/12
sum(table(x)/12)
sum
source("/Users/Markus/Research/SoftwareNature/Menzies/R-SCRIPTS/loadmerged.r")
source("/Users/Markus/Research/SoftwareNature/Menzies/R-SCRIPTS/Gini.r")
loadmerged()
source("/Users/Markus/Research/SoftwareNature/Menzies/R-SCRIPTS/loadmerged.r")
loadmerged()
ls()
ant-1.3
"ant-1.3"
"aa" + "bb"
rm("ant-1.3")
rm("ant-1.4")
rm("ant-1.5")
rm("ant-1.6")
rm("ant-1.7")
ls()
rm("lucene-2.0")
rm("lucene-2.2")
rm("lucene-2.4")
rm("xalan-2.4")
rm("xalan-2.5")
rm("xalan-2.6")
rm("xalan-2.7")
ls()
source("/Users/Markus/Research/SoftwareNature/Menzies/R-SCRIPTS/loadmerged.r")
loadmerged()
ls()