R语言里的幂律拟合，请高手帮忙

请R语言高手指教帮忙~~
以上是数据做出来h与j的关系图见第一张图，符合幂律分析，我想用在R语言里做幂律分布的直线拟合，具体应该是对数据进行LOG计算，后做双log图，然后用最小二乘进行双直线的拟合,做来出来以下类似第二张的效果图。另外还需要直线的拟合结果公式及R的平方值，和两条直接的交点的值。数据如下：

  h         j
50  298.93877
100  125.24578
150  78.65552
200  55.00881
250  44.20180
300  37.06160
350  29.82777
400  27.38242
450  24.45392
500  21.45078
550  18.57950
600  16.15758
650  16.10872
700  14.23703
750  15.47879
800  14.80924
850   NA