join was slower, so used R