Wien's Displacement Law for real bodies It is known that for perfect blackbodies,
$$\lambda T= c$$
where $\lambda= \text{peak wavelength}$
      $T= \text{Absolute temperature}$
      $c= \text{Wien's constant}$                                                   
But this is for perfect blackbodies only, which have no theoretical existence. Does a similar formula exist for real bodies, which expresses $\lambda T$ in terms of its emissitivity $\epsilon$? I googled it, but found no relevant results.       
 A: 
But this is for perfect blackbodies only, which have no theoretical existence.

Just for clarification, the Wien displacement law is not valid for emission spectrum "of blackbodies only". It is valid for any radiation in thermodynamic equilibrium with matter (equilibrium radiation), which is a more general concept that does not require presence of any black body. This equilibrium radiation has energy density spectrum given by the Planck function
$$
\rho_{Planck}(\lambda, T) = \frac{8\pi h c}{\lambda^5}\frac{1}{e^{\frac{hc}{\lambda kT}}-1}
$$
and can be produced and maintained by non-black bodies, for example by enclosure of any material brought to thermodynamic equilibrium at temperature $T$. Wien's displacement law says that the product of $T$ and $\lambda_{max}$ for which the energy density spectrum is maximum is the same for all temperatures.

Does a similar formula exist for real bodies, which expresses λT in terms of its emissitivity ϵ?

Let me rephrase your question : for emission intensity spectrum $I(\lambda, T)$ of non-black body with known emissivity $\epsilon(\lambda, T)$, is there a simple relation between the wavelength of maximum and temperature, similar to above Wien's displacement law ?
The answer is that in general, most probably no, because the emission radiation of real body has intensity spectrum proportional to
$$
I(\lambda, T) \propto \rho_{Planck}(\lambda, T) \epsilon(\lambda, T),
$$
which means that the number, positions and height of peaks are modified by the emissivity $\epsilon(\lambda, T)$. This depends on the body, its material and thickness and can be fairly general function of $\lambda$ and $T$. There may be some general restriction on this function due to physical laws, but Planck's formula does not imply any, apart from the condition that $0<\epsilon \leq 1$ for all frequencies. One would have to study more detailed model of matter explaining value of $\epsilon(\lambda, T)$ to find such relations.
A: See page 4 of Smith's Computing Plank Function, in the section titled "5. Non-Black Body radiation." My answer is merely a simplified, paraphrased version of what Ronald B. Smith wrote!
You'll have to use Planck's Law rather than Wien's Law (Wien's law can be derived from Planck's Law):
$$B_{\lambda}(T)=\frac{2hc^2\lambda^{-5}}{e^{hc/k\lambda T}-1}$$
where T is temperature, B is radiance, lambda is wavelength, and the others are constants.
$$I=\epsilon B$$
where I is intensity, B is radiance, and epsilon is emissivity. The emissivity of a 60 W incandescent lightbulb is 0.4 according to this study.
Replacing $B$ with $I/\epsilon$, and rearranging for T, we get:
$$T=\frac{hc}{k\lambda\ln(\frac{2hc^2\epsilon*\lambda^{-5}}{I}+1)}$$
Hope this helps :D
