ClearAll
dataN=234
datafather=List[
165,169,173,174,175,174,175,171,175,178,163,170,175,172,174,173,175,173,169,180,172,161,168,172,178,165,163,160,177,169,175,170,167,168,170,163,170,168,171,174,171,170,164,175,175,171,155,175,171,165,175,175,168,160,165,170,166,164,161,175,175,160,163,170,175,167,170,172,172,170,174,170,168,172,169,175,169,178,189,175,164,175,155,175,173,165,160,167,175,176,178,180,167,168,165,168,167,164,170,155,174,160,174,170,171,173,170,168,165,175,175,175,165,170,160,170,179,176,176,168,175,165,166,170,180,170,173,165,173,169,174,176,170,180,173,178,170,170,164,176,175,173,178,170,175,170,170,175,169,187,160,173,175,165,170,159,171,174,175,170,173,182,168,170,167,170,175,178,165,172,172,165,169,164,171,168,159,171,160,183,173,167,178,171,160,173,175,170,170,175,172,181,167,179,165,182,167,183,170,167,170,168,170,164,178,166,168,175,165,172,168,166,173,168,172,175,169,169,165,171,171,172,176,180,165,183,165,167,179,178,165,176,180,170,171]


N[Mean[datafather]]
N[Variance[datafather]]
Sqrt[%]

datamother=List[
156,163,159,161,158,150,155,163,158,148,163,155,167,161,161,160,165,152,152,165,150,157,150,162,155,154,158,150,156,150,155,158,158,160,155,155,165,165,164,157,164,150,160,155,158,170,160,150,165,153,156,163,160,165,150,164,153,166,153,155,152,155,154,154,165,153,165,147,168,160,170,161,158,165,162,165,162,159,170,165,166,160,145,164,163,160,150,153,154,150,163,165,155,160,159,150,165,150,164,165,160,165,153,160,158,156,150,154,155,163,150,150,150,170,162,160,170,165,160,148,160,154,160,155,147,146,160,165,162,157,162,156,160,165,154,158,165,160,160,160,160,155,164,165,160,160,150,165,158,164,154,156,163,155,162,160,156,167,160,160,164,162,150,155,156,169,158,156,158,174,155,155,161,165,163,165,157,156,150,160,163,164,159,155,167,158,160,163,150,160,155,158,158,165,155,155,166,162,165,168,157,162,150,152,155,165,155,158,153,163,156,157,163,157,155,162,155,155,155,168,159,155,155,165,155,165,150,145,162,158,155,158,152,160,156]

N[Mean[datamother]]
N[Variance[datamother]]
Sqrt[%]

databoy=List[
168,166,174,173,172,171,160,169,170,165,161,165,171,170,176,175,175,163,167,174,162,171,170,174,172,173,172,170,172,169,168,172,171,175,164,158,178,179,177,168,178,170,180,170,173,172,177,178,160,172,174,179,173,170,167,167,163,173,165,168,168,165,173,167,173,160,169,169,160,171,185,169,168,164,172,171,180,175,179,169,172,178,159,164,171,170,166,172,171,174,173,175,170,163,160,165,176,160,170,152,175,180,170,163,172,178,165,172,176,175,170,178,166,175,171,171,182,174,180,169,175,169,161,166,160,170,175,170,186,173,178,168,170,170,172,172,177,168,171,172,170,160,186,173,174,172,170,176,165,180,166,180,165,170,172,167,166,177,162,163,173,170,166,176,171,170,180,173,161,170,170,164,172,180,181,170,165,179,165,168,170,173,170,162,168,177,178,173,165,165,175,172,166,178,168,183,177,187,170,170,167,170,170,160,168,171,173,171,159,175,164,178,174,174,176,168,168,159,182,170,169,174,173,184,165,169,161,157,178,176,172,164,169,177,173]

N[Mean[databoy]]
N[Variance[databoy]]
Sqrt[%]

datagirl=List[
154,155,159,158,162,153,155,164,155,164,152,150,161,164,158,162,165,162,153,165,153,154,153,160,166,163,160,150,156,160,160,162,162,162,152,155,170,168,170,157,159,160,163,157,175,155,154,170,158,154,164,167,166,165,155,165,153,165,149,166,157,160,160,152,160,136,156,145,153,157,165,155,158,157,165,162,163,160,169,165,163,163,150,158,165,165,153,165,150,168,165,165,158,161,155,158,160,160,167,160,162,155,156,150,153,152,150,154,160,159,159,170,160,152,152,185,165,167,160,160,158,148,152,152,155,155,170,150,163,160,160,158,158,170,156,156,160,157,161,165,163,155,166,160,161,165,160,165,155,171,158,156,140,155,160,164,154,158,164,153,165,161,150,160,157,167,160,163,153,157,160,150,162,160,153,160,160,160,145,164,152,160,160,156,153,160,160,165,155,155,160,170,160,166,160,155,155,160,160,168,160,162,160,156,159,165,149,164,156,160,158,155,163,170,168,168,155,148,155,163,162,164,160,180,150,163,150,149,170,164,157,164,166,160,153]
N[Mean[datagirl]]
N[Variance[datagirl]]
Sqrt[%]

datafm=Table[{datafather[[i]],datamother[[i]]},{i,1,dataN}]
datafb=Table[{datafather[[i]],databoy[[i]]},{i,1,dataN}]
datafg=Table[{datafather[[i]],datagirl[[i]]},{i,1,dataN}]
datamb=Table[{datamother[[i]],databoy[[i]]},{i,1,dataN}]
datamg=Table[{datamother[[i]],datagirl[[i]]},{i,1,dataN}]
databg=Table[{databoy[[i]],datagirl[[i]]},{i,1,dataN}]

N[Correlation[datafm]]

N[Correlation[datafb]]

N[Correlation[datafg]]

N[Correlation[datamb]]

N[Correlation[datamg]]

N[Correlation[databg]]

MenMin=150;MenMax=185;WomenMin=140;WomenMax=180;
modelfm=LinearModelFit[datafm,x,x]
Show[ListPlot[datafm, PlotRange->{{MenMin,MenMax},{WomenMin,WomenMax}}, AxesOrigin->{MenMin,WomenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]]]

modelfb=LinearModelFit[datafb,x,x]
Show[ListPlot[datafb, PlotRange->{{MenMin,MenMax},{MenMin,MenMax}}, AxesOrigin->{MenMin,MenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]]]
Show[ListPlot[datafb, PlotRange->{{MenMin,MenMax},{MenMin,MenMax}}, AxesOrigin->{MenMin,MenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]],Plot[modelfb["BestFit"],{x,0,200}]]

modelmg=LinearModelFit[datamg,x,x]
Show[ListPlot[datamg, PlotRange->{{WomenMin,WomenMax},{WomenMin,WomenMax}}, AxesOrigin->{WomenMin,WomenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]]]
Show[ListPlot[datamg, PlotRange->{{WomenMin,WomenMax},{WomenMin,WomenMax}}, AxesOrigin->{WomenMin,WomenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]],Plot[modelmg["BestFit"],{x,0,200}]]

modelfg=LinearModelFit[datafg,x,x]
Show[ListPlot[datafg, PlotRange->{{MenMin,MenMax},{WomenMin,WomenMax}}, AxesOrigin->{MenMin,WomenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]]]
Show[ListPlot[datafg, PlotRange->{{MenMin,MenMax},{WomenMin,WomenMax}}, AxesOrigin->{MenMin,WomenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]],Plot[modelfg["BestFit"],{x,0,200}]]

modelmb=LinearModelFit[datamb,x,x]
Show[ListPlot[datamb, PlotRange->{{WomenMin,WomenMax},{MenMin,MenMax}}, AxesOrigin->{WomenMin,MenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]]]
Show[ListPlot[datamb, PlotRange->{{WomenMin,WomenMax},{MenMin,MenMax}}, AxesOrigin->{WomenMin,MenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]],Plot[modelmb["BestFit"],{x,0,200}]]

modelbg=LinearModelFit[databg,x,x]
Show[ListPlot[databg, PlotRange->{{MenMin,MenMax},{WomenMin,WomenMax}}, AxesOrigin->{MenMin,WomenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]]]
Show[ListPlot[databg, PlotRange->{{MenMin,MenMax},{WomenMin,WomenMax}}, AxesOrigin->{MenMin,WomenMin}, AspectRatio -> 1,PlotStyle->PointSize[Large]],Plot[modelbg["BestFit"],{x,0,200}]]

datafmb=Table[{datafather[[i]],datamother[[i]],databoy[[i]]},{i,1,dataN}]

datafmg=Table[{datafather[[i]],datamother[[i]],datagirl[[i]]},{i,1,dataN}]

modelfmb=LinearModelFit[datafmb,{x,y},{x,y}]

modelfmg=LinearModelFit[datafmg,{x,y},{x,y}]

plot1=ListPointPlot3D[datafmb, PlotRange->{{MenMin,MenMax},{WomenMin,WomenMax},{MenMin,MenMax}}, AspectRatio -> 1,PlotStyle->PointSize[Large],Filling->Bottom]
plot2=Plot3D[modelfmb["BestFit"],{x,MenMin,MenMax},{y,WomenMin,WomenMax}, MeshShading ->{{None,None},{None,None}}, Mesh -> 10]
Show[plot1,plot2]