I'm trying to fit this [data](https://drive.google.com/file/d/1D5k1xaTE6ferPpfOomhOEBio8oEIAOvF/view?usp=sharing)

[![data](https://i.stack.imgur.com/hujdb.jpg)](https://i.stack.imgur.com/hujdb.jpg)

As you can see there's two peaks. I want to fit this curve with two normal distribution with the same variance but different mean.

So far I've got this:

```
nlm = NonlinearModelFit[data, {PDF[MixtureDistribution[{w, 1 - w}, {NormalDistribution[m1, o], NormalDistribution[m2, o]}], x], {0 < w < 1, m1 > m2 + 1, m2 > 5, 3 > o > 0}}, {m1, m2, w, o}, x]
Show[Plot[Normal[nlm], {x, 0, 40}], ListPlot[data], PlotRange -> All]

```

[![fit](https://i.stack.imgur.com/h4cMV.jpg)](https://i.stack.imgur.com/h4cMV.jpg)

You can see that the fitting seems about right for the means, but the areas are wrong. Integration of the area gives 1 and the sum of the points also. There's seem to be I'm missing something, probably in relation to the PDF