Is the fundamental group of every subset of $\mathbb{R}^2$ torsion-free?

I don't have a copy of the 1998 paper by Eda from which the result apparently follows. However, the question can be answered using results from the paper The fundamental groups of subsets of closed surfaces inject into their first shape groups by Fischer & Zastrow, which is freely available on the arXiv. They show that the fundamental group of a compact proper subset of a closed surface is locally free (among other things). A group is called locally free if every finitely generated subgroup is free, which implies that it is torsion-free (the subgroup generated by an element $g\not=1$ is free, hence infinite, so $g$ does not have finite order). From this, we can show the following.

1) If $X$ is a proper subset of a closed surface and $x\in X$, then $\pi_1(X,x)$ is torsion-free.

In particular, $\mathbb{R}^2$ is homeomorphic to the sphere $S^2$ with a point removed, so every $X\subseteq\mathbb{R}^2$ is homemorphic to a proper subset of $S^2$ and $\pi_1(X,x)$ is torsion-free for $x\in X$.

To show that (1) does indeed follow from the stated result, consider a closed curve $\gamma\colon[0,1]\to X$ with $\gamma(0)=\gamma(1)=x$ such that $[\gamma]\in\pi_1(X,x)$ has finite order, so that $[\gamma]^{n}=1$ for some $n\ge1$. It needs to be shown that $[\gamma]=1$. First, writing $\gamma^{n}=\gamma\ast\cdots\ast\gamma$ for the n-fold product of $\gamma$, we have that $\gamma^{n}$ is null-homotopic in $X$. Letting $H\colon[0,1]^2\to X$ be a null-homotopy for $\gamma^{(n)}$ then $Y\equiv{\rm Im}(H)\subseteq X$ is compact and $\gamma^n$ in null-homotopic in $Y$. By the result from the paper, $\pi_1(Y,x)$ is torsion-free, so $\gamma$ is null-homotopic in $Y$ and, hence, in $X$.

Update: Reading through the paper, it is apparent that the following is true:

If $x\in X\subseteq\mathbb{R}^2$ and $\gamma\colon[0,1]\to X$ is a curve with $\gamma(0)=\gamma(1)=x$ which is not null-homotopic in $X$, then there is a finite $F\subseteq\mathbb{R}^2\setminus X$ such that $\gamma$ is not null-homotopic in $\mathbb{R}^2\setminus F$.

This implies the result asked for, since if $\gamma$ is not null-homotopic in $X$, then it is not null-homotopic in $\mathbb{R}^2\setminus F$ for some finite set $F\subseteq\mathbb{R}^2\setminus X$. However, the fundamental group of $\mathbb{R}^2\setminus F$ is freely generated by $\vert F\vert$ elements, so is torsion-free. Therefore, $\gamma^n$ is not null-homotopic in $X\subseteq\mathbb{R}^2\setminus F$.

To prove the quoted result, let $X^\prime\subseteq X$ be the union of $\gamma([0,1])$ and the bounded connected components of $\mathbb{R}^2\setminus\gamma([0,1])$ which are subsets of $X$. Also let $U_1,U_2,\ldots$ be the bounded connected components of $\mathbb{R}^2\setminus\gamma([0,1])$ which do not lie entirely in $X$. If we choose $r > 0$ such that $\gamma([0,1])\subseteq B_r(0)$ and let $B_{r_i}(z_i)\subseteq U_i$ be open balls about points $z_i\not\in X$ then the Riemann mapping theorem can be used to show that $X^{\prime\prime}=\bar B_r(0)\setminus\bigcup_iB_{r_i}(z_i)$ retracts onto $X^\prime$. So, $\gamma$ is not null-homotopic in $X^{\prime\prime}$.

It then needs to be shown that $\gamma$ is not null-homotopic in $X_n=\bar B_r(0)\setminus\bigcup_{i=1}^nB_{r_i}(z_i)$ for a finite $n$. This is the difficult part, and is established in the appendix of the linked paper by using the fact that sets such as $X^{\prime\prime}$ are homeomorphic to something like a Sierpinski carpet (at least, the Sierpinski carpet is the most complex that it can get). A lengthy argument is used, counting crossings of horizontal and vertical strips in the square to show that $\gamma$ cannot be null-homotopic in $X_n$ for every $n$. Then, setting $F=\{z_1,z_2,\ldots,z_n\}$, $\mathbb{R}^2\setminus F$ retracts onto $X_n$, so $\gamma$ is not null-homotopic in $\mathbb{R}^2\setminus F$.

I'll also give a quick explanation of where the result stated in the paper comes from, although I haven't worked through all of the details (yet).

The idea is to embed the fundamental group $\pi_1(X,x)$ into another group $\check{\pi}_1(X,x)$ called the first shape homotopy group, which is calculated as an inverse limit of fundamental groups of simplicial complexes. If $\mathcal{U}$ is a finite open cover of $X$ (or, just locally finite if $X$ is not compact) then you can define the abstract simplicial complex $N(\mathcal{U})$ consisting of subsets $\Delta\subseteq\mathcal{U}$ for which $\bigcap\Delta\not=\emptyset$. Then you can define natural homomorphisms of fundamental groups $\pi_1(X,x)\to\pi_1(N(\mathcal{U}),\{U\})$. Taking the inverse limit over refinements of the partition gives the first shape group and a natural homomorphism $\pi_1(X,x)\to\check{\pi}_1(X,x)$. In general, this does not have to be either injective or surjective. However, for subsets of a closed surface then it is injective, as stated in the title of the paper. It can fail to be injective in three dimensions, and an example of this failure is provided in the paper.

In the case where $X$ is compact, then it can be written as the intersection of a sequence of piecewise-linear submanifolds of the surface with nonempty boundary, which each have finitely generated and free fundamental groups. The first shape group can then be written as the inverse limit of a sequence of finitely generated and free groups. You can show that subgroups of inverse limits of free groups are locally free.