I think ACuriousMind's answer is unfortunately misleading. It has some truth to it, but it is not really the correct answer from a group theoretic standpoint. In particular it is true and important to know that the Hodge star operator anti-commutes with parity inversion. But as Jagerber48 points out, the Hodge star argument cannot ever simultaneously account for the properties of both $\vec{E}$ and $\vec{B}$ under the action of $\mathcal{P}$ and $\mathcal{T}$.
The problem here is that when we begin to care about reflections, we are not merely looking at representations of the proper orthochronus Lorentz group $SO_+^{\uparrow}(1,3)$ as we are so often used to, but of the full Lorentz group $O(1,3)$. The entire Lorentz group has four distinct connected components of group elements, which can be reasonably named the identity component, the parity inverted component, the time reversed component, and the time-parity reversed component. $$ \mathcal{I}, \>\>\mathcal{P},\>\> \mathcal{T},\>\> \mathcal{PT}.$$
While these names are physically meaningful, the fact that all these operators are involutions, and that $\mathcal{P}$ and $\mathcal{T}$ are reflecting exactly opposite compact subspaces means we can equally understand these labels as $$\mathcal{I},\>\> \mathcal{P},\>\> -\mathcal{P},\>\> -\mathcal{I}$$
The elements $\mathcal{I}$ and $-\mathcal{I} = \mathcal{PT}$, being proportional to the identity, will commute with every element of the Lorentz group; the set of those group elements with this property is called the center of the Lorentz group, often denoted $\mathcal{Z}\left[O(1,3)\right] = Z_2$. I.e. the center here is a discrete subgroup of two elements.
It seems obvious to state that the two elements of the center act differently when applied to $4$-vectors. A more surprising fact is that acting upon $F^{\mu\nu}$ with the $\mathcal{PT}$ from our vector representation, is the same as doing nothing at all! (Try it yourself)
This is because the second tensor power of the vector representation of $O(1,3)$ does not inherit a faithful representation of the full Lorentz group, but instead forms a representation of $PO(1,3)$, known as the projective orthogonal group (this is exactly analogous to how a spinor is nontrivially affected under a rotation by $2\pi$, yet vectors are left unaffected, $O$ doubly covers $PO$ in even dimensions). The Projective orthogonal groups are precisely defined as the 'centerless' versions of the orthgonal groups. What this means abstractly is that in $PO(1,3)$ any Lorentz transforms which are related by application of an element of the center, are identified. So out of the four connected components of $O(1,3)$, both $\mathcal{I}$ and $\mathcal{PT}$ become $\mathcal{I}$ in $PO(1,3)$, and both $\mathcal{P}$ and $\mathcal{T} (= -\mathcal{P})$ become $\mathcal{P}$ in $PO(1,3)$. This is precisely what your argument shows, and this is how we know we are in this representation of the Projective Orthogonal group: if an object transforms under $PO(1,3)$, there are only two kinds of Lorentz transforms for it, not four.
But this is bad, because we know $\vec{E}$ and $\vec{B}$ fields have distinct time reversal properties. The reason for the 'error' in our ability to perform this operation is in the assumption that when our representation of the proper orthochronus Lorentz group (the component connected to the identity) is faithful, that so to must be our representation of $O(1,3)$.
There are two ways out of this problem: the more cumbersome one, is the option of realizing the adjoint representation (where $F$ lives) as a subgroup of the vector representation of $O(3,3) \supset O(1,3)$. This explicitly $6$ dimensional representation of the Lorentz group has a non-trivial center, as we desire.
The second option, is to perform a discrete central extension on our 'accidentally' projective Lorentz transforms. I.e. for every Lorentz transform we can apply to $F$: $$\Lambda^T F \Lambda$$
we introduce the possibility of being transformed instead by $$-\left(\Lambda^T F \Lambda\right).$$ Thus when acting on rank two tensors, transforms from $\mathcal{I}$ and $\mathcal{P}$ are given as usual, and the transforms from $\mathcal{T}$ and $\mathcal{PT}$ acquire an additional minus sign. With this upgraded transformation property the second tensor power becomes a faithful representation of the full Lorentz group! Where before we found we could only end up transforming $F$ by elements of the two components $\mathcal{I}$ and $\mathcal{P}$, by the inclusion of the minus sign we have found our way back to Lorentz transforms in the other two discrete connected components $\mathcal{T}=-\mathcal{P}$ and $\mathcal{PT} = -\mathcal{I}$
Thus the answer to your question is that abstractly, for any particular real-representation $r$ we know that, $$\mathcal{T}_r=-\mathcal{P}_r.$$ However under the tensor product of vector representations we have seen $$\mathcal{P}_v \otimes \mathcal{P}_v = \mathcal{P}_F$$ but of course then $$\mathcal{T}_v \otimes \mathcal{T}_v = -\mathcal{P}_v \otimes -\mathcal{P}_v = \mathcal{P}_F$$ and so we find we pass to a representation of the Projective orthogonal group. Then, through a central extension of tagging $\mathcal{T}$ and $\mathcal{PT}$ transformations with an overall negative sign, we find ourselves once again in a faithful representation of the full Lorentz group $O(1,3)$ and regain the ability to apply the time reversal operator:
$$\mathcal{T}_F = -\mathcal{P}_F$$
which when acting on $F$ appropriately reverses the magnetic field and leaves invariant the electric field, as expected.