Calculating total number of infected from R0

The first piece of code solves for the final number of susceptibles given a value of R0.

Solve[1==SFraction-1/R0*Log[SFraction],SFraction]

f[R0_]:=-((ProductLog[-(E^(-R0)) R0])/R0);

I am now going to plot this function to indicate how much of the population is predicted to be infected.

plot1=Plot[1-f[R0],{R0,0,20},PlotRange->{-.01,1.1},PlotStyle->{Black,Thick},AxesLabel->{"(Subscript[R, 0])","Total Fraction Infected"},AxesStyle->Directive[Black, 20]]

I am now going to put on this graph the values of (Subscript[R, 0])for various diseases.

DataPoints={{2,1-f[2]},{1.5,1-f[1.5]},{17,1-f[17]},{12,1-f[12]},{5,1-f[5]},{3,1-f[3]}};
plot2=ListPlot[DataPoints->{"Covid","H1N1","Measels","Chicken Pox","Polio","Smallpox"},PlotStyle->{Red,PointSize[Large]},LabelStyle->{16}];
Show[plot1,plot2]