Just add a third column with link strength to the association matrix plotCophylo2 [Read more...]

WelshPerm [Read more...]

mat.sort [Read more...]

R only does natural (neperian) log scales by default, and this is lame. Here is a simple code to do decimal log scale, pretty much a requirement for scientists… The force(x/y)lim options works for natural and log scales (for the later case, you need to specify the ... [Read more...]

Producing this kind of graphs (below) in R can be a pain in the a*s. Here is a simple code that requires that data are presented in lists (see the example below). multbar [Read more...]

The GenEstim function presented here uses a very simple genetic algorithm to estimate parameters. The function returns the best estimated set of parameters ($estim), the AIC ($information) at each generation, and the cost of the best model ($bestcost) at each generation. Results of running the program with a logistic function : ... [Read more...]

A quick workaround to have a filled.contour plot with natural log10-log10 scale (instead of the default natural log scale) plotmat [Read more...]

Here is a script that colors a 3D plot according to the z value (height) of each point. persp.withcol [Read more...]

I spend some time looking for an algorithm to find local extrema in a vector (time series). The solution I used is to “walk” through the vector by step larger than 1, in order to retain only one value even when the values are very noisy (see the picture at the ... [Read more...]