1 Load data

Here we are going to analyse the Experimental data packages of Bioconductor. See the home of the analysis here were we already transormated the data. From that stats we are going to analyse the software category:

load("stats.RData")
stats <- stats[Category == "Experimental", ]
yearly <- yearly[Category == "Experimental", ]
stats
##                     Package                Date Month Year     Category
##     1:              SNAData 2017-02-01 01:00:00    02 2017 Experimental
##     2:              SNAData 2017-04-01 02:00:00    04 2017 Experimental
##     3:              SNAData 2017-05-01 02:00:00    05 2017 Experimental
##     4:              SNAData 2017-08-01 02:00:00    08 2017 Experimental
##     5:              SNAData 2016-02-01 01:00:00    02 2016 Experimental
##     6:              SNAData 2016-03-01 01:00:00    03 2016 Experimental
##     7:              SNAData 2016-04-01 02:00:00    04 2016 Experimental
##     8:              SNAData 2016-05-01 02:00:00    05 2016 Experimental
##     9:              SNAData 2016-09-01 02:00:00    09 2016 Experimental
##    10:              SNAData 2016-11-01 01:00:00    11 2016 Experimental
##    ---                                                                 
## 19590: sampleClassifierData 2017-09-01 02:00:00    09 2017 Experimental
## 19591: sampleClassifierData 2017-10-01 02:00:00    10 2017 Experimental
## 19592: sampleClassifierData 2017-11-01 01:00:00    11 2017 Experimental
## 19593: sampleClassifierData 2017-12-01 01:00:00    12 2017 Experimental
## 19594:         topdownrdata 2017-10-01 02:00:00    10 2017 Experimental
## 19595:         topdownrdata 2017-11-01 01:00:00    11 2017 Experimental
## 19596:         topdownrdata 2017-12-01 01:00:00    12 2017 Experimental
## 19597:           vulcandata 2017-10-01 02:00:00    10 2017 Experimental
## 19598:           vulcandata 2017-11-01 01:00:00    11 2017 Experimental
## 19599:           vulcandata 2017-12-01 01:00:00    12 2017 Experimental
##        Nb_of_distinct_IPs Nb_of_downloads
##     1:                 22              25
##     2:                 24              32
##     3:                 39              53
##     4:                 12              12
##     5:                 37              71
##     6:                 61              83
##     7:                 77             131
##     8:                 57              70
##     9:                 25              47
##    10:                 24              26
##    ---                                   
## 19590:                  6               9
## 19591:                 18              42
## 19592:                 13              57
## 19593:                 16              61
## 19594:                  6              18
## 19595:                  5              41
## 19596:                 15              45
## 19597:                  7              22
## 19598:                  5              10
## 19599:                 11              16

There have been 351 Experimental packages in Bioconductor. Some have been added recently and some later.

2 Packages

2.1 Number

First we explore the number of packages being downloaded by month:

theme_bw <- theme_bw(base_size = 16)
scal <- scale_x_datetime(date_breaks = "3 months")
theme <- theme(axis.text.x = element_text(angle = 60, hjust = 1))
ggplot(stats[, .(Downloads = .N), by = Date], aes(Date, Downloads)) +
  geom_bar(stat = "identity") + 
  theme_bw +
  ggtitle("Packages downloaded") +
  theme + 
  scal + 
  xlab("")
Packages in Bioconductor with downloads

Figure 1: Packages in Bioconductor with downloads

The number of packages being downloaded is increasing with time almost exponentially. Partially explained with the incorporation of new packages

ggplot(stats[, .(Number = sum(Nb_of_downloads)), by = Date], aes(Date, Number)) +
  geom_bar(stat = "identity") + 
  theme_bw +
  ggtitle("Downloads") +
  scal +
  theme +
  xlab("")
Downloads of packages

Figure 2: Downloads of packages

Even if the number of packages increase exponentially, the number of the downloads from 2011 grows linearly with time. Which indicates that each time a software package must compete with more packages to be downloaded.

pd <- position_dodge(0.1)
ggplot(stats[, .(Number = mean(Nb_of_downloads),
                 sem = sd(Nb_of_downloads)/sqrt(.N)),
              by = Date], aes(Date, Number)) +
  geom_errorbar(aes(ymin = Number - sem, ymax = Number + sem), width = .1, position = pd) +
  geom_point() + 
  geom_line() +
  theme_bw +
  ggtitle("Downloads") +
  ylab("Mean download for a package") +
  scal +
  theme + 
  xlab("")
Downloads of packages per package. The error bar indicates the 95% confidence interval.

Figure 3: Downloads of packages per package
The error bar indicates the 95% confidence interval.

Here we can apreciate that the number of downloads per package hasn’t changed much with time. If something, now there is less dispersion between packages downloads.

2.2 Incorporations

This might be due to an increase in the usage of packages or that new packages bring more users. We start knowing how many packages has been introduced in Bioconductor each month.

today <- base::date()
year <- substr(today, 21, 25)
month <- monthsConvert(substr(today, 5, 7))
incorporation <- stats[ , .SD[which.min(Date)], by = Package, .SDcols = "Date"]
histincorporation <- incorporation[, .(Number = .N), by = Date, ]
ggplot(histincorporation, aes(Date, Number)) + 
  geom_bar(stat="identity") + 
  theme_bw + 
  ggtitle("Packages with first download") +
  scal +
  theme +
  xlab("")
New packages

Figure 4: New packages

We can see that there were more than 60 packages before 2009 in Bioconductor, and since them occasionally there is a raise to 10 new downloads (Which would be new packages being added).

ggplot(histincorporation, aes(Date, Number)) + 
  geom_bar(stat="identity") + 
  theme_bw + 
  ggtitle("Packages with first download") +
  scal +
  theme +
  xlab("") + 
  ylim(c(0, 20))
## Warning: Removed 1 rows containing missing values (position_stack).
New packages

Figure 5: New packages

Close view to the new packages not previously downloaded. ## Removed

Using a similar procedure we can approximate the packages deprecated and removed each month. In this case we look for the last date a package was downloaded, excluding the current month:

deprecation <- stats[, .SD[which.max(Date)], by = Package, .SDcols = c("Date",  "Year", "Month")]
deprecation <- deprecation[Month != month & Year == Year, , .SDcols = "Date"] # Before this month
histDeprecation <- deprecation[, .(Number = .N), by = Date, ]
ggplot(histDeprecation, aes(Date, Number)) + 
  geom_bar(stat = "identity") + 
  theme_bw + 
  ggtitle("Packages without downloads") +
  scal +
  theme + 
  ylab("Last seen packages") +
  xlab("")
Date where a package was last downloaded. Aproximates to the date when packages were removed from Bioconductor.

Figure 6: Date where a package was last downloaded
Aproximates to the date when packages were removed from Bioconductor.

Here we can see the packages whose last download was in certain month, assuming that this means they are deprecated. It can happen that a package is no longer downloaded but is still in Bioconductor repository, this would be the reason of the spike to 3000 packages as per last month. In total there are 328 packages downloaded. We further explore how many time between the incorporation of the package and the last download.

df <- merge(incorporation, deprecation, by = "Package")
timeBioconductor <- unclass(df$Date.y-df$Date.x)/(60*60*24*365) # Transform to years
hist(timeBioconductor, main = "Time in Bioconductor", xlab = "Years")
abline(v = mean(timeBioconductor), col = "red")
abline(v = median(timeBioconductor), col = "green")
Time of packages between first and last download

(#fig:time.package)Time of packages between first and last download

Packages tend to stay up to 10 years. Not surprisingly the number of packages incorporated before 2009 and still in the repository are of 0 packages. But those packages not removed how do they do in Bioconductor?

3 Packages downloads

3.1 Ratio downloads per IP

We can start comparing the number of downloads (different from 0) by how many IPs download each package.

pd <- position_dodge(0.1)
ggplot(stats[, .(Number = mean(Nb_of_downloads/Nb_of_distinct_IPs),
                 sem = sd(Nb_of_downloads/Nb_of_distinct_IPs)/sqrt(.N)), 
              by = c("Date")], aes(Date, Number)) +
  geom_point() +
  geom_errorbar(aes(ymin = Number - sem, 
                    ymax = Number + sem),
                width = .1, position = pd) +
  geom_line() +
  theme_bw +
  ggtitle("Downloads per IP") +
  ylab("Mean downloads per IP") +
  xlab("") + 
  theme +
  scal
Downloads and distinct IPs of all months and packages. Each color is a package, the black line represents 2 downloads per IP.

Figure 7: Downloads and distinct IPs of all months and packages
Each color is a package, the black line represents 2 downloads per IP.

Not surprisingly most of the package has two downloads from the same IP, one for each Bioconductor release (black line). However, there are some packages where few IPs download many times the same package, which may indicate that these packages are mostly installed in a few locations.

ratio <- stats[, .(Mean = mean(Nb_of_downloads/Nb_of_distinct_IPs),
                   sem = sd(Nb_of_downloads/Nb_of_distinct_IPs)/sqrt(.N),
                   max = max(Nb_of_downloads/Nb_of_distinct_IPs),
                   min = min(Nb_of_downloads/Nb_of_distinct_IPs)), 
              by = c("Package")]
ratio <- ratio[order(Mean, decreasing = TRUE), ]
ratio$Package <- as.character(ratio$Package)
ratio
##                          Package       Mean        sem        max   min
##   1:          HarmonizedTCGAData 106.457143 33.8354280 170.142857 54.80
##   2:               restfulSEData  12.114286  4.6380525  21.214286  2.00
##   3:                topdownrdata   4.733333  1.7333333   8.200000  3.00
##   4:                methyvimData   3.723856  1.3083869   6.333333  2.25
##   5:                msPurityData   3.602444  1.6449769  36.000000  1.00
##   6:                PtH2O2lipids   3.219956  0.8444034  18.500000  1.00
##   7:                     plasFIA   2.593908  0.3930521   6.222222  1.00
##   8:                   mvoutData   2.293581  0.3303490  22.222222  1.00
##   9:             OnassisJavaLibs   2.263312  0.4123213   3.857143  1.20
##  10:                  hapmapsnp6   2.254495  0.4074776  34.229730  1.00
##  ---                                                                   
## 342:                         all   1.000000         NA   1.000000  1.00
## 343:           breastcancermainz   1.000000         NA   1.000000  1.00
## 344:        breastcancertransbig   1.000000         NA   1.000000  1.00
## 345:            eatonetalchipseq   1.000000         NA   1.000000  1.00
## 346:                ecolileucine   1.000000         NA   1.000000  1.00
## 347:                  encodnasei   1.000000         NA   1.000000  1.00
## 348:                    neve2006   1.000000         NA   1.000000  1.00
## 349:             rnainteractmapk   1.000000         NA   1.000000  1.00
## 350: rnaseqdata.hnrnpc.bam.chr14   1.000000         NA   1.000000  1.00
## 351:                  HMP16SData   1.000000         NA   1.000000  1.00

We can see that the package with more downloads from the same IP is HarmonizedTCGAData, followed by, restfulSEData, topdownrdata and the forth one is methyvimData.

Now we explore if there is some seasons cycles in the downloads, as in figure 2 seems to be some cicles.

3.2 By date

First we can explore the number of IPs per month downloading each package:

ggplot(stats, aes(Date, Nb_of_distinct_IPs, col = Package)) + 
  geom_line() + 
  theme_bw +
  ggtitle("IPs") +
  ylab("Distinct IP downloads") +
  scal +
  theme + 
  guides(col = FALSE)
Distinct IP per package

Figure 8: Distinct IP per package

As we can see there are two groups of packages at the 2009 years, some with low number of IPs and some with bigger number of IPs. As time progress the number of distinct IPs increases for some packages. But is the spread in IPs associated with an increase in downloads?

ggplot(stats, aes(Date, Nb_of_downloads, col = Package)) + 
  geom_line() + 
  theme_bw +
  ggtitle("Downloads per IP") +
  ylab("Downloads") +
  scal +
  theme + 
  guides(col = FALSE)
Downloads per year

Figure 9: Downloads per year

Surprisingly some package have a big outburst of downloads to 400k downloads, others to just 100k downloads. But lets focus on the lower end:

ggplot(stats, aes(Date, Nb_of_downloads, col = Package)) + 
  geom_line() + 
  theme_bw +
  ggtitle("Downloads per package every three months") +
  ylab("Downloads") +
  scal +
  ylim(0, 10000)+
  theme + 
  guides(col = FALSE)
Downloads per year

Figure 10: Downloads per year

There are many packages close to 0 downloads each month, but most packages has less than 10000 downloads per month:

ggplot(stats, aes(Date, Nb_of_downloads, col = Package)) + 
  geom_line() + 
  theme_bw+
  ggtitle("Downloads per package every three months") +
  ylab("Downloads") +
  scal +
  ylim(0, 2500)+
  theme +
  guides(col = FALSE)
## Warning: Removed 39 rows containing missing values (geom_path).
Downloads per year

Figure 11: Downloads per year

As we can see, in general the month of the year also influences the number of downloads. So we have that from 2010 the factors influencing the downloads are the year, and the month.

Maybe there is a relationship between the downloads and the number of IPs per date

ggplot(stats, aes(Date, Nb_of_downloads/Nb_of_distinct_IPs, col = Package)) + 
  geom_line() + 
  theme_bw +
  ggtitle("IPs") +
  ylab("Ratio") +
  scal +
  theme +
  guides(col = FALSE)
Ratio downloads per IP per package

Figure 12: Ratio downloads per IP per package

We can see some packages have ocasional raises of downloads per IP. But for small ranges we miss a lot of packages:

ggplot(stats, aes(Date, Nb_of_downloads/Nb_of_distinct_IPs, col = Package)) + 
  geom_line() + 
  theme_bw +
  ggtitle("IPs") +
  ylab("Ratio") +
  scal +
  theme +
  guides(col = FALSE) +
  ylim(1, 5)
Ratio downloads per IP per package

Figure 13: Ratio downloads per IP per package

But most of the packages seem to be more or less constant and around 2. ## By year To measure if the same IP has downloaded the same package in several months we can compare if the IPs per year are closer to the IP per month.

staty <- stats[, .(Nb_of_distinct_IPs = sum(Nb_of_distinct_IPs)), by = c("Year", "Package")]
year <- staty[yearly, , on = c("Package", "Year")]
year <- year[, Repeated_IP := Nb_of_distinct_IPs-i.Nb_of_distinct_IPs, by = c("Package", "Year")]
year <- year[, Repeated_IP_per := Repeated_IP/Nb_of_distinct_IPs*100, by = c("Package", "Year")]
ggplot(year, aes(Year, Repeated_IP_per)) + 
  geom_violin() +
  xlab("") +
  ylab("Percentatge of repeated IPs") +
  theme_bw +
  guides(col = FALSE)

This gives an idea that most of the IPs don’t update all their packages to the newest version of each package, or it could be that as most of the IPs are dynamic they change the measure is not completely reliable

year2 <- year[, .(m = mean(Repeated_IP_per), sem = sd(Repeated_IP_per)/sqrt(.N)), by = "Year"]
ggplot(year2, aes(Year, m)) +
  geom_bar(stat = "identity") +
  geom_errorbar(aes(ymin = m - sem, ymax = m + sem), 
                width = .1, position = pd) + 
  ylab("Percentatge") +
  xlab("") +
  ggtitle("Mean download per ") + 
  theme_bw
Mean percentage of updates of the packages. The error bar are the SEM.

Figure 14: Mean percentage of updates of the packages
The error bar are the SEM.

We can see however that the percentatge of IPs which download the package quite similar for the packages in Bioconductor. It remains also around 17% of the total IPs.

year$Year <- as.numeric(year$Year)
d <- date()
thisYear <- as.numeric(substr(d, nchar(d)-3, nchar(d)))
ggplot(year, aes(Year, Repeated_IP_per, col = Package)) +
  geom_point() +
  geom_line() +
  guides(col = FALSE) + 
  scale_x_continuous(breaks = seq(2009, thisYear, 1)) + 
  ylab("Percentatge") +
  ggtitle("Update from the same IP") +
  theme_bw
Percentatge of updates from the same IP by package.

Figure 15: Percentatge of updates from the same IP by package

We can see that some packages ocassionaly raise above the mean.

4 Models

4.1 Position in Bioconductor

We can observe if the packages has been consistently the most downloaded package of Bioconductor:

PercDate <- stats[, .(Package, Downloads = Nb_of_downloads/sum(Nb_of_downloads)), by = Date]
PercDate <- PercDate[order(Date, order(Downloads)), ]
ggplot(PercDate, aes(Date, Downloads, col = Package)) + 
  geom_line() + 
  theme_bw + 
  ggtitle("Downloads of packages by dates") +
  xlab("") + 
  ylab("% of downloads") +
  guides(col = FALSE) +
  scal + 
  theme
Downloads of packages in Bioconductor.

Figure 16: Downloads of packages in Bioconductor

We can see that usually the percentatge of downloads of a packages doesn’t reach 5% of downloads of the category. And that there is a huge differences between packages in the number of downloads. Let’s compare the most downloaded package to the other package:

OrdDate <- PercDate[, .(Package, Ord = Downloads/max(Downloads)), by = Date]
ggplot(OrdDate, aes(Date, Ord, col = Package)) + 
  geom_line() + 
  theme_bw + 
  ggtitle("Downloads of packages by dates") +
  xlab("") + 
  ylab("% of the Package more downloaded") +
  guides(col = FALSE) +
  scal + 
  theme
Position of packages in Bioconductor

Figure 17: Position of packages in Bioconductor

We can observe that usually there are few package that are closer to the most downloaded packages and lots of them are far from it. We can see the rank between them to be able to see the evolution of downloads along time of a package:

rankDate <- OrdDate[, .(Package, rank = rank(Ord)/.N), by = Date]
ggplot(rankDate, aes(Date, rank, col = Package)) + 
  geom_line() + 
  theme_bw + 
  ggtitle("Rank of packages by downloads") +
  xlab("") + 
  ylab("Position by downloads") +
  guides(col = FALSE) +
  scal + 
  theme
% of packages in Bioconductor. Closer to 1 indicates the top downloaded packages.

Figure 18: % of packages in Bioconductor
Closer to 1 indicates the top downloaded packages.

Only if we select to follow some package we can track them:

packages <- c("ALL", "airway", "Affyhgu133A2Expr", "GSVAdata", 
              "TCGAWorkflowData", "affydata", "RTCGA.rnaseq", 
              "SomaticCancerAlterations", "Fletcher2013a")
ggplot(rankDate[Package %in% packages, ], aes(Date, rank, col = Package))+ 
  geom_line() + 
  theme_bw + 
  ggtitle("Relative rank of packages by downloads") +
  xlab("") + 
  ylab("Position by downloads") +
  scal + 
  theme
Evolution of downloads of packages in Bioconductor

Figure 19: Evolution of downloads of packages in Bioconductor

Here we can see that ALL and affydata are on the top downloaded packages since 2009. Recently airway has grown to reach the top downloaded packages, while other packages haven’t reached the top 25% packages by downloads on a month.

4.2 Package cycle

We can also observe when did a package reach the maximum number of downloads:

PercPack <- stats[, .(Date, Downloads = Nb_of_downloads/sum(Nb_of_downloads)), by = Package]
ggplot(PercPack, aes(Date, Downloads, col = Package)) + 
  geom_line() + 
  theme_bw + 
  ggtitle("Growth of the packages") +
  xlab("") + 
  ylab("Downloads/max(Downloads)") +
  guides(col = FALSE) +
  scal + 
  theme
Cycle of packages. Percentatge of downloads of a package along time.

Figure 20: Cycle of packages
Percentatge of downloads of a package along time.

Here we can see the date when a package reached the highest number of downloads. For most packages the higher downloads are on the recent months.

OrdPack <- PercPack[, .(Date, rank = Downloads/max(Downloads)), by = Package]
ggplot(OrdPack, aes(Date, rank, col = Package)) + 
  geom_line() + 
  theme_bw + 
  ggtitle("Growth of the packages") +
  xlab("") + 
  ylab("Downloads/max(Downloads)") +
  guides(col = FALSE) +
  scal + 
  theme

Here we can see when the package hast the most downloads along time. As usuall we need to focus on fewer packages to be able to distinguish between package to see their evolution:

ggplot(OrdPack[Package %in% c(packages, "RTools4TB", "SemSim"), ], 
       aes(Date, rank, col = Package)) + 
  geom_line() + 
  theme_bw + 
  ggtitle("Growth of the packages") +
  xlab("") + 
  ylab("Downloads/max(Downloads)") +
  scal + 
  theme
Cycle of few packages. Position of package downloads respect the maximum downloads of the packages along time.

Figure 21: Cycle of few packages
Position of package downloads respect the maximum downloads of the packages along time.

As expected the packages that keep up with Bioconductor growth have a peak near the end of the serie.

We can combine both models into a single one to have a timeless comparison of the packages. So we will know if the max position reached in Bioconductor is done when there are more downloads.

model <- rankDate[OrdPack, on = c("Package", "Date")]
setnames(model, "i.rank", "rank.P")
setnames(model, "rank", "rank.B")
ggplot(model[Package %in% packages], aes(rank.P, rank.B, col = Package)) +
  geom_point() + 
  theme_bw + 
  geom_line() +
  ylab("Position in Bioconductor") + 
  xlab("Position in package")
Relative position of the package. Position in Bioconductor downloads and position of package downloads.

Figure 22: Relative position of the package
Position in Bioconductor downloads and position of package downloads.

This is a timeless comparison of the growth of packages in Bioconductor and against themselfs. The higher on the y-axis the higher they have been on Bioconductor, and the higher on the x-axis the most downloads they had.

Next we add the time as a factor so we standarize by the time a package has been in Bioconductor. We have omitted months with 0 downloads of a package so this might bias a bit the ranking.

model[, rank.D := rank(as.numeric(Date))/.N, by = Package]
ggplot(model[Package %in% packages], aes(rank.D, rank.B, col = Package)) +
  geom_point() + 
  theme_bw + 
  geom_line() +
  xlab("By date") + 
  ylab("Position in Bioconductor")
Relative position of the packages. Position relative to Bioconductor and itself compared to date.

Figure 23: Relative position of the packages
Position relative to Bioconductor and itself compared to date.

  
ggplot(model[Package %in%  packages], aes(rank.D, rank.P, col = Package)) +
  geom_point() + 
  theme_bw + 
  geom_line() +
  xlab("By date") + 
  ylab("Position in Package")
Relative position of the packages. Position relative to Bioconductor and itself compared to date.

Figure 23: Relative position of the packages
Position relative to Bioconductor and itself compared to date.

We can observe that some packages climb to the top of Bioconductor sooner than others and some climb and then decrease the number of downloads proportional to the the number of downloads in Bioconductor. # SessionInfo {.unnumbered}

sessionInfo()
## R version 3.4.3 (2017-11-30)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 16.04.3 LTS
## 
## Matrix products: default
## BLAS: /usr/lib/libblas/libblas.so.3.6.0
## LAPACK: /usr/lib/lapack/liblapack.so.3.6.0
## 
## locale:
##  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
##  [3] LC_TIME=es_ES.UTF-8        LC_COLLATE=en_US.UTF-8    
##  [5] LC_MONETARY=es_ES.UTF-8    LC_MESSAGES=en_US.UTF-8   
##  [7] LC_PAPER=es_ES.UTF-8       LC_NAME=C                 
##  [9] LC_ADDRESS=C               LC_TELEPHONE=C            
## [11] LC_MEASUREMENT=es_ES.UTF-8 LC_IDENTIFICATION=C       
## 
## attached base packages:
## [1] stats     graphics  grDevices utils     datasets  methods   base     
## 
## other attached packages:
## [1] data.table_1.10.4-3 ggplot2_2.2.1       BiocStyle_2.6.1    
## 
## loaded via a namespace (and not attached):
##  [1] Rcpp_0.12.14     knitr_1.17       magrittr_1.5     munsell_0.4.3   
##  [5] colorspace_1.3-2 rlang_0.1.4      highr_0.6        stringr_1.2.0   
##  [9] plyr_1.8.4       tools_3.4.3      grid_3.4.3       gtable_0.2.0    
## [13] htmltools_0.3.6  yaml_2.1.16      lazyeval_0.2.1   rprojroot_1.2   
## [17] digest_0.6.13    tibble_1.3.4     bookdown_0.5     evaluate_0.10.1 
## [21] rmarkdown_1.8    labeling_0.3     stringi_1.1.6    compiler_3.4.3  
## [25] scales_0.5.0     backports_1.1.2